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ABSTRACT 



Lower confidence limit estimation procedures for the reliability of several 
systems are developed and their accuracies evaluated using computer simulation. 
The procedures use test data on components of the system which can have failure 
times with either exponential or Weibull distributions or both. Testing scenarios for 
the components can be truncated by number of failures or by planned test times. 

Although the evaluation effort was focussed on series systems in this thesis, the 
procedures readily apply to other systems as described in the thesis. The evaluations 
demonstrate the procedures to be quite accurate when sufficient component testing 
is performed. 

Two FORTRAN computer programs were written to perform the evaluation. 
They are annotated in Users’ Guides and can be used to determine the accuracy of 
these approximate lower confidence limit procedures for a given specific system and 
associated set of input parameters. 



6 .! 

THESIS DISCLAIMER 

The reader is cautioned that the computer programs developed in this 
research may not have been exercised for all cases of interest. While every effort has 
been made, within the time available, to ensure that the programs are free of 
computational and logic errors, they cannot be considered validated. Any application 
of these programs without additional verification is at the risk of the user. 
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I. INTRODUCTION 



This thesis develops approximate lower confidence interval procedures for the 
reliability of complex systems using test data on components of the system. The 
accuracies of these procedures are also assessed using computer simulation. The 
procedures can be used for any complex system whose reliability does not decrease 
when the reliability of any one of the components is increased. 

The failure times of the continuously operating components are assumed to 
have either an exponential or Weibull distribution. Parameters in both distributions 
are assumed to be unknown. The Weibull distribution is used to model the lifetime 
probability distribution of electronic components with non-constant failure rate 
functions. It is also used to model the lifetime probability distributions of mechanical 
devices, since their failure rate functions are usually increasing with operating time. 

Lower confidence limit estimation procedures for system reliability are needed 
during the development phase of systems to provide indications of a contractor’s 
ability to meet a stated system reliability goal as development progresses and the 
results of test programs become available. These procedures are also needed to 
assess the reliability of systems that have been operating in the field for some time 
and have accumulated histories of failure data and unique configurations of modified 
or repaired components. 

Few textbooks on reliability treat the problem of system reliability interval 
estimation. Those that do usually limit the discussion to series or parallel systems. 
Moreover, the procedures they present are not adaptable to other more complex 
systems. Mann, Shafer, and Singpurwalla [Ref.l pp 487-524] provide one of the better 
treatments of a variety of these methods in Chapter 10 of their book. This chapter 
provides an excellent summative discussion of the many procedures that were 
developed from 1954 to 1974. However, none of the procedures reviewed in their 
book can accomodate the use of test data from a mix of components with both 
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exponential and Weibull failure time distributions. The procedures presented in this 
thesis does accomodate this type of system with a mixture of different component 
types. In addition, the procedures presented in this thesis can accomodate an 
additional mix of components for which only attribute data has been collected. 

Procedures developed in this thesis are extensions of a procedure developed by 
Myhre, Sanders and Rosenfeld [Ref.2], In their paper, they assume the failure times 
of continuously operating components have exponential distributions with associated 
failure rates, X { . The test data on the remaining components, the number of observed 
failures/;- in n { tests, are assumed to have Poisson distributions with associated means 
n t q t . They assume the ratios X/Xj , q/qj and X/qj are known and develop confidence 
interval estimation procedures for system reliability that use this information. They 
also show that the accuracy of their procedure is not very sensitive to moderate 
inaccuracies of these ratios. This suggests that it might be possible to estimate the 
ratios from the data as part of the interval estimation process and not suffer 
significant loss of accuracy in the interval estimates. Estimating these ratios is part 
of the procedures developed in this thesis. 

This thesis also provides an annotated computer program that can be used to 
assess the accuracies of the lower confidence limit procedures when applied to any 
specific system. Sufficient annotations are provided throughout the program in 
Appendix C. This program provides the user with a means for verifying the accuracy 
of these proposed lower confidence limit procedures for his specific system and 
testing program, that is, sample sizes and type of truncation. This capability will 
allow the user to answer many "what if type of questions. 
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II. THEORY 



A. Interval Estimation Procedure for Exponential Failure Times 

A system is defined to be quasi-coherent if an increase in reliability of any one 
of its components does not cause a decrease in system reliability. The components 
of a quasi-coherent system do not need to be statistically independent. However, 
throughout this thesis, it is assumed that all components are statistically independent. 

Suppose a quasi-coherent system has k components and the distribution of the 
failure time of component i is exponential with failure rate A,-. Then the system 
reliability R s can be written as a function of A,- , j = 1, 2, ..., k as follows: 

*,(0 = g( *1*^2 f'dk ) 



where t t is the operating time for component i. Let m be any one of the k 
components and r ( - = , for i - 1, 2, ..., k. Then equation (2.1) may be viewed 

as 



R s ( 0 = §( x »,’ r v r 2’—’ r k’ V* •••«'* ) - 

If the r/s are known and ^ o) were an upper 100(l-a)% confidence limit for A m , 

the corresponding lower confidence limit for R s (t) would be: 

■^(OlO) = s( ^ m,U(ay r v r 2’—’ r k’ { V { 2 ) *" ^-3) 

Specifically, if we have a series system of independent components, so that 



R s (t) = exp{ -2 A tU ) 



k 

L 

i-l 



... (2.4) 



= ex P { } 

i-l 
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then, 



^j(OlO) eX P^ ^m,U(a)^l r i t i ^ 



... (2.5) 



i=i 



If n i items of component i are tested until/) failures occur, 7) denotes the total test 

k 

time accumulated on all the n, items, and F-Zf, then the expression 

i=i 



2K,Zr<T, 

/* 1 



has a Chi-square distribution with 2 F degrees of freedom. See Bain and Engelhardt 
[Ref.3], The corresponding 100(l-a)% upper confidence limit for A m is 






F 



/*! 



... ( 2 . 6 ) 



where X 2 «, 2 f i s th e 100(l-o:)th percentile point of a Chi-square distribution with 2 F 
degrees of freedom. 

If the testing on component i is terminated when a total test time of 7) has 
been accumulated by all n, items, then the equation for X ni t/(a ) becomes 



X c,2(l +F) 

— ... (2.7) 

22> i 7 ', 

/*! 



In this case /) is random and so is F. 

If testing on each of the n, items of component i are tested until a planned test 
time or failure, and failed items are replaced immediately, then equation (2.7) will 
be the exact expression for < k m . If failures are not replaced, then equation (2.7) 

is approximate. See Lee, Bain and Englehardt [Ref.3 pp 486-495]. Department of 
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Defense document NAVSEA OD29304B "Reliability and Availability Evaluation 
Program Manual" [Ref.4 p 5-42] provides nearly exact procedures for X m when 

testing is terminated by planned test time for each item tested and failures are not 
replaced. 

The values of the r/s are assumed to be unknown in this thesis. When testing 
is terminated by the number of failures, a nearly unbiased estimator for r f is 




where = (/■ - 1 ) / T i and the index m denotes the component with largest value 

of A ; . . The ratio, (f t ~l )/Ti , is an unbiased estimator for A, (see Appendix A). If 1 / 
were unbiased for 1 /A m then r. would be an unbiased estimator for r ( . Replacing X m 

with % m f m / (f m ~ 1 ) in equation (2.8) will yield an unbiased estimator r ( . for r { . 

Multiplying by this constant is nullified by a cancellation with the same 

constant in the final equation for the system reliability lower confidence limit, so 
equation (2.8) is used to estimate r,-. Using estimator r i for r, , equation (2.6) 

becomes 



* 0 , 2 F 

m ,u(a) — ... ( 2 . 9 ) 

/=! 



It is important to note that the index m denotes the component for which 

= (f r l )/T i is the largest among all the components in the system. The corresponding 
equation for the 100(l-a)% lower confidence limit on the reliability of a series 
system is 
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... ( 2 . 10 ) 



k 

4(0 l(b) = ex p { -IwEVi } 

i-i 

The corresponding lower confidence limit for the reliability of any quasi- 
coherent system is given by equation (2.3) with r f replaced by r t . 

In this thesis, equation (2.8) with = f i /T i will also be used to estimate r { 

under exponential assumptions when testing is terminated after an accumulated test 
time is achieved (truncated), and when at least two components have at least one 
failed test item. This is done because failures will not be replaced in the time 
truncated test plans that are simulated in this thesis. In this type of testing both /• 
and T, are random. Under time truncation, it is possible that no failures will occur 
on any component tested in which case equation (2.8) is undefined. Also, if only one 
component has one failure and the remaining components have zero failure, equation 
(2.8) would be zero for all i except the case when i = m. 

All of the confidence limit procedures in this thesis have a common special 
method for computing the lower confidence limit of system reliability in the two 
cases of zero or one failure. This feature amounts to a modification to equations 
(2.3) and (2.10). 

When either zero failures or one failure have occured among all components, 
the test data is examined for each component to determine the total number, N { , of 
equivalent component mission tests (/ = 1, 2, ..., k). These N- s and the system 
configuration are analyzed to determine the equivalent number of mission tests, N, 
for the system that would have occurred if Nj , N 2 , ..., N k of these k components 
were assembled into systems. 

For a series system, this N will be equivalent to min { N lt N 2 , ..., A* }■ The 
100(l-a)% lower confidence limit of system reliability if zero failures occurred is then 
computed directly as follows: 
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- % - < 211 > 

If exactly one failure occurred among all k components, w hl t> e ^e 



solution for p in the equation 



p N + Np N ~* (l-p) = a 



... ( 2 . 12 ) 



These confidence limit equations are the standard binomial lower confidence limit 
equations. Equations (2.11) and (2.12) are part of the set of equations used to 
compute ^(0 l( O ) f° r a U °f the time truncated interval estimation procedures in 



It is important to remember that the symbol 7) in equations (2.7), (2.8) and 
(2.9) denote total accumulated test time for component i ; that is 



where is the test time accumulated on the ;th test item of component i and n, is 
the number of test items of component i being tested. 

B. Interval Estimation Procedure for Weibull Failure Times 

Consider a series system with k components. Let the time to failure, X t , of 
component / have a Weibull distribution with density 



this thesis. 



T, - £ T a 



... (2.13) 




Then 



K/(0 = exp{ -(V,) P ‘ } > l i > 0 



... (2.15) 
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and 



R s(0 = n exp{-A,V‘} 

«-l 

= exp{-^A,V} ...(2.16) 

i= 1 

= exp{-A;^r//'} 
i-1 

where A ; ‘ = A ; Pl , A m * is any one of the \ ’ , i = 1, 2, k , and r, = \ If the 

ySj’s are known, then Xf l will have a constant failure rate A, Pl and the procedures 
described in Section A can be used to obtain j$ 5 (0l(«) with 7E replaced by T in 

equation (2.13). 

Suppose jSj is unknown and X i(I) , X i>(2) , ..., X i(fi) are the ordered failure times 
under either type of truncated testing for component i in the system. Solutions 
and for fa and \ in the two equations given in equation (2.17) are the 

maximum likelihood estimates for fa and See Mann and others [Ref.l pp 189-191]. 
These equations are used for both types of test truncation. If for component i, 
testing is terminated on the fl h failure, then t u = X i(fi) in equation (2.17). The 
solution, , is a biased estimator for fa. Bain [Ref.5 pp 220] provides a table of 

constants B(nJ which depends on number of test items such that /),* = faB(n ( ) 



is a nearly unbiased estimator for fa. 

f, 



E (0 . (".-/iXiV* 

j - 1 



ft 



f, 



E */« . 

1=1 



P, 



1 = ir \nX i(n 

/>,- fib 0) 



... (2.17a) 
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and 




/, 

E-xiS, * 



... (2.17b) 



If the testing for component i is terminated at failure or at a given time t oi for 
each of the n, items on test, then t ^ = t oi in equations (2.17a) and (2.17b). 

Now, let 






with i = 1, 2, ..., k 
j = 1, 2, ..., n i 



... (2.18) 



In this thesis, the distribution of 7^ is approximated by the exponential distribution 

P * 

with failure rate A, ' s and procedures similar to those in Section A are used 



to obtain the lower confidence limit on system reliability. Define 




* 

where 7\ =Y, T ij > 1 = 2 > •••> k . Let X' m = max^,. £*. 

Note that this defines the index m. Define r ( . by 



... (2.19) 



X* T 

- 7f - *;•<£> 

•/ m 



x: 



... ( 2 . 20 ) 



for both types of test truncation plans. 

Then an approximate 100(l-a)% upper confidence limit for X^, is given by 
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m,U{ o) 




... ( 2 . 21 ) 



where 



F" 




‘ ’ for all components 



for all components 



if test till specified time 



if test till f th failure 



... ( 2 . 22 ) 



The corresponding approximate 100(l-a)% lower confidence limit R s (t) L(a) for the 
reliability R s (t) of a series system is given by 



when at least two components have at least one failure. Equations (2.11) and (2.12) 
also apply here when the total failures over all components is either zero or one. 

The accuracies of these approximate confidence interval procedures were 
evaluated by using computer simulations which are described in the next chapter. 
During this evaluation process, the degrees of freedom in the expressions X 2 b,2F* 

in equation (2.21), X 2 ^ 2 f in equation (2.9) and X 2 0 , 2 (i *f) in equation (2.7) were 

increased and decreased from the defined values of F* and F given by these 
equations. The purpose of these modifications was to find more accurate lower 
confidence limit procedures. The specific increases and decreases are described in 
Chapter III. The results show that for some cases the procedures with modified 
degrees of freedom are more accurate. 




... (2.23) 
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III. COMPUTER SIMULATION 



A. Test Plan 1 : Testing w, Until / Failures (RETP1) 

RETP1 is a program written in FORTRAN, on the Amdahl mainframe 
computer, which performs the computer simulation of the random failure times of 
the different types of components in the system. A documentation of this program 
and its associated subroutines is included in Appendix B. 

The program accepts input parameters via an input file IN1.DAT. For each 
replication, it generates the failure times for all the component items included in the 
test plan using a uniform random number generating subroutine LRNDPC. A quick 
evaluation of LRNDPC (see Appendix D) by plotting U(n + 1) vs U(n) illustrates the 
uniformity of the routine. The program determines the total test time accumulated 
for each component in the system and computes the estimates of the key parameters 
and the consequent lower confidence limit for system reliability for that replication. 
The process is repeated 1000 times. When all replications are done, the routine 
EVAL processes the lower confidence limit estimates from all 1000 replications and 
determines the two measures of accuracy for the run, namely RSLOW and LEVEL. 

RSLOW is the 100(l-a) percentile of the ordered set of lower confidence limits 
from the 1000 replications computed in a run. The true reliability of the system is 
RS. The closer RSLOW is to RS, the greater the accuracy of the procedure under 
evaluation in the run. If the procedure is exact, RSLOW will be equivalent to RS. 
To be conservative, RSLOW should always be lower than RS. 

LEVEL measures the proportion of 1000 lower confidence limits, from a run 
with 1000 replications, which are lower than the true system reliability RS. The 
closer LEVEL is to the specified confidence level for the procedure, 1-a, the better 
the procedure. Values of LEVEL greater than (1-a) reflect an under-estimation of 
RS which is conservative. Values of LEVEL lesser than 1-a signal an over- 
estimation of RS which may be undesirable. 
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Simulation runs are performed using RETP1 for all combinations of failure 
time distributions and levels of key input parameters listed below. 



(a) System. - 8 Exponential components in Series (Case 1) 

- 8 Weibull components in Series (Case 2) 

- 4 Exp and 4 Wei (Mixed) components in Series (Case 3) 

(b) True System Reliability (RS). 

- Hi (greater than 0.9) (Type A) 

- Lo (greater than 0.8) (Type B) 

(c) Level of Significance (a). 



- 0.1 
- 0.2 



(d) Degrees of Freedom for x 2 statistic (DF) as a function of the total number of 
failed test components (NFC) and total number of system components (NCOMP). 

- DF = 2 * NFC 

- DF = 2 * ( NFC + NCOMP ) 

- DF = 2 * ( NFC - NCOMP ) 

- DF = 2 * NFC - NCOMP 

e) Test Plan for each component. 

- Test 5 until 5 failures 

- Test 15 until 15 failures 

- Test 15 until 11 failures 

- Test 15 until 7 failures 

- Test 15 until 3 failures 



For the 8 exponential components in Case 1, the mission time for each of the 
component is chosen to be 10 hrs. The program will accomodate different 
component mission times. The chosen values of the scale parameters, \ , were 
different depending on whether the system is highly reliable (Type A) or one with 
a lower reliability (Type B). The ratios between the largest and the smallest failure 
rate was chosen to be 8 and 4.5 respectively for Type A and Type B systems. 
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For the 8 Weibull components in Case 2, the mission time for each of the 
components was chosen to be 10 hrs. The chosen values of the scale parameters, \ iy 
were different depending on whether the system is highly reliable (Type A) or one 
with a lower reliability (Type B). The ratio between the largest and smallest failure 
rate was chosen to be 8 for both system types. The shape parameter is chosen to be 
2 for all cases. The program will accomodate any value greater than zero for the 
shape parameter. 

A mixture of exponential and Weibull components with those parameters 
described in the last two paragraphs is chosen for the Type A and Type B systems 
of Case 3. 

Each simulation run of 1000 replication results in an output file OUTl.DAT. 
The raw output from all the RETP1 runs are summarized in tabular form and placed 
in Appendix E. Each table corresponds to a specific run case and system type 
combination. 
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B. Test Plan 2 : Testing for a Specified Planned Test Time (RETP2) 

RETP2 is another program written in FORTRAN, on the Amdahl mainframe 
computer, which performs the computer simulation of the random failure times of 
the different types of components in the system. A documentation of this program 
and its associated subroutines is included in Appendix C. 

The structure of this program is quite similar to that of RETP1 described in 
Section A of this chapter. The program accepts input parameters via an input file 
IN2.DAT. For each replication, it generates the failure times for all the component 
items included in the test plan using LRNDPC. The program then determines the 
number of failed test components for each component in the system and computes 
the estimates of the key parameters and the consequent lower confidence limit for 
system reliability for that replication. The process is repeated 1000 times. When all 
replications are done, the routine EVAL processes the lower confidence limit 
estimates from all the 1000 replications and determines the two measures of accuracy 
for the run, namely RSLOW and LEVEL. The definitions of these two measures 
were discussed in the Section A. 

Simulation runs are performed using RETP2 for all combinations of failure 
time distributions and levels of key input parameters listed below. 

(a) System. - 8 Exponential components in Series (Case 4) 

- 8 Weibull components in Series (Case 5) 

- 4 Exp and 4 Wei (Mixed) components in Series (Case 6) 

(b) True System Reliability (RS). 

- Hi (greater than 0.9) (Type A) 

- Lo (greater than 0.8) (Type B) 

(c) Level of Significance (a). 

- 0.1 
- 0.2 
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(d) Degrees of Freedom for x 2 statistic (DF) as a function of the total number of 
failed test components (NFC) and total number of system components (NCOMP). 

- DF = 2 * ( 1 + NFC ) 

- DF = 1.3 * 2 * ( 1 + NFC ) 

(e) 10 values of K = 0.25, 0.5, 1, 2, 3, 4, 5, 10, 20, 30 

where AT is a factor such that the expected number of failures for an exponential 
component during the specified total test time for that component is 0.6 times K. 
This K factor determines the bounds of the expected total number of failed test items, 
E[NFC]. The accuracy of the lower confidence limit procedures are highly correlated 
with E[NFC]. 

For the 8 exponential components in Case 4, the mission time for each of the 
component is chosen to be 5 hrs. The program can accomodate different component 
mission times. The chosen values of the scale parameters, were different 
depending on whether the resultant system is highly reliable (Type A) or one with 
a lower reliability (Type B). The failure rate was chosen to be 0.001 and 0.005 
failures/hr for all the components respectively. Total test time to be accumulated 
by each component is computed according to the following method. For each 
component i , r, represents the amount of operating time required to result in a 40% 
survival probability, that is, an expected failure of 0.6 component with an exponential 
failure time distribution. The computation for r, is as follows: 

R.ih) = exp(-X/ J .) 

= 0.4 

= -(|)ln(0.4) 

T t , the total amount of test time to be accumulated for component i would be 
K times r, which will result in an expected number of 0.6 times K failed items for this 
component. E[NFC] will then be 8 times that number, since there are 8 such 
components in the system. 
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For the 8 Weibull components in Case 5, the mission time for each of the 
component is chosen to be 15 hrs. The chosen values of the scale parameters, \ , 
were 0.005 failures/hr for the Type A system and 0.01 failures/hr for a Type B 
system. The shape parameter is chosen to be 2. For each Weibull component, a 
maximum of 20 test items are tested to failure. Estimation of E[NFC] and thus the 
total test time to be accumulated for each component is based on the exponential 
failure time model described in the earlier paragraphs. 

A mixture of exponential and Weibull components with those parameters 
described in the last two paragraphs is chosen for the Type A and Type B systems 
of Case 6. 

Each simulation run of 1000 replication results in an output file OUT2.DAT. 
The raw output from all the RETP2 runs are summarized in tabular form and placed 
in Appendix F. Each table corresponds to a specific run case and system type 
combination. 
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IV. RESULTS AND DISCUSSION 

The results of the simulation runs are summarized and discussed in this section. 
Tables 1A, IB, 2A, 2B, 3A and 3B in Appendix E, and Tables 4A, 4B, 5A, 5B, 6A 
and 6B in Appendix F present the accuracy results in tabular form for all run cases 
that were simulated: a few of these tables appear in this section to facilitate 
discussion of the results. 

A. Test Plan 1 : Testing n, Until Failures (RETP1) 

Table 1A displays the simulation results for Case 1 Type A. In this case, the 
system is comprised of 8 components in series. The failure time of each component 
has an exponential distribution. The failure rates of the 8 components range from 
0.0002 failures/hour to 0.0016 failures/hour. The mission time of the system is 10 
hours and the mission operating time of each component was also set to 10 hours. 
The component mission times do not need to be equal to the system mission time 
for the procedures evaluated in this thesis. This was discussed in Chapter II and is 
allowed for in all of the lower confidence limit equations. System reliability, RS, in 
Table 1A is 0.931 . Throughout this case, all of these parameters remain fixed. 

In simulation number 1 (S/N: 1) in Table 1A, five items for each of the 8 
components in the system are tested until they fail. Thus the number of failed 
components (NFC) is 40. One set of this 40 failure times is randomly generated for 
each simulation run. For this set of data, four 90% lower confidence limits and four 
80% lower confidence limits are computed. The four limits correspond to the values 
assigned to the degrees of freedom parameter F in the symbol X 2 „,2f which is a factor 
in the upper confidence limit equation for X U( ^ . The four different methods for 

computing this degree of freedom appear in the ”Deg of Freedom" column. NCOMP 
denotes the number of components in the system. Thus for each simulation run, 
eight lower confidence limits are computed. After 1000 replications of this 
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Table 1A : 



8 Exp in Series, RS = 0.931 (Hi) 

min \ = 0.0002 f/hr, max \ = 0.0016 f/hr, UT = 10 hrs 



S/N 


Test 

Plan 


Deg of 
Freedom 


a 


Measures of Accuracy 


RSLOW 


LEVEL 


1 


Test 5 until 
5 failed. 

NFC=40 


2*NFC 

(80) 


0.1 


0.919 


0.982 


0.2 


0.919 


0.960 


2*(NFC+ 

NCOMP) 

(%) 


0.1 


0.906 


1.000 


0.2 


0.906 


0.999 


2*NFC- 

NCOMP 

(72) 


0.1 


0.927 


0.949 


0.2 


0.929 


0.880 


2*(NFC- 

NCOMP) 

(64) 


0.1 


0.934 


0.8*1 


0.2 


0.933 


0.702 


2 


Test 15 until 
15 failed. 

NFC = 120 


2*NFC 

(240) 


0.1 


0.928 


0.955 


0.2 


0.927 


0.908 


2*(NFC + 
NCOMP) 
(256) 


0.1 


0.923 


0.990 


0.2 


0.923 


0.975 


2*NFC- 

NCOMP 

(232) 


0.1 


0.930 


0.916 


0.2 


0.932 


0.833 


2*(NFC- 

NCOMP) 

(224) 


0.1 


0.932 


0.844 


0.2 


0.932 


0.747 


3 


Test 15 until 
11 failed. 

NFC = 88 


2*NFC 

(176) 


0.1 


0.927 


0.955 


0.2 


0.926 


0.916 


2*(NFC + 
NCOMP) 
(192) 


0.1 


0.921 


0.996 


0.2 


0.920 


0.988 


2*NFC- 

NCOMP 

(168) 


0.1 


0.930 


0.916 


0.2 


0.929 


0.843 


2*(NFC- 

NCOMP) 

(160) 


0.1 


0.933 


0.843 


0.2 


0.932 


0.735 
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Table 1A : 8 Exp in Series, RS = 0.931 (Hi) (Cont...) 

min \ = 0.0002 f/hr, max \ = 0.0016 f/hr, UT = 10 hrs 



S/N 


Test 


Deg of 




Measures of Accuracy 


Plan 


Freedom 


a 


RSLOW 


LEVEL 


4 


Test 15 until 


2*NFC 


0.1 


0.924 


0.970 




7 failed. 


(112) 


0.2 


0.923 


0.931 




NFC =56 


2*(NFC + 
NCOMP) 
(128) 


0.1 


0.915 


0.994 






0.2 


0.913 


0.994 






2*NFC- 

NCOMP 

(104) 


0.1 


0.929 


0.919 






0.2 


0.928 


0.753 






2*(NFC- 

NCOMP) 

(96) 


0.1 


0.934 


0.835 






0.2 


0.933 


0.720 


5 


Test 15 until 


2*NFC 


0.1 


0.915 


0.97S 




3 failed. 


(48) 


0.2 


0.912 


0.975 




NFC =24 


2*(NFC+ 

NCOMP) 

(64) 


0.1 


0.891 


1.000 






0.2 


0.888 


1.000 






2*NFC- 

NCOMP 

(40) 


0.1 


0.927 


0.944 






0.2 


0.926 


0.860 






2*(NFC- 

NCOMP) 

(32) 


0.1 


0.939 


0.753 






0.2 


0.939 


0.634 



simulation are run, the 2 measures of accuracy RSLOW and LEVEL are computed. 
The lower confidence limit procedures are exact if RSLOW = RS in which case 
LEVEL = l-o. 

Table 1A displays the accuracy results for 5 different sampling plans which are 
described in S/N: 1, 2, 3, 4 and 5. 

A comparison of the four values of RSLOW for each of these five sampling 
plans reveal that the lower confidence limit procedure with degrees of freedom equal 
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to 2*NFC-NCOMP is the most accurate lower confidence limit procedure. In S/N 
1, for example, the RSLOW value of 0.927 is the largest such value below the RS 
value of 0.931. Values of RSLOW above RS are optimistic and not as desirable as 
values of RSLOW which are equi-distant below RS. 

The values of RSLOW and LEVEL are based on 1000 replications and their 
accuracy merit should roughly be measured to the nearest one hundredth. That is 
we should round 0.927 to 0.93 and compare it with RS = 0.93. It is evident that the 
lower confidence limit procedure with degrees of freedom equal to 2*NFC-NCOMP 
is very accurate for all cases simulated. 

Table 2 A displays the accuracy results of Case 2 for Type A systems. In this 
case, the 8 components connected in series have the shape parameter P = 2 and the 
scale parameters \ varying between 0.001 and 0.008 failures/hr. Mission time is 10 
hours and each component has this same mission or utilization time (UT). 
Inspection of Table 2A reveals the following: 

(1) More than 5 items of each component should be tested until failure for 
any of these procedures to be reasonably accurate. 

(2) If 15 items of each component are tested until all fail, then these 
procedures will be reasonably accurate when the degrees of freedom is either 
2*NFC-NCOMP or 2*(NFC-NCOMP). 

(3) The procedures are reasonably accurate for 80% confidence level when 
the truncation is not below 7 out of 15 items. 

(4) The procedures are slightly conservative at the 90% confidence level when 
the truncation is 7 out of 15 items or 11 out of 15 items. 

There are numerous ways to modify these lower confidence limit procedures 
to effect improvements in their accuracy. One avenue is to modify the estimate for 
the shape parameter, 0. Some very recent work in the literature provides a method 
for estimating p that differs greatly from the maximum likelihood estimator (MLE) 
and does not require computer iteration. 
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Table 2A : 



8 Wei in Series, RS = 0.980 (Hi) 

min X = 0.001 f/hr, max X = 0.008 f/hr, UT = 10 hrs 



S/N 


Test 

Plan 


Deg of 
Freedom 


a 


Measures of Accuracy 


RSLOW 


LEVEL 


1 


Test 5 until 
5 failed. 

NFC = 40 


2*NFC 

(80) 


0.1 


0.947 


0.992 


0.2 


0.930 


0.989 


2*(NFC+ 

NCOMP) 

(96) 


0.1 


0.937 


0.994 


0.2 


0.918 


0.993 


2*NFC- 

NCOMP 

(72) 


0.1 


0.951 


0.989 


0.2 


0.937 


0.986 


2*(NFC- 

NCOMP) 

(64) 


0.1 


0.956 


0.985 


0.2 


0.983 


0.981 


2 


Test 15 until 
15 failed. 

NFC =120 


2*NFC 

(240) 


0.1 


0.978 


0.918 


0.2 


0.974 


0.913 


2*(NFC+ 

NCOMP) 

(256) 


0.1 


0.977 


0.931 


0.2 


0.978 


0.924 


2*NFC- 

NCOMP 

(232) 


0.1 


0.979 


0.914 


0.2 


0.975 


0.901 


2*(NFC- 

NCOMP) 

(224) 


0.1 


0.980 


0.904 


0.2 


0.975 


0.889 


3 


Test 15 until 
11 failed. 

NFC = 88 


2*NFC 

(176) 


0.1 


0.982 


0.876 


0.2 


0.977 


0.860 


2*(NFC+ 

NCOMP) 

(192) 


0.1 


0.980 


0.981 


0.2 


0.979 


0.882 


2*NFC- 

NCOMP 

(168) 


0.1 


0.983 


0.861 


0.2 


0.978 


0.989 


2*(NFC- 

NCOMP) 

(160) 


0.1 


0.983 


0.840 


0.2 


0.979 


0.819 
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Table 2A : 



8 Wei in Series, RS = 0.980 (Hi) (Cont...) 

min X = 0.001 f/hr, max X = 0.008 f/hr, UT = 10 hrs 



S/N 


Test 

Plan 


Deg of 
Freedom 


a 


Measures of Accuracy 


RSLOW 


LEVEL 




Test 15 until 
7 failed. 

NFC =56 


2*NFC 

(112) 


0.1 


0.987 


0.800 


0.2 


0.981 


0.779 


2*(NFC+ 

NCOMP) 

(128) 


0.1 


0.985 


0.839 


0.2 


0.978 


0.824 


2*NFC- 

NCOMP 

(104) 


0.1 


0.988 


0.776 


0.2 


0.982 


0.753 


2*(NFC- 

NCOMP) 

(96) 


0.1 


0.989 


0.746 


0.2 


0.983 


0.732 


5 


Test 15 until 
3 failed. 

NFC =24 


2*NFC 

(48) 


0.1 


0.994 


0.621 


0.2 


0.991 


0.584 


2*(NFC+ 

NCOMP) 

(64) 


0.1 


0.993 


0.705 


0.2 


0.988 


0.685 


2*NFC- 

NCOMP 

(40) 


0.1 


0.995 


0.548 


0.2 


0.992 


0.514 


2* (NFC- 
NCOMP) 
(32) 


0.1 


0.996 


0.468 


0.2 


0.993 


0.417 



Table 3A displays the results of Case 3 for Type A systems. In this case 8 
components are connected in series. Four of them have failure times with 
exponential distributions and the remaining four have failure times with Weibull 
distributions each with shape parameter, p, equal to 2. 
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Table 3A : 4 Exp and 4 Wei (Mixed) in Series, RS = 0.980 (Hi) 

min A = 0.002 f/hr, max A = 0.008 f/hr, UT = 10 hrs 



S/N 


Test 

Plan 


Deg of 
Freedom 


a 


Measures of Accuracy 


RSLOW 


LEVEL 


1 


Test 5 until 
5 failed. 

NFC =40 


2*NFC 

(80) 


0.1 


0.979 


0.942 


0.2 


0.978 


0.905 


2*(NFC+ 

NCOMP) 

(96) 


0.1 


0.975 


0.987 


0.2 


0.974 


0.976 


2*NFC- 

NCOMP 

(72) 


0.1 


0.981 


0.881 


0.2 


0.980 


0.805 


2*(NFC- 

NCOMP) 

(64) 


0.1 


0.983 


0.771 


0.2 


0.982 


0.684 


2 


Test 15 until 
15 failed. 

NFC =120 


2*NFC 

(240) 


0.1 


0.981 


0.863 


0.2 


0.980 


0.800 


2*(NFC+ 

NCOMP) 

(256) 


0.2 


0.979 


0.941 


0.2 


0.979 


0.898 


2*NFC- 

NCOMP 

(232) 


0.1 


0.981 


0.881 


0.2 


0.980 


0.89* 


2* (NFC- 
NCOMP) 
(224) 


0.1 


0.98* 


0.725 


0.2 


0.981 


0.631 


3 


Test 15 until 
11 failed. 

NFC=88 


2*NFC 

(176) 


0.1 


0.981 


0.864 


0.2 


0.980 


0.801 


2*(NFC + 
NCOMP) 
(192) 


0.1 


0.974 


0.951 


0.2 


0.978 


0.907 


2*NFC- 

NCOMP 

(168) 


0.1 


0.982 


0.802 


0.2 


0.981 


0.898 


2*(NFC- 

NCOMP) 

(160) 


0.1 


0.982 


0.702 


0.2 


0.982 


0.591 
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Table 3 A : 



4 Exp and 4 Wei (Mixed) in Series, RS = 0.980 (Hi) (Cont...) 
min X - 0.002 f/hr, max X = 0.008 f/hr, UT - 10 hrs 



S/N 


Test 

Plan 


Deg of 
Freedom 


a 


Measures of Accuracy 


RSLOW 


LEVEL 


4 


Test 15 until 
7 failed. 

NFC=56 


2*NFC 

(112) 


0.1 


0.981 


0.865 


0.2 


0.980 


0.787 


2*(NFC+ 

NCOMP) 

(128) 


0.1 


0.9*2 


0.952 


0.2 


0.978 


920 


2*NFC- 

NCOMP 

(104) 


0.1 


0.982 


0.769 


0.2 


0.982 


0.676 


2*(NFC- 

NCOMP) 

(96) 


0.1 


0.983 


0.644 


0.2 


0.983 


0.523 


5 


Test 15 until 
3 failed. 

NFC =24 

1 


2*NFC 

(48) 


0.1 


0.982 


0.843 


0.2 


0.981 


0.762 


2*(NFC+ 

NCOMP) 

(64) 


0.1 


0.976 


0.970 


0.2 


0.975 


0.941 


2*NFC- 

NCOMP 

(40) 


0.1 


0.984 


0.644 


0.2 


0.984 




0.580 


2*(NFC- 

NCOMP) 

(32) 


0.1 


0.987 


0.459 


0.2 


0.987 


0.356 



Inspection of Table 3A reveals the following: 

(1) The two procedures corresponding to degrees of freedom equal to 2*NFC 
and 2*NFC-NCOMP are reasonably accurate for all 5 simulation cases. 

(2) The procedures appears to be nearly equally accurate for both 80% and 
90% confidence levels. 
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B. Test Plan 2 : Testing for a Specified Planned Test Time (RETP2) 

In the simulations for test plan 2, components were tested until failure or until 
some planned test time scenario. Failed items were not replaced. Components 
whose failure times had exponential distributions were tested until a pre-determined 
total test time was accumulated for their type of component. Components whose 
failure times had Weibull distributions were tested until failure or a pre-determined 
planned test time for that test item. The latter truncation plan is needed for Weibull- 
type items in order to use the maximum likelihood estimates [as in equation (2.17)] 
to solve for $ . 

Inspection of Table 4A reveals that the lower confidence limit procedure for 
degrees of freedom equal to 2*(1 + NFC) is quite accurate when enough testing is 
done to make the expected number of failures, E[NFC], greater than or equal to 4.8. 
This testing constraint is well within the domain of constraints set on testing in 
development programs for major systems within the Department of Defense. 

Examination of Table 5A reveals that the lower confidence limit procedure for 
degrees of freedom equal to 2*(1 + NFC) is moderately accurate when enough testing 
is done to make E[NFC] greater than or equal to 9. The accuracy diminishes slightly 
as E[NFC] increases. This could be corrected by decreasing the degrees of freedom 
slightly to make RSLOW slightly larger. 

The results displayed in Table 6A show that the lower confidence limit 
procedure for degrees of freedom equal to 2*(1 + NFC) is quite accurate when 
enough testing is done to make E[NFC] greater than or equal to 9.6. 
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Table 4A : 



8 Exp in Series, RS = 0.961 (Hi) 
A = 0.001 f/hr, UT = 5 hrs 



S/N 


Degrees 

of 

Freedom 


K / E[NFC] 
(TT) 


a 


Measures of Accuracy 


RSLOW 


LEVEL 


1 


2*(1 + NFC) 

L 


0.25 / 1.2 
(225) 


0.1 


0.95? 


0.851 


0.1 


0.935 


0.851 


0.5 / 2.4 
(450) 


0.1 


0.957 


0.857 


0.2 


• 

0.954 


0.857 


1.0 / 4.8 
(900) 


0.1 


0.957 


0.914 


0.2 


0.957 


0.850 


2.0 / 9.6 
(1800) 


0.1 


0.958 


0.916 


0.2 


0.960 


0.850 


3.0 / 14.4 
(2700) 


0.1 


0.959 


0.916 


0.2 


0.959 


0.809 


4.0 / 19.2 
(3600) 


0.1 


0.959 


0.937 


0.2 


0.960 


0.843 


5.0 / 24 
(4500) 


0.1 


0.960 


0.926 


0.2 


0.960 


0.814 


10.0 / 48 
(9000) 


0.1 


0.960 


0.924 


0.2 


0.960 


0.809 


20.0 / 96 
(18000) 


0.1 


0.960 


0.914 


0.2 


0.961 


0.820 


30.0 / 144 
(27000) 


0.1 


0.961 


0.906 


0.2 


0.961 


0.804 
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Table 5A : 



8 Wei in Series, RS = 0.956 (Hi) (*) 
\ = 0.005 f/hr, UT = 15 hrs 



S/N 


Degrees 

of 

Freedom 


K / E[NFC] 
(TT) 


a 


Measures of Accuracy 


RSLOW 


LEVEL 


1 


2*(1 + NFC) 


0.25 / 1.2 
(45) 


0.1 


1.000 


0.186 


0.2 


1.000 


0.158 


0.5 / 2.4 
(90) 


0.1 


0.9«« 


0.501 


0.2 


0.979 


0.458 


1.0 / 4.8 
(180) 


0.1 


0.967 


0.767 ' 


0.2 


0.960 


0.732 


2.0 / 9.6 
(360) 


0.1 


0.957 


0.879 


0.2 


0.952 


0.854 


3.0 / 14.4 
(540) 


0.1 


0.952 


0.934 


0.2 


0.946 


0.922 


4.0 / 19.2 
(720) 


0.1 


0.952 


0.940 


0.2 


0.946 


0.928 


5.0 / 24 
(900) 


0.1 


0.952 


0.940 


0.2 


0.946 


0.928 


10.0 / 48 
(1800) 


0.1 


0.952 


0.940 


0.2 


0.946 


0.928 


20.0 / 96 
(3600) 


0.1 


0.952 


0.940 


0.2 


0.946 


0.928 


30.0 / 144 
(5400) 


0.1 


0.952 


0.940 


0.2 


0.946 


0.928 



(*) 20 test items for each Weibull component. 
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Table 6 A : 4 Exp and 4 Wei (Mixed) in Series, RS = 0.958 (Hi) (*) 

A(exp) = 0.001 f/hr, UT(exp) = 5 hrs 
A(wei) = 0.005 f/hr, UT(wei) = 15 hrs 



S/N 


Degrees 

of 

Freedom 


K / E[NFC] 
TT(exp) 
TT(wei) 


a 


Measures of Accuracy 


RSLOW 


LEVEL 


1 


2*(1 + NFC) 


0.25 / 1.2 
(225) 
(45) 


0.1 


1.000 


0.6*4 


0.2 


0.995 


0.451 


0.5 / 2.4 
(450) 
(90) 


0.1 


0.957 


0.8*3 


0.2 


0.975 


0.174 


1.0 / 4.8 
(900) 
(180) 


0.1 


0.971 


0.736 


0.2 


0.965 


0.684 


2.0 / 9.6 
(1800) 
(360) 


0.1 


0.964 


0.803 


0.2 


0.960 


0.765 


3.0 / 14.4 
(2700) 
(540) 


0.1 


0.964 


0.874 


0.2 


0.956 


0.684 


4.0 / 19.2 
(3600) 
(720) 


0.1 


0.960 


0.873 


0.2 


0.957 


0.839 


5.0 / 24 
(4500) 
(900) 


0.1 


0.95# 


0.887 


0.2 


0.956 


0.861 


10.0 / 48 
(9000) 
(1800) 


0.1 


0.95# 


0.891 


0.2 


0.956 


0.862 


20.0 / 96 
(18000) 
(3600) 


0.1 


0.956 


0.892 


0.2 


0.956 


0.867 


30.0 / 144 
(27000) 
(5400) 


0.1 


0.960 


0.877 


0.2 


0.956 


0.862 1 
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The accuracy results of simulations performed in this thesis cannot be extended 
to systems that differ significantly from those simulated here. However, it can be 
said that the procedures which are accurate for series systems are also usually 
accurate for 1-out-of-k parallel systems because system reliability, R s , can be written 
in terms of component reliabilities, R h as 



«, - 1 - n (i - «,> 

i-l 

k 

Thus, upper confidence limits on (1 - R.) will yield a lower confidence 



1-1 



limit on R s . The accuracy of the upper confidence interval procedures for 
* 

Q ( 1 - Rj) , a series-type problem, should be nearly the same as those obtained 

i=i 



in this thesis because equations like (2.9) would be replaced with equations for the 
lower confidence limit on A m and would look like 



m,L(a) 



l-o, 2F 



JEW 



... (4.1) 



1-1 



If the degrees of freedom parameter F is large, as it is in the cases simulated in this 
thesis, the associated Chi-Square distribution with 2 F degrees of freedom is nearly 
symmetric about its mean so the lower tail has nearly the same shape as the upper 
tail. This characteristic should yield very similar accuracies for as was 

obtained for u( ^ a) in this thesis. These accuracy comparisons translate directly 

to similar comparisons about the accuracies of the associated lower confidence limit 
procedures for the system reliability of series and parallel systems. 
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Table 7 : Summary of Procedure Accuracy by Simulation Category (RETP1) 



S/N 


Key Parameters and 
their Levels 


Observation and Discussion 
based on RSLOW and LEVEL 


1 


Run Cases 

■ All Exponential (1) 

■ All Weibull (2) 

■ Mixed (3) 


■ For all exponential systems, accurate 
procedures were developed for component 
sample sizes > 5. 

■ For all Weibull systems, component 
sample sizes should be > 15 with 






truncation at r > 7 failures. 


2 


System Reliabilities 
. Hi (> 0.9) (Type A) 
■ Lo ( > 0.8) (Type B) 


Accurate procedures were developed for 
both cases if sample sizes are adequate. 


3 


Levels of Significance 

■ a = 0.1 

■ a = 0.2 


Accuracy varied slightly depending on 
system type and test plan, but accurate 
procedures exist for both levels. 


4 


Degrees of Freedom 

■ DF = 2*NFC 

. DF = 2*(NFC+NCOMP) 

■ DF = 2*(NFC-NCOMP) 

■ DF = 2*NFC-NCOMP 


Greatest accuracy as follows: 

■ All exponential : 2*NFC-NCOMP 

■ All Weibull : 2*NFC-NCOMP 

or 2*(NFC-NCOMP) 


5 


Test Plan 

■ Test 5 until 5 failures 

■ Test 15 until 15 failures 

■ Test 15 until 11 failures 

■ Test 15 until 7 failures 

■ Test 15 until 3 failures 


Accurate procedures existed for all run 
cases for all system types except for the all 
Weibull system where number of test 
items n > 15 and number of failures r 
should be > 7. 
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Table 8 : Summary of Procedure Accuracy by Simulation Category (RETP2) 



S/N 


Key Parameters and 
their Levels 


Observation and Discussion 
based on RSLOW and LEVEL 


1 


Run Cases 

■ All Exponential (4) 

■ All Weibull (5) 

■ Mixed (6) 


Procedures were accurate for DF = 
2*(1 + NFC) when enough testing was 
done to make the expected number of 
failed components E[NFC] > 9. 


2 


System Reliabilities 

■ Hi (> 0.9) (Type A) 

■ Lo (> 0.8) (Type B) 


Same as S/N(l). 


3 


Levels of Significance 

■ a = 0.1 

■ a = 0.2 


Same as S/N(l). 


4 


Degrees of Freedom 
. DF = 2*(1 + NFC) 

■ DF = 1.3*2*(1 + NFC) 


DF = 2*(1 + NFC) was the most accurate 
procedure. 


5 


Test Plan 
■ K factors of 
0.25, 0.5, 1, 2, 3, 
4, 5, 10, 20, 30 


K should be chosen so that E[NFC] > 9. 



Tables 7 and 8 provide cursory summaries of some constraints needed to assure 
the existence of one or more accurate lower confidence limit procedures among the 
procedures that were evaluated. The simulation scenarios are divided into five 
categories for this summarization. 
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V. APPLICATION EXAMPLES 



Based on the procedures evaluated by the RETP1 and RETP2 runs, four 
different test plans and failure time data were constructed to illustrate the use of the 
procedures in providing a lower lOO(l-a) % confidence limit for the system reliability 
of a series system with different types of components. 



CASE 1 : 8 Exponential Components in Series 

TEST PLAN 1 - Test 15 until 7 fails for each component 

I . Raw Data 



Comp Ordered Failure Times (h) 



i 


T( 1) 


T(2) 


T(3 ) 


T(4 ) 


T( 5) 


T(6) 


T(7) 


1 


300.0 


400.0 


500.0 


600.0 


700.0 


800.0 


900.0 


2 


350.0 


450.0 


550.0 


650.0 


750.0 


850.0 


950.0 


3 


400.0 


500.0 


600.0 


700.0 


800.0 


900.0 


1000.0 


4 


450.0 


550.0 


650.0 


750.0 


850.0 


950.0 


1050.0 


5 


500.0 


600.0 


700.0 


800.0 


900.0 


1000.0 


1100.0 


6 


550.0 


650.0 


750.0 


850.0 


950.0 


1050.0 


1150.0 


7 


600.0 


700.0 


800.0 


900.0 


1000.0 


1100.0 


1200.0 


8 


650.0 


750.0 


850.0 


950.0 


1050.0 


1150.0 


1250.0 



II. Data Summary 
Comp 

i UT(i) 


NC(i) 


NF(i) 


TT(i ) 


ELM(i) 


ER(i ) 


ER(i)* 
TT(i ) 


1 


5.0 


15 


7 


11400.0 


0.00053 


1.00000 


11400.0 


2 


5.0 


15 


7 


12150.0 


0.00049 


0.93827 


11400.0 


3 


5.0 


15 


7 


12900.0 


0.00047 


0.88372 


11400.0 


4 


5.0 


15 


7 


13650.0 


0.00044 


0.83516 


11400.0 


5 


5.0 


15 


7 


14400.0 


0.00042 


0.79167 


11400.0 


6 


5.0 


15 


7 


15150.0 


0.00040 


0.75248 


11400.0 


7 


5.0 


15 


7 


15900.0 


0.00038 


0.71698 


11400.0 


8 


5.0 


15 


7 


16650.0 


0.00036 


0.68468 


11400.0 



III. Estimation Procedure for RSLOW 



Parameter df Value 



ALPHA 0 . 1 

NFC 56 

CHISQD 112 131.56 (df = 2 * NFC , CHISQD from tables) 

LMU 0.00072 



RSLOW 0 . 97647 
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CASE 2 : 8 Weibull Components in Series 

TEST PLAN 1 - Test 15 until 7 fails for each component 



I . Raw Data 



Comp 

i 


T ( 1) 


T(2) 


Ordered 

TC3) 


Failure 
T(4 ) 


Times (h) 
T(5) 


T(6) 


T(7 ) 


1 


10.0 


20.0 


30.0 


40.0 


50.0 


60.0 


70.0 


2 


20.0 


30.0 


40.0 


50.0 


60.0 


70.0 


80.0 


3 


30.0 


40.0 


50.0 


60.0 


70.0 


80.0 


90.0 


4 


40.0 


50.0 


60.0 


70.0 


80.0 


90.0 


100.0 


5 


50.0 


60.0 


70.0 


80.0 


90.0 


100.0 


110.0 


6 


60.0 


70.0 


80.0 


90.0 


100.0 


110.0 


120.0 


7 


70.0 


80.0 


90.0 


100.0 


110.0 


120.0 


130.0 


8 


80.0 


90.0 


100.0 


110.0 


120.0 


130.0 


140.0 



II. Data Summary 



Comp ERCi)* 



i 


UT(i) 


NC( i ) 


NF (i ) 


TTCi) 


ELM(i) 


ER(i) 


TT(i) 


1 


5.0 


15 


7 


5.3E+04 


1.32E-04 


1.00000 


5.3E+04 


2 


5.0 


15 


7 


7.2E+04 


9.79E-05 


0.74406 


5.3E+04 


3 


5.0 


15 


7 


9.3E+04 


7.54E-05 


0.57328 


5 . 3E+04 


4 


5.0 


15 


7 


1.2E+05 


5.98E-05 


0.45431 


5.3E+04 


5 


5.0 


15 


7 


1.4E+05 


4.85E-05 


0.36842 


5.3E+04 


6 


5.0 


15 


7 


1.7E+05 


4.01E-05 


0.30452 


5.3E+04 


7 


5.0 


15 


7 


2.1E+05 


3.37E-05 


0.25577 


5 . 3E+04 


8 


5.0 


15 


7 


2.4E+05 


2.87E-05 


0.21777 


5.3E+04 



III. Estimation Procedure for RSLOW 



Parameter 


df 


Value 


BETA 




2.0 


ALPHA 




0.1 


NFC 




56 



CHISQD 

LMU 


112 131.56 (df « 2 * NFC , CHISQD from tables) 

1.55E-04 


RSLOW 


0.98497 



IV. Workarea 



Comp Ordered Failure Times (h) raised to the power of BETA 



i 


T ' ( 1 ) 


T'(2) 


T’(3) 


T ' (4 ) 


T'(5) 


T* (6) 


T* (7) 


1 


1.0E+02 


4.0E+02 


9.0E+02 


1.6E+03 


2.5E+03 


3.6E+03 


4.9E+03 


2 


4.0E+02 


9.0E+02 


1.6E+03 


2.5E+03 


3.6E+03 


4.9E+03 


6.4E+03 


3 


9 . 0E+02 


1.6E+03 


2.5E+03 


3.6E+03 


4 . 9E+03 


6.4E+03 


8.1E+03 


4 


1.6E+03 


2.5E+03 


3.6E+03 


4.9E+03 


6.4E+03 


8.1E+03 


1.0E+04 


5 


2.5E+03 


3.6E+03 


4.9E+03 


6.4E+03 


8.1E+03 


1.0E+04 


1.2E+04 


6 


3 . 6E+03 


4 . 9E+03 


6.4E+03 


8.1E+03 


1.0E+04 


1.2E+04 


1.4E+04 


7 


4.9E+03 


6.4E+03 


8.1E+03 


1.0E+04 


1.2E+04 


1.4E+04 


1.7E+04 


8 


6.4E+03 


8.1E+03 


1.0E+04 


1.2E+04 


1.4E+04 


1.7E+04 


2.0E+04 
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CASE 3 : 8 Exponential Components in Series 
TEST PLAN 2 ~ Test until TT(i) for each component 

I . Raw Data 



No data except noting the number of failures for each component (NF(i)) 

II. Data Summary 



Comp 

i 


UT(i) 


NF( i ) 


TT(i ) 


ELM(i) 


ER(i) 


ER(i)* 

TT(i) 


ER(i)* 

UT(i) 


1 


5.0 


6 


2000.0 


0.00300 


0.85714 


1714.3 


4.3 


2 


5.0 


6 


2000.0 


0.00300 


0.85714 


1714.3 


4.3 


3 


5.0 


5 


2000.0 


0.00250 


0.71429 


1428.6 


3.6 


4 


5.0 


7 


2000.0 


0.00350 


1.00000 


2000.0 


5.0 


5 


5.0 


5 


2000.0 


0.00250 


0.71429 


1428.6 


3.6 


6 


5.0 


5 


2000.0 


0.00250 


0.71429 


1428.6 


3.6 


7 


5.0 


4 


2000.0 


0.00200 


0.57143 


1142.9 


2.9 


8 


5.0 


7 


2000.0 


0.00350 


1.00000 


2000.0 


5.0 



III. Estimation Procedure for RSLOW 



Parameter df Value 



ALPHA 0 . 1 

NFC 45 

CHISQD 92 109.76 (df *= 2 * ( 1 + NFC ) , CHISQD from tables) 

LMU 0.00427 



RSLOW 0.87180 
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CASE A ; 8 Weibull Components in Series 

TEST PLAN 2 - Test 20 until TT(i) for each component 



I . Raw Data 



Comp 

i 


TT(i ) 


T(l) 


Ordered 

T(2) 


Failure 
T(3 ) 


Times (h) 
T( A ) 


T(5) 


T(6 ) 


1 


60.0 


10.0 


20.0 


30.0 


A0 . 0 


50.0 


60.0 


2 


60.0 


15.0 


25.0 


35.0 


A5.0 


55.0 




3 


60.0 


20.0 


30.0 


A0 . 0 


50.0 


60.0 




A 


60.0 


25.0 


35.0 


A5.0 


55.0 






5 


60.0 


30.0 


A0 . 0 


50.0 


60.0 






6 


60.0 


35.0 


A5.0 


55.0 








7 


60.0 


A0 . 0 


50.0 


60.0 








8 


60.0 


A5.0 


55.0 











II. Data Summary 



Comp 

i 


UT(i) 


NC(i ) 


NF(i) 


ET(i ) 


ELM(i) 


ER(i) 


ER(i)* 

ET(i) 


1 


5.0 


20 


6 


6.0E+0A 


1.01E-0A 


1.00000 


6.0E+0A 


2 


5.0 


20 


5 


6.1E+0A 


8.18E-05 


0.81118 


5 . 0E+0A 


3 


5.0 


20 


5 


6.3E+0A 


7.9AE-05 


0.7870A 


5.0E+0A 


A 


5.0 


20 


A 


6.5E+0A 


6.20E-05 


0.61A99 


A . 0E+0A 


5 


5.0 


20 


A 


6.6E+0A 


6.0AE-05 


0.59919 


A . 0E+0A ' 


6 


5.0 


20 


3 


6.7E+0A 


A.A5E-05 


0.AA090 


3.0E+0A 


7 


5.0 


20 


3 


6.9E+0A 


A.35E-05 


0.A3179 


3.0E+0A 


8 


5.0 


20 


2 


7.0E+0A 


2.86E-05 


0.2839A 


2.0E+0A 



III. Estimation Procedure for RSLOW 



irameter 


df 


Value 




BETA 




2.0 




ALPHA 




0.1 




NFC 
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CHISQD 


66 


81.09 (df “ 2 * ( 1 + NFC ) 


, CHISQD from tables) 


LMU 




1.28E-0A 




RSLOW 




0.98A25 





IV. Workarea 



Comp Ordered Failure Times (h) raised to the power of BETA 



i 


TT' (1) 


T ' (1) 


T ' (2) 


T’ (3 ) 


T ' ( A ) 


T'(5) 


T ' (6) 


1 


3.6E+03 


1.0E+02 


A . 0E+02 


9.0E+02 


1.6E+03 


2 . 5E+03 


3.6E+03 


2 


3.6E+03 


2.3E+02 


6.3E+02 


1.2E+03 


2.0E+03 


3.0E+03 




3 


3.6E+03 


A.0E+02 


9.0E+02 


1.6E+03 


2.5E+03 


3.6E+03 




A 


3.6E+03 


6.3E+02 


1.2E+03 


2.0E+03 


3.0E+03 






5 


3.6E+03 


9.0E+02 


1.6E+03 


2.5E+03 


3.6E+03 






6 


3.6E+03 


1.2E+03 


2 . 0E+03 


3.0E+03 








7 


3.6E+03 


1.6E+03 


2.5E+03 


3.6E+03 








8 


3.6E+03 


2.0E+03 


3.0E+03 
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VI. CONCLUSION 



Some of the lower confidence limit procedures developed and evaluated in this 
thesis are reasonably accurate for the series systems simulated for test plans with 
sample sizes and truncation scenarios that are usually experienced in DoD aquisition 
programs. The accuracy of these methods can be varied by modifying the degrees 
of freedom parameter, F, in the x 2 b ,f term in equation (2.9), namely: 

c _ X 2 0j2 f 

*nt,U( a) “ 7 

2EW 

1=1 

The computer program can be modified with modest effort to accomodate specific 
complex coherent systems so long as the components have failure time distributions 
that are exponential or Weibull. This means that the computer program provided in 
this thesis can be used to develop a reasonably accurate lower confidence limit for 
the system reliability of a specific complex quasi-coherent system with independent 
components. This can be done by choosing the failure distribution and associated 
parameters of the components, the corresponding test plan parameters and the 
desired level of confidence. The simulation can then be run for this set of 
parameters for various equations for the degrees of freedom parameter, F, to 
determine an equation for F that yields a lower confidence limit with a satisfactory 
degree of accuracy. 

When testing is truncated on the number of failures, a reasonably accurate 
procedure existed for all 3 cases of systems (all exponential, all Weibull and mixed ) 
that were simulated when (a) the sample size of the components was 10 or larger, 
and, (b) the ratio of the number of failures to the sample size was at least 0.5. When 
testing was truncated by planned test time, reasonably accurate procedures were 
found for cases where the expected number of failures was at least 7. 
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VII. RECOMMENDATIONS 

The computer program developed in this thesis facilitates the development of 
lower confidence limit procedures for explicit quasi-coherent systems. Systems other 
than series systems with large numbers of components (eg. 30) should be simulated 
to test the versatility of the general lower confidence limit methods used here. 

Modified estimates for the /J (shape) parameter in the Weibull failure time 
distribution and the parameter r = \/\ m should be explored in an attempt to find 
more accurate procedures. 
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APPENDIX A : Derivation of Formula Used 



Suppose T has an exponential distribution with failure rate X. Suppose T (» 
T (3) ,... , T (r) are the first r ordered statistics in a random sample of size n from this 
exponential distribution. Let 5 be defined by 

1-1 

It is well known that 2X5 has a Chi-Square distribution with 2 r degrees of freedom 
(See Ref.3 p 488). The maximum likelihood estimator for X is given by 

X = L 

s 

We seek an unbiased estimator for i . Suppose X has a Chi-Square distribution 
with 2 r degrees of freedom, then 

Ux ’ r)s Tl 

and the integral of this function from zero to infinity equals 1 (See Ref.3). 



oo 

E{— ) = f — i — x (r ‘b~ 1 exp(-— )dx 
* J 0 2T(r) 2 

= rCr-!)? 1 (r- i)-i flvp( .^y y 

1 2 r_1 r(r-l) ^ 2 } 



2W J 
1 

2(r-l) 
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Then 



£ <*> ■ - r2X£( 2^5 ) ' 



r\ 

Fi) 



Therefore 



£(— *) - E(Q - * 

r S 
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APPENDIX B : Users’ Guide for RETP1 



Reliability Estimation Test Plan 1 (RETP1). 
by YEE, Kah-Chee 
July 91 



1. Brief Description. 

RETP1 is a computer program written in FORTRAN that runs on the Amdahl 
mainframe at NPGS. It allows the user to simulate exponential and Weibull failure 
times of component items being tested to evaluate the accuracy of a confidence limit 
estimation procedure based on Type II data censoring (that is, testing n^ items of 
component i until fj of them failed). 

2. Program Input. (IN1.DATA) 

The input of the program are specified to the program via an input file called 
IN1.DAT. A sample input file is shown below. 



This file contains the inputs required by the RETP1 model. 

Update only the numerical values between dotted lines as appropriate. 
Do not delete any of the comment lines. (INI. DAT) 14 May 91 



Value 


Type 


Units 


Description 


Variable 


16807.0 


Real 


- 


initial random seed 


ISEED 


8 


Int 


- 


total # of components in system 


NCOMP 


4 


Int 


- 


# of Exponential components 


NEXP 


4 


Int 


- 


# of WEIbull components 


NWEI 


0 


Int 


- 


# of GEOmetric components 


NGEO 


0.01 


Real 


- 


tolerance for MLE 


TOL 


0.20 


REAL 


- 


DESIRED SIGNIFICANCE LEVEL 


ALPHA 


1000 


Int 


- 


# of replications desired 


NREP 


3 


Int 




test case number 

1 - all EXP 

2 * all WEI 

3 * EXP + WEI 

4 = EXP + WEI + GEO 


TCN 


8 


Int 


- 


number of cut sets 


NCS 
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C TEST PLAN : Testing NC(I) items of component i 

C until NF(I) of them fails. (Use REAL numbers ONLY!!!) 



Comp 


Comp 


Comp Parameters 


Util 


Test Plan Inputs 


Number 


Type 


Scale 


Shape 


Time/Cycle 


# Comp 


# Failed 


I 


TY(I) 


PARMl(I) 


PARM2 ( I ) 


UT(I) 


UC(I) 


NC(I) 


NF(I) 


Int 


Int 


Real 


Real 


(hrs) 


Int 


Int 


Int 


1.0 


1.0 


0.0020 


1.0 


10.0 




15.0 


3.0 


2.0 


1.0 


0.0040 


1.0 


10.0 




15.0 


3.0 


3.0 


1.0 


0.0060 


1.0 


10.0 




15.0 


3.0 


4.0 


1.0 


0.0080 


1.0 


10.0 




15.0 


3.0 


5.0 


2.0 


0.0020 


2.0 


10.0 




15.0 


3.0 


6.0 


2.0 


0.0040 


2.0 


10.0 




15.0 


3.0 


7.0 


2.0 


0.0060 


2.0 


10.0 




15.0 


3.0 


8.0 


2.0 


0.0080 


2.0 


10.0 




15.0 


3.0 


Note : 


TY(I ) 


=1 EXPONENTIAL P(surv) = 


exp(-PARMl)*T) 






TY(I) 


=2 WEI BULL P(surv) = 


exp ( - ( PARM1*T ) **PARM2 ) 




TY(I) 


-3 GEOMETRIC P(surv) - 


PARM1**N 






SYSTEM 


CONFIGURATION : 


Identification of 


CUT SETs 












- min groups 


of components 


that have 


to fail 








for the system to fail. 






Cut Set 


# in Set 


List of Components in Cutset 




J 


C0MP(J ,1) 


C0MP(J ,2) . 


. . up 


to C0MP(J,1) components 


i 




1 


10 0 0 


0 0 


0 0 0 






2 




1 


2 0 0 0 


0 0 


0 0 0 






3 




1 


3 0 0 0 


0 0 


0 0 0 






4 




1 


4 0 0 0 


0 0 


0 0 0 






5 




1 


5 0 0 0 


0 0 


0 0 0 






6 




1 


6 0 0 0 


0 0 


0 0 0 






7 




1 


7 0 0 0 


0 0 


0 0 0 






8 




1 


8 0 0 0 


0 0 


0 0 0 







3. Program Flow and Logic. (NAME1.DEF, PARM1.DEF and RETP1.FOR) 

Input parameters are first read in by the program by calling the INPUT 
subroutine. The program then evoke the SIM subroutine which generates the 
random failure times and compute the key statistics required in the procedure. The 
next subroutine EVAL determines the measures of accuracy for run. REPORT is 
the subroutine which generates the output file for the run OUTl.DAT. 

The variables in the program RETP1.FOR are described in the file 
NAME1.DEF as listed below. 



41 



c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 



This file contains the declaration for input and output variables 
used in the the RETP1 model. (NAME1.DEF) 14 May 91 



Input Variables. 



ISEED = initial random seed selected. 

SEED * current random seed. 

RS = true overall series system reliability. 

ALPHA = level of significance desired. 

NREP - number of replications desired for the simulation. 
TPN - test plan number (1). 

TCN — test case number (1, 2, 3 or 4). 

NCOMP - total number of components in the system. 

NEXP - number of components with EXP failure times. 

NWEI = number of components with WEI failure times. 

NGEO — number of components with GEO failure times. 

TOL — desired tolerance for MLE of WEI shape parameter. 



Distribution 
TY(I) - type 
PARM (1,1) 
P ARM (2,1) 
UT(I) 

UC(I) 



Exponential 

1 

Scale (1/hr) 



WEIbull GEOmetric 

2 3 

Scale (1/hr) Prob 

Shape 

utilization time (hrs) for component i (EXP and WEI), 
utilization cycles for component i (GEO only) . 



NC(I) * number of test samples (sample size) for cmponent i. 
NF(I) * desired number of failures in test for component i. 

NCS = number of cut-sets for the system. 

COMP(J,K) * kth parameter of cut-set j (first being the no. of 
components belonging to the cut-set). 



Assumed Variables. 



MAXCOMP -= maximum number of components allowed in the system. 
MAXREP = maximum number of replications permitted. 

MAXCUT = maximum number of cut-sets. 

Program and Output Variables. 



RS = true overall system reliability. 

TT(I) = total accumulated failure time (hr) for component i 
(EXP and WEI only). 

TC(I) = total accumulated cycles to failure (incl. failure cycle) 
for component i (GEO only) . 

EBETA(I) * estimate for shape parameter of component i (if Weibull) . 
REL1(J) -= actual reliability for cut-set j. 

REL2(J) “ computed reliability for cut-set j for current replication. 
ELM(I) = estimated component failure rate (1/hrs) for component i. 
ELMAX(M) = max estimated component failure rate for rep m (1/hrs). 
ER(I) = ratio of estimated failure rate to LMAX. 

NFC(M) = total number of failed test components. 
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C LMU(M) = upper confidence bound for failure rate (1/hrs) , 

C RSL(M) = lower confidence limit estimated for system reliability 
C ... for the mth replication. 

C ORSL(M) = ordered RSL(M) (ascending) . 

C RSLOW = (1 -ALPHA) xlOO percentile of set of RSL(M) . 

C LEVEL «= achieved confidence level, ie. proportion of RSL(M) that 
C are lesser than RS (conservative estimate). 

C 

C-- END OF NAME1.DEF 

C 



Together with the main program in RETP1.FOR are the other subroutines 
needed in the simulation. The declaration of variables is done in the file 
PARM1.DEF. Relevant descriptions are included as comment lines in the source 
code to help explain the program segments. A listing of PARM1.DEF and 
RETP1.FOR is given below. 

c 

C This file contains the declaration for input and output variables 
C used in the the RETP1 model. (PARM1.DEF) 14 May 91 

C 

INTEGER MAXCOMP, MAXREP 

PARAMETER( MAXCOMP - 100 , MAXREP = 1000 , MAXCUT - 20 ) 

REAL*8 ISEED, SEED 

INTEGER NREP , TCN, NCOMP, NEXP, NWEI , NGEO, NCS , 

* NC (MAXCOMP) , NF (MAXCOMP) , TY ( MAXCOMP ) , NFC (MAXREP ) , 

* UC (MAXCOMP) , TC (MAXCOMP) , COMP (MAXCUT, MAXCOMP) 

REAL* 8 RS, ALPHA, UT (MAXCOMP ) , TT (MAXCOMP ) , 

* PARM(2 .MAXCOMP) , ELM ( MAXCOMP ) , ER (MAXCOMP ) , 

* LMU(MAXREP) , RSL(MAXREP) , ORSL(MAXREP) , 

* ELMAX (MAXREP) , RSLOW, LEVEL, T0L, EBETA(MAXCOMP) , 

* REL1 (MAXCUT) , REL2 (MAXCUT) 

C 

C0MM0N/BL0CK1/ISEED , SEED, NREP, TCN, NCOMP, NC, NF, NEXP, NWEI, 

* NGEO, NCS, TY, NFC, UC , TC , COMP 
C0MM0N/BL0CK2/RS , ALPHA, UT, TT, PARM, ELM, ER, LMU , 

* RSL, ORSL, ELMAX, RSLOW, LEVEL, TOL, EBETA, 

* REL1 , REL2 
C 

C-- END OF PARM1.DEF 

C 

C 

C This file contains the main program and the subroutines 
C for the Reliability Estimation Test Plan 1 (RETP1) model. 

C (RETP1.FOR) - runs on a IBM PC Compatible. 

C 

C MAINFRAME VERSION. . . 

C 

C Test Plan 1 : Testing NC(I) items for component i 
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C Until NF(I) of them fails. 

C 

C by Yee Kah-Chee SMC 2802. 

C 14 May 91. 

C 

C 1. Main Program (RETP1) . 

C 

PROGRAM RETP1 

Include the declaration files. 

INCLUDE 'NAME1 DEF' 

INCLUDE ' PARM1 DEF' 

Read in input data. 

CALL INPUT 

Activate simulation. 

CALL SIM 

Process and evaluate output data. 

CALL EVAL 

Generate simulation report. 

CALL REPORT 

STOP 
END 



2. Input Initialisation Subroutine (INPUT). 

SUBROUTINE INPUT 

This subroutine reads in the inputs for the RETP1 model. 

Include the declaration file. 

INCLUDE ' PARM1 DEF' 

INTEGER I, J, K, DUM2(11) 

REAL* 8 DUM1 ( 7 ) 

Read data from 'INI. DAT' designated as logic unit 1. 

OPEN (UNIT=1) 
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c 



c 



c 



c 



c 



READ(l.lO) 

10 FORMAT (IX ,/////) 

READ ( 1 , * ) I SEED 
READ(1 , *) NCOMP 
READ(1 ,*) NEXP 
READ(1 , *) NWEI 
READ(1 , *) NGEO 
READ(1 , *) TOL 
READ(1 ,*) ALPHA 
READ(1 , *) NREP 
READ(1 , *) TCN 
READ(1 , 20) 

20 FORMAT (IX,///) 

READ(1 , *) NCS 
READ(1 , 30) 

30 FORMAT (IX, ////////) 

DO 50 K = 1, NCOMP 
READ(1,*) DUM1 

I = NINT(DUM1(1) ) 

TY(I) = NINT(DUM1(2) ) 

IF (TY(I).EQ.l) THEN 
PARM(1 , I ) - DUM1(3) 
PARM(2 , I ) - DUM1 (4) 
UT(I) = DUM1 ( 5 ) 

NC(I) - NINT(DUM1(6) ) 
NF(I ) = NINT(DUM1(7) ) 
EBETA(I) = 0 

ELSEIF (TY(I) . EQ. 2) THEN 
PARM (1,1) = DUM1 ( 3 ) 
PARM(2 , I) - DUM1 ( 4 ) 
UT(I) - DUM1 ( 5 ) 

NC(I) - NINT(DUM1(6) ) 
NF(I) = NINT(DUM1(7) ) 
ELSEIF ( TY ( I ) . EQ . 3 ) THEN 
PARM (1,1) = DUM1 ( 3 ) 
PARM (2 , I) - DUM1(4) 
UC(I) = NINT(DUM1(5) ) 
UT(I) = DUM1 ( 5 ) 

NC(I) - NINT(DUMl(6) ) 
NF(I) - NINT(DUM1(7) ) 
END IF 

50 CONTINUE 
READ(1 , 60) 

60 FORMAT (IX,///////////) 

DO 80 I - 1, NCS 
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READ(1 , *) DUM2 
J = DUM2(1) 

COMP(J , 1) = DUM2 ( 2 ) 

DO 70 K = 1, C0MP(J , 1) 
C0MP(J ,K+1) = DUM2(K+2) 
70 CONTINUE 
80 CONTINUE 

CLOSE (UNIT-1) 

RETURN 

END 



3. Subroutine (SIM) 



SUBROUTINE SIM 

This subroutine simulates NREP possible outcomes of the test plan 
desired in order to obtain the raw estimates of LMU(M) and RSL(M) 
for each of the replication. 

Include the declaration file 
and declare local variables. 

INCLUDE 'PARM1 DEF' 

INTEGER I, J, K, M, ISUM, KEY 

REAL* 8 UNI 

REAL* 8 SUM, PROD, 

* FT(MAXCOMP) , OFT (MAXCOMP) 

SEED = ISEED 

Compute overall true system reliability RS . 

RS - 1.0 
DO 30 J = 1, NCS 
PROD =1.0 

DO 20 I = 1, C0MP(J , 1) 

K = C0MP(J , 1+1) 

PROD = PROD* ( 1 - SURV (TY (K) , PARM(1 , K) , PARM(2 , K) , UT (K) ) ) 
20 CONTINUE 

REL1(J) - 1.0 - PROD 
RS = RS * REL1(J) 

30 CONTINUE 

Start of Simulation. 

(Initialize replication counter M) . 

M = 1 

DO WHILE (M.LE.NREP) 

PRINT 33, M 



46 



n n n 



35 FORMAT( IX, 'Replication ',14) 



C 

C 

C 

C 

C 

C 

C 

C 

C 



C 



Test Plan : Sample and determine unknown TT(I) 

with known NC(I) until NF(I) fails. 

Generate NC(I) failure times, put them in ascending order 
with the smallest failure time on the top of the list. 

DO 70 I - 1, NCOMP 

DO 40 K = 1, NC(I) 

CALL LRNDPC (SEED, UNI ,1) 

IF(TY(I) .EQ.l) THEN 

FT(K) - -LOG (UNI) /PARM(1 , I) 

ELSEIF(TY(I) . EQ. 2) THEN 

FT(K) - (1 . 0/PARM(l , I) )*( -LOG (UNI) )**(1 . 0/PARM(2 , I) ) 
ELSEIF(TY(I) . EQ. 3) THEN 
FT(K) = 1.0 

DO WHILE (UNI .LT. PARM(1 , I) ) 

FT(K) - FT(K) + 1.0 
CALL LRNDPC (SEED, UNI, 1) 

ENDDO 
END IF 

40 CONTINUE 



Bubble Sort the failure times in ascending order. 



CALL BUBBLE(NC(I) , FT, OFT) 

C 

C Compute the total time accumulated in the test and the estimate 
C for the failure rate of the component as in the procedure. 

C 

IF (TY(I) .NE.2) THEN 
SUM - 0.0 

DO 50 K = 1, NF(I) 

SUM = SUM + OFT(K) 

50 CONTINUE 

TT(I) - FLOAT (NC ( I ) -NF(I) )*OFT (NF(I ) ) + SUM 
IF (TY(I).EQ.l) THEN 

ELM(I) -= FLOAT (NF(I)-l) /TT ( I ) 

ELS El F (TY(I) . EQ . 3) THEN 

ELM(I) - FL0AT(NF(I)-1)/TT(I) 

END IF 

ELS El F (TY(I) . EQ. 2) THEN 
C 

C PRINT 55, M, I 

C 55 FORMAT (IX, 'REP - ',13, ' Comp - ',13,/) 

C 

CALL MLESHAPE(OFT,NC(I) ,NF(I) , TOL , 1 . DO , EBETA( I ) ) 
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c 

c 



c 



60 



C 



70 



EBETA(I) = BN (NC ( I ) ) *EBETA( I ) 

SUM =0.0 

DO 60 K - 1, NF(I) 

SUM = SUM + OFT(K)**EBETA(I) 

CONTINUE 

TT(I) - FLOAT (NC(I) -NF(I) )*0FT(NF(I) )**EBETA(I) + SUM 
ELM(I) = FL0AT(NF(I))/TT(I) 

ENDIF 

CONTINUE 



Determine the total number of failed test items. 



I SUM = 0 

DO 80 I - 1, NC0MP 
ISUM = I SUM + NF(I) 
80 CONTINUE 

NFC(M) = ISUM 



Determine the maximum failure rate estimate 
and identify that component. 



ELMAX(M) = 0.0 
KEY = 0 

DO 90 I = 1, NCOMP 

IF (ELM(I) .GT.ELMAX(M)) THEN 
ELMAX(M) = ELM(I) 

KEY = I 
ENDIF 

90 CONTINUE 

Compute the ratios of the failure rate estimates to their maximum. 

DO 100 I = 1, NCOMP 

ER(I) = ELM(I) / ELMAX(M) 

100 CONTINUE 



Determination of LMU(M) 

SUM = 0.0 

DO 110 I - 1, NCOMP 

SUM = SUM + (ER(I)*TT(I) ) 

110 CONTINUE 

LMU(M) = CHISQD(1 -ALPHA, 2*(NFC(M) -NCOMP) )/(2*SUM) 



Compute estimate of overall reliability RSL(M) for the system. 
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RSL(M) = 1.0 
DO 130 J - 1, NCS 
PROD -1.0 

DO 120 I - 1, C0MP(J , 1) 

K - C0MP( J , 1+1) 

IF (TY(K).EQ.l) THEN 

PROD - PROD* ( 1 - SURV(TY(K) , LMU(M)*ER(K) ,EBETA(K) ,UT(K) ) ) 
ELSEIF (TY(K) . EQ . 2) THEN 

PROD - PROD* ( 1 - SURV(TY(K) , (LMU(M)*ER(K) )**(1 . /EBETA(K) ) , 
* EBETA(K) ,UT(K) ) ) 

ELSEIF (TY(K) . EQ. 3) THEN 

PROD = PROD* ( 1 - SURV(TY(K) , 1 . D0-LMU(M)*ER(K) , 0 . DO ,UT(K) ) ) 
ENDIF 

120 CONTINUE 

REL2(J) = 1.0 - PROD 
RSL(M) = RSL(M) * REL2(J) 

130 CONTINUE 



Increment replication counter. 
M = M + 1 



ENDDO 

END 



4. Subroutine (EVAL). 



SUBROUTINE EVAL 

This subroutine calls BUBBLE to sort the array RSL(NREP) in 
ascending order to get an ordered array ORSL(NREP) . It also 
determine the estimate for RSLOW at the specified significance 
level ALPHA and the value of LEVEL in which ORSL(LEVEL) is closest 
to the true reliability RS . 

Include the declaration files 
and declare the local variables. 

INCLUDE ' PARM1 DEF' 

INTEGER INDEX 
REAL* 8 DIFF 

Order the array RSL(NREP) in ascending order. 

DO 10 M - 1, NREP 
ORSL(M) = RSL(M) 

10 CONTINUE 

Bubble Sort. Sink the larger of the pair. 
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c 

CALL BUBBLE (NREP , RSL, ORSL) 

C 

C Determine the (1-ALPHA)% lower confidence bound for the system 
C reliability. 

C 

RSLOW - 0RSL(NINT(NREP*(1-ALPHA) ) ) 

C 

C Finding the % confidence level for the true reliability RS. 

C (ie. the proportion of RSL(M) lesser than RS) 

C 

DIFF - 1.0 
INDEX = 0 

DO 200 M = 1, NREP 

IF (ABS(ORSL(M)-RS) .LT.DIFF) THEN 
DIFF - ABS(0RSL(M) -RS) 

INDEX = M 
ENDIF 

200 CONTINUE 
C 

LEVEL = FLOAT (INDEX) /NREP 

Record evaluated parameters in RAW1.DAT (unit 2). 

0PEN(UNIT=2) 

WRITE(2 , 300) 

300 FORMAT (IX, ' M LMU(M) ELMAX(M) RSL(M) ' , 

* ' ORSL(M) NFC(M) ' ) 

DO 500 M = 1, NREP 

WRITE(2 ,400) M,LMU(M) , ELMAX(M) ,RSL(M) ,0RSL(M) ,NFC(M) 

400 F0RMAT(1X, 16 , 2F12 . 7 , 2F12 . 7 , 110) 

500 CONTINUE 

CLOSE (UNIT=2) 

RETURN 
END 



5. Report Generation Subroutine (REPORT). 



SUBROUTINE REPORT 
C 

C This subroutine record the simulation results into the 'OUTl.DAT' 
C file as logic unit 3. 

C 

C Include the declaration files 
C and declare local variables. 

C 

INCLUDE ' PARM1 DEF' 

INTEGER I, J, K, DUM(10) 

C 
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Write to output file '0UT1.DAT' designated as logic unit 3. 



OPEN (UNIT— 3) 

C 

WRITE( 3 , 10) 

WRITE(3 , 20) NREP 
WRITE( 3 ,25) 

WRITE(3 , 26) 

WRITE( 3 , 30) 

WRITE(3 ,40) 

WRITE (3,50) I SEED , NCOMP , ALPHA , TOL , NCS , TON 
C 



WRITE( 3 , 60) 

DO 200 I - 1, NCOMP 

WRITE(3 ,70) I , TY(I) , PARM(1 , I) ,PARM(2,I) ,UT(I) ,NC(I) ,NF(I) 

200 CONTINUE 
WRITE( 3 , 80) 

WRITE(3 , 90) 

DO 300 I = 1, NCOMP 

WRITE(3 , 100) I ,NF(I) ,TT(I) ,ELM(I) , ER(I) , EBETA(I) 

300 CONTINUE 

WRITE(3 , 110) 

WRITE(3 , 120) 

DO 500 J - 1, NCS 
DO 400 K - 1, 10 
DUM(K) = COMP(J,K) 

400 CONTINUE 

WRITE(3 , 130) J ,DUM,REL1(J) ,REL2(J) 

500 CONTINUE 

WRITE(3 , 140) 

WRITE(3 , 150) RS , ELMAX(NREP) , LMU(NREP) .RSLOW, LEVEL 

10 FORMAT (IX, '0UT1.DAT : Output File of the RETP1 simulation') 

20 FORMAT ( IX , ' after ',15,' replications’,/) 

25 FORMAT( IX, 'COMMENTS : 8 COMPONENTS IN SERIES ') 

26 FORMAT (IX, ' DF = 2 * (NFC - NCOMP) ' ,/) 

30 FORMAT ( IX, ' Input Parameters:',/) 

40 FORMAT (IX, ' I SEED NCOMP ALPHA TOL NCS TCN' ,/) 

50 FORMAT (IX, F10.1,I8,F8.4,F8.5,2I6,/) 

60 FORMAT (IX, ' I TY(I) PARMl(I) PARM2(I) UT(I) NC(I) NF(I)',/) 
70 F0RMAT(1X, 13,16, 2F9 . 5 , F8 . 2 , 218) 

80 FORMAT (IX,/, 'Output Parameters for the LAST Replication:',/) 

90 FORMAT (IX, ' INF(I) TT(I) ELM(I) ER(I)', 

* ' EBETA(I) ' ,/) 

100 FORMAT (IX, 13 , 16 , E16 . 7 , 2F14 . 7 , F14 .7) 

110 F0RMAT(1X,/, 'Cut -Set Data: ' ,/) 

120 FORMAT (IX, ' J NUM Component List 

* ' REL1 REL2(M) ' ,/) 

130 FORMAT (IX, 13, 15, 913, 2 F12. 9) 

140 FORMAT ( IX , / , ' RS ELMAX(M) LMU(M) ' , 
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RSLOW 



LEVEL' ,/) 



* 1 
150 FORMAT (IX, 5F12 . 7 , / ) 

C 

CLOSE (UNIT-3) 

C 

RETURN 
END 

C 

C This portion of the file contains functions and subroutines 
C used in the RETP1 model. 

C - 14 May 91 

C - by Yee Kah-Chee SMC 2802 

C 

C A. Random Number Generating Subroutine (LRNDPC) . 

C (Courtesy of Mr. David Lim Hung-Heng) 

C 

SUBROUTINE LRNDPC (DSEED,U,N) 

INTEGER N, I 

REAL* 8 U(N) 

REAL* 8 D31M1 , DSEED, D31 

D31M1=2**31 - 1 
D31 =2**31 

DATA D31M1/2147483647 . DO/ 

DATA D31 /2147483648 . DO/ 

DO 5 1=1, N 

DSEED = DM0D( 950706376 . D0*DSEED , D31M1) 

DSEED = DM0D( 16807 . D0*DSEED, D31M1) 

5 U(I) = DSEED / D31 
RETURN 
END 



B. Survivability Function. 



FUNCTION SURV (TYPE, PARI, PAR2, UTIL) 

C 

C This function returns the survival probability of the component of 
C different types (TYPE) with scale (PARI) and shape (PAR2) parameters 
C given the specified utilization times or cycles (UTIL) . 

C 

INTEGER TYPE, N 
REAL* 8 PARI, PAR2 , UTIL 
C 

IF (TYPE. EQ. 1) THEN 

SURV = EXP( - (PAR1*UTIL) ) 

ELSEIF (TYPE.EQ. 2) THEN 

SURV = EXP ( - ( ( PAR1*UTIL) **PAR2 ) ) 

ELSE 

N = NINT(UTIL) 

SURV = PAR1**N 
END IF 
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c 



END 



C. Bubble Sort Routine in ASCENDING Order. 



SUBROUTINE BUBBLE (N , LIST , OLIST) 

This subroutine performs a bubble sort in increasing order (ie. sink 
the greater numeral) for the first N terms in an array LIST and 
returns the result in OLIST. 



LOGICAL DONE 
INTEGER N, K, PAIR 
REAL* 8 LIST(*) , OLIST(*) 



Sink the larger of the pair. 

DO 50 K = 1, N 

OLIST(K) = LIST (K) 

50 CONTINUE 

PAIR - N - 1 
DONE - .FALSE. 

DO WHILE (.NOT. DONE) 

DONE = .TRUE. 

DO 100 K - 1, PAIR 

IF (OLIST(K) . GT. OLIST (K+l)) THEN 
TEMP = OLIST (K) 

OLIST (K) - OLIST (K+l) 
0LIST(K+1) - TEMP 
DONE - .FALSE. 

END IF 

100 CONTINUE 

PAIR = PAIR - 1 
ENDDO 
END 



D. Unbiasing Factor for Biased MLE for Weibull Shape Parameter. 



FUNCTION BN(I) 

C 

C This functon returns the value of the unbiased factor for the biased 
C maximum likelihood estimate of the shape parameter of a Weibull 
C distribution with a sample size of N. 

C 

INTEGER I 
IF (I.LE.5) THEN 

BN - (1*0 . 699)/(5 . 0) 

ELSEIF (I.EQ.6) THEN 
BN - 0.752 

ELSEIF (I.EQ.7) THEN 



53 



non 



BN - 0.786 

ELS El F (I.EQ.8) THEN 
BN = 0.82 

ELSEIF (I.EQ.9) THEN 
BN = 0.8395 
ELSEIF (I.EQ.10) THEN 
BN - 0.859 

ELSEIF (I.EQ.ll) THEN 
BN - 0.871 

ELSEIF (I.EQ.12) THEN 
BN - 0.883 

ELSEIF (I.EQ.13) THEN 
BN - 0.892 

ELSEIF (I.EQ.14) THEN 
BN - 0.901 

ELSEIF (I.EQ.15) THEN 
BN = 0.9075 
ELSEIF (I.EQ.16) THEN 
BN - 0.914 

ELSEIF (I.EQ.17) THEN 
BN = 0.9185 
ELSEIF (I.EQ.18) THEN 
BN = 0.923 

ELSEIF (I.EQ.19) THEN 
BN - 0.927 

ELSEIF (I.EQ.20) THEN 
BN - 0.931 

ELSEIF (I.LE.25) THEN 

BN = 0.931+(I-20)*0. 014/5.0 
ELSEIF (I.LE.30) THEN 

BN = 0.945+(I-25)*0. 01/5.0 
ELSEIF (I.LE.40) THEN 

BN = 0.955+(I-30)*0. 011/10.0 
ELSEIF (I.LE.60) THEN 

BN = 0.966+(I-40)*0. 012/20.0 
ELSEIF (I.LE.80) THEN 

BN = 0.978+(I-60)*0. 006/20.0 
ELSEIF (I.LE.100) THEN 

BN = 0.984+(I-80)*0. 003/20.0 
ELSEIF (I.LE.120) THEN 

BN = 0.987+(I-100)*0. 003/20.0 
ELSE 

BN - 1.0 
END IF 
RETURN 
END 



E. Biased MLE of Weibull Shape Parameter. 



SUBROUTINE MLESHAPE (T , N , R , DEL , B , BNEW) 
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c 

C This subroutine returns a biased estimator (BNEW) for a Weibull 
C shape parameter using the Newton-Raphson' s Method of Successive 
C Approximation. The data parameters consist of an ascending ordered 
C list of failure times (T) , sample size (N) , number of failed samples 
C (R) , tolerance for convergence (DEL) and an initial estimate of the 
C shape parameter (B) . 

C 

LOGICAL DONE 
INTEGER N, R, I 

REAL* 8 GFUNCT, G PRIME, B, BOLD, BNEW, T(*) , DEL, 

* TERM1 , TERM2 , TERM3 , SUM1, SUM2 , SUM3 , SUM4 , STEP 
C 

BNEW - B 
DONE = .FALSE. 

C 

DO WHILE (.NOT. DONE) 

C 

DONE = .TRUE. 

TERM1 - FLOAT (N-R)*(T(R)**BNEW) 

TERM2 = FLOAT (N-R)*(T(R)**BNEW)*LOG(T(R) ) 

TERM3 - FLOAT (N - R) * (T (R) **BNEW) *LOG ( T(R) )*LOG(T(R) ) 

SUM1 = 0.0 

SUM2 =0.0 

SUM3 = 0.0 

SUM4 - 0.0 

C 

DO 50 I = 1, R 

SUM1 - SUM1 + T ( I ) **BNEW 

SUM2 = SUM2 + (T(I)**BNEW)*L0G(T(I) ) 

SUM3 = SUM3 + (T(I)**BNEW)*LOG(T(I) )*LOG(T(I) ) 

SUM4 = SUM4 + L0G(T(I)) 

50 CONTINUE 
C 

GFUNCT = ( SUM2+TERM2 ) / ( SUM1+TERM1 ) - (1.0/BNEW) 

* - (1 . 0/FLOAT(R) )*SUM4 

C 

GPRIME = ( 1 . 0/ ( SUMl+TERMl ) **2 ) * ( (SUM1+TERM1)*(SUM3+TERM3) 

* - ( SUM2+TERM2 ) **2 ) 

* + (1.0/BNEW** 2) 

C 

C PRINT 60, GFUNCT, GPRIME, BNEW 

C 60 FORMAT (IX, ' GFUNCT =', F8 . 3 , ' GPRIME =', F8 . 3 , ' BNEW«',F8.3) 

C 

C Control magnitude of the marching step towards convergence 
C as no more than 0.1. 

C 

IF ( (GFUNCT. LT.0) .AND. (GPRIME. GT. 0) ) THEN 
STEP = VMAX(- .1D0, ( GFUNCT/G PRIME ) ) 

ELSEIF ( (GFUNCT. GT.0) .AND. (GPRIME. LT.0)) THEN 
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STEP = VMAX( - . IDO , (GFUNCT/GPRIME) ) 

ELSE 

STEP - VMIN( .IDO, (GFUNCT/GPRIME)) 

ENDIF 

BOLD = BNEW 
BNEW - BNEW - STEP 

Check for convergence of the MLE for the shape parameter B. 



IF (ABS (BOLD -BNEW) .GT.DEL) THEN 
DONE - .FALSE. 

ENDIF 

Avoid overflow error due to large MLE value caused by small 
GPRIME (slope) as GFUNCT approaches to near zero. 

Stop when magnitude of BNEW exceeds 7. 

IF (BNEW. GT. 7.0) THEN 
BNEW = BOLD 
DONE - .TRUE. 

ENDIF 

ENDDO 

RETURN 

END 



E. Chi-Square Quantile Function. 



C 

C 

C 

C 

C 

C 

C 

C 



FUNCTION CHISQD(P.N) 



Modified version of Algorithm 451 from Comunications of the ACM 
Aug 1977 Vol.16 No. 8 . 

This function evaluates the quantile at the probability level P 
(left tail area) for the Chi-square distribution with 
N degrees of freedom. 

REAL*8 P 
REAL X 
INTEGER IF 

DIMENSION C(21) , A(19) 

DATA C/ 1 . 565326E- 3 , 

* 1 . 060438E- 3 , 

* -6 . 950356E-3 , 

* -1 . 323293E-2 , 

* 2 . 277679E-2 , 

* -8 . 986007E-3, 

* - 1 . 513904E- 2 , 

* 2 . 530010E-3 , 
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-1.450117E-3, 


* 


5 . 169654E-3 , 


* 


- 1 . 153761E-2 , 


* 


1 . 128186E-2 , 


* 


2 . 607083E-2 , 


* 


-0.2237368, 


* 


9 . 780499E- 5 , 


* 


-8.426812E-4, 


* 


3 . 125580E-3 , 


* 


-8 . 553069E-3 , 


* 


1 . 348028E-4 , 


* 


0.4713941, 


* 


1.0000886 / 


DATA 


A/ 1 . 264616E-2 , 


* 


-1.425296E-2, 


* 


1 . 400483E-2 , 


* 


- 5 . 886090E-3 , 


* 


- 1 . 091214E- 2 , 


* 


-2 . 304527E-2 , 


* 


3 . 135411E-3 , 


* 


-2 . 728484E-4 , 


* 


-9 . 699681E-3 , 


* 


1 . 316872E-2 , 


* 


2 . 618914E- 2 , 


* 


-0.2222222, 


* 


5.406674E-5, 


* 


3 . 483789E-5 , 


* 


-7 . 274761E-4 , 


* 


3 . 292181E-3 , 


* 


-8.729713E-3, 


* 


0.4714045, 


* 


1. / 



IF (N-2) 10, 20, 30 



10 CALL XFROMP( . 5*(1 . -P) ,X, IF) 

CHISQD = X 

CHISQD - CHI SQD*CHI SQD 
RETURN 

20 CHISQD « -2.*LOG(l. -P) 

RETURN 
30 F = N 

FI = l./F 

CALL XFROMP ( P , X , I F) 

T - X 

F2 = SQRT(F1)*T 

IF (N.GE. (2+INT(4.*ABS(T)))) GO TO 40 

CHISQD = (((((( (C(1)*F2+C(2) )*F2+C(3) )*F2+C(4) )*F2 

* +C(5))*F2+C(6))*F2+C(7))*F1+((((((C(8)+C(9)*F2)*F2 

* +C ( 10 ) ) *F2+C ( 11 ) ) *F2+C ( 12 ) ) *F2+C ( 13 ) ) *F2+C ( 14) ) ) *F1+ 

* ( ( ( ( (C(15)*F2+C(16) )*F2+C(17) )*F2+C(18) )*F2 

* +C(19) )*F2+C(20) )*F2+C(21) 
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GO TO 50 

40 CHISQD=(((A(1)+A(2)*F2)*F1+(((A(3)+A(4)*F2)*F2 

* +A(5))*F2+A(6)))*F1+(((((A(7)+A(8)*F2)*F2+A(9))*F2 

* +A(10))*F2+A(11))*F2+A(12)))*F1+(((((A(13)*F2 

* +A ( 14 ) ) *F2+A ( 15 ) ) *F2+A ( 16 ) ) *F2+A ( 17 ) ) *F2*F2 

* +A(18))*F2+A(19) 

50 CHISQD - CHISQD*CHISQD*CHISQD*F 
RETURN 
END 



C 

C 

C 

C 

C 

C 



C 

C 

C 

C 



F. Standard Normal Variate Computation Subroutine. 

SUBROUTINE XFROMP ( P , X , I FAULT ) 

Algorithm AS 24 J .R.STAT.SOC. C. (1969) Vol.18. No. 3. 

This subroutine computes the standard normal deviate X for 
the specified left tail area P. 

REAL* 8 P 
DIMENSION A(5) 

DIMENSION CONNOR (17), HSTNGS(6) 

DATA CONNOR/ 8 . 0327350124E- 17 , 



* 


1 . 4483264644E-15 , 


* 


2.4668270103E-14, 


* 


3 . 9554295164E-13 , 


* 


5 . 9477940136E-12 , 


* 


8 . 3507027951E-11 , 


* 


1.0892221037E-9, 


* 


1 . 3122532964E-8 , 


* 


1.4503852223E-7, 


* 


1.4589169001E-6, 


* 


1 . 3227513228E- 5 , 


* 


1 . 0683760684E-4 , 


* 


7 . 5757575758E-4 , 


* 


4 . 6296296296E-3 , 


* 


2 . 3809523810E-2 , 


* 


0.1, 


★ 


0.33333333333 / 



DATA RTHFPI / 1.2533141373 / 
DATA RRT2PI / 0.3989422804 / 
DATA TERMIN / 1.0E-11 / 

DATA HSTNGS / 2.515517, 

* 0.802853, 

* 0.010328, 

* 1.432788, 
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c 



* 

* 



0.189269, 
0.001308 / 



I FAULT = 1 

IF ((P.LE.0.0) .OR. (P.GE.1.0)) GO TO 100 
I FAULT - 0 

Get first approximation XO to deviate by Hastings' formula 
B - P 

IF(B.GT.0.5) B - 1.0 - B 

F - - LOG(B) 

E - SQRT(F+F) 

XO - -E + ( (HSTNGS (3 )*E+HSTNGS (2) )*E+HSTNGS (1) ) / 

* ( ( (HSTNGS (6)*E+HSTNGS (5) )*E+HSTNGS (4) )*E+1 . 0) 

IF (XO.LT.0.0) GO TO 1 
XO = 0.0 
PO = 0.5 
XI = -RTHFPI 
GO TO 7 



Find the area PO corresponding to XO 

1 Y - X0**2 

IF (XO . LE. -1.9) GO TO 3 
Y = -0 . 5*Y 

(1) series approximation 

PO = CONNOR (1) 

DO 2 L=2 , 17 

2 PO - P0*Y + CONNOR (L) 

PO - ( P0*Y+1 . 0 ) *X0 
XI = - (PO+RTHFPI)*EXP( -Y) 

PO = PO*RRT2PI +0.5 
GO TO 7 

(2) continued fraction approximation 

3 Z - 1.0/Y 
A(2) = 1.0 
A(3) - 1.0 
A(4) - Z + 1.0 
A(5) - 1.0 
W = 2.0 

C 

4 DO 6 L=1 ,3,2 
DO 5 J-1,2 

K = L + J 
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KA - 7 - K 
C 

5 A(K) - A(KA) + A(K)*W*Z 
C 

6 W = W + 1.0 
APPRXU = A(2)/A(3) 

APPRXL - A(5)/A(4) 

C “ APPRXU - APPRXL 
IF (C.GE.TERMIN) GO TO 4 
XI = APPRXL/XO 
PO - -X1*RRT2PI*EXP( -0 . 5*Y) 

C 

C Get accurate value of deviate by Taylor Series 
C (XI, X2 , X3 are derivatives for the Taylor Series 
C 

7 D = F + LOG ( PO ) 

X2 = X0*X1*X1 -XI 

X3 = Xl**3 + 2 . 0*XO*X1*X2 -X2 
X = ( (X3*D/3 . 0+X2)*D/2 . 0+Xl)*D + XO 
IF (P.LE.0.5) GO TO 100 
X = -X 
100 RETURN 
END 



G. Maximum Function. 



FUNCTION VMAX(X , Y) 
REAL* 8 X, Y 
IF(X.GT.Y) THEN 

VMAX - X 
ELSE 

VMAX = Y 
END IF 
RETURN 
END 



H. Minimum Function. 



FUNCTION VMIN(X,Y) 

REAL* 8 X, Y 
IF(X.LT.Y) THEN 
VMIN - X 
ELSE 

VMIN - Y 
END IF 
RETURN 
END 

4. Program Output. (OUTl.DAT) 
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The result for the simulation run based on the input parameters specified in 
IN1.DAT are computed and written to the file OUTl.DAT. A sample of this file is 
as follows. 



OUTl.DAT : Output File of the RETP1 simulation 
after 1000 replications 



COMMENTS : 8 COMPONENTS IN SERIES 

DF - 2 * (NFC - NCOMP ) 

Input Parameters: 

ISEED NCOMP ALPHA TOL NCS TCN 

16807.0 8 0.2000 0.01000 8 3 



I 


TY(I) 


PARMl(I) PARM2(I) 


UT(I) NC(I) 


NF(I) 




1 


1 


0.00200 1.00000 


10.00 


15 


3 




2 


1 


0.00400 1.00000 


10.00 


15 


3 




3 


1 


0.00600 1.00000 


10.00 


15 


3 




4 


1 


0.00800 1.00000 


10.00 


15 


3 




5 


2 


0.00200 2.00000 


10.00 


15 


3 




6 


2 


0.00400 2.00000 


10.00 


15 


3 




7 


2 


0.00600 2.00000 


10.00 


15 


3 




8 


2 


0.00800 2.00000 


10.00 


15 


3 




Output Parameters for the LAST Replication: 






i 


NF(I) 


TT(I) 


ELM(I) 




ER(I) 


EBETA(I) 


i 


3 


0.19 6909 1E+04 


0.0010157 




0.1703872 


0.0000000 


2 


3 


0 . 3355078E+03 


0.0059611 




1.0000000 


0.0000000 


3 


3 


0 . 5999981E+03 


0.0033333 




0.5591815 


0.0000000 


4 


3 


0 . 6423594E+03 


0.0031135 




0.5223055 


0.0000000 


5 


3 


0 . 2563836E+16 


0.0000000 




0.0000000 


6.3524983 


6 


3 


0 . 8186629E+07 


0.0000004 




0.0000615 


2.9571536 


7 


3 


0 . 5675037E+09 


0.0000000 




0.0000009 


4.2418658 


8 


3 


0 . 5675519E+05 


0.0000529 




0.0088672 


2.3471863 



Cut -Set Data: 



J NUM Component List 



REL1 REL2 (M) 



1 

2 

3 

4 

5 

6 



110 0 0 
1 2 0 0 0 

1 3 0 0 0 

1 4 0 0 0 

1 5 0 0 0 

1 6 0 0 0 



0 0 0 0 
0 0 0 0 
0 0 0 0 
0 0 0 0 
0 0 0 0 
0 0 0 0 



0 0.980198622 
0 0.960789382 
0 0.941764534 
0 0.923116326 
0 0.999600053 
0 0.998401225 



0.990279973 

0.944286704 

0.968452990 

0.970502377 

0.999999940 

0.999680758 



61 



7 

8 



1 7 0 0 0 0 

1 8 0 0 0 0 



0000 0.996406436 0.999911249 
0000 0.993620396 0.988757312 



RS 


ELMAX(M) 


LMU(M) 


RSLOW 


0.8089645 


0.0059611 


0.0057325 


0.8691733 



LEVEL 

0.3380000 
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APPENDIX C : Users’ Guide for RETP2 



Reliability Estimation Test Plan 2 (RETP2). 
by YEE, Kah-Chee 
July 91 



1. Brief Description. 

RETP2 is a computer program written in FORTRAN that runs on the Amdahl 
mainframe at NPGS. It allows the user to simulate exponential and Weibull failure 
times of component items being tested to evaluate the accuracy of a confidence limit 
estimation procedure based on Type I data censoring (that is, testing items of 
component i until a specified total testing time is achieved for each of them). 

2. Program Input. (IN2.DATA) 

The input of the program are specified to the program via an input file called 
IN2.DAT. A sample input file is shown below. 

This file contains the inputs required by the RETP2 model. 

Update only the numerical values between dotted lines as appropriate. 

Do not delete any of the comment lines. (IN2.DAT) 20 Jun 91 



Value 


Type 


Units 


Description 


Variable 


16807.0 


Real 


- 


initial random seed 


ISEED 


8 


Int 


- 


total # of components in system 


NCOMP 


4 


INT 


- 


# OF EXPONENTIAL COMPONENTS 


NEXP 


4 


INT 


- 


# OF WEIBULL COMPONENTS 


NWEI 


0 


INT 


- 


# OF GEOMETRIC COMPONENTS 


NGEO 


0.01 


Real 


- 


tolerance for MLE 


TOL 


0.20 


REAL 


- 


DESIRED SIGNIFICANCE LEVEL 


ALPHA 


1000 


INT 


- 


# OF REPLICATIONS DESIRED 


NREP 


3 


INT 




TEST CASE NUMBER 

1 = all EXP 

2 - all WEI 

3 - EXP + WEI 

4 = EXP + WEI + CYC 


TCN 


8 


Int 


- 


number of cut sets 


NCS 
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TEST PLAN : Testing until TT(I) (total test time) is accumulated. 
(Use REAL numbers ONLY!!!) 



Comp 

Number 

I 

Int 


Comp 

Type 

TY(I) 

Int 


Comp Parameters 
Scale Shape 

PARMl(I) PARM2 ( I ) 

Real Real 


Util 

Time/Cycle 

UT(I) 

(hrs) (eyes) 


Test Plan Inputs 
Accumulated 
TT(I) NC(I) 

(hrs) (eyes) Int 


1.0 


1.0 


0.005 


1.0 




5.0 




5400.0 


20.0 


2.0 


1.0 


0.005 


1.0 




5.0 




5400.0 


20.0 


3.0 


1.0 


0.005 


1.0 




5.0 




5400.0 


20.0 


4.0 


1.0 


0.005 


1.0 




5.0 




5400.0 


20.0 


5.0 


2.0 


0.010 


2.0 




15.0 




2700.0 


20.0 


6.0 


2.0 


0.010 


2.0 




15.0 




2700.0 


20.0 


7.0 


2.0 


0.010 


2.0 




15.0 




2700.0 


20.0 


8.0 


2.0 


0.010 


2.0 




15.0 




2700.0 


20.0 


Note : 


TY(I) 


-1 EXPONENTIAL P(surv) 


“ exp( 


-PARM2)*T) 






TY(I) 


=2 WEI BULL P(surv) 


- exp( 


- ( PARM1*T ) **PARM2 ) 




TY(I) 


=3 GEOMETRIC P(surv) 


= PARM1**T 




SYSTEM 


CONFIGURATION : 


Identification 


of CUT 


SETs 










- min groups 


of 


components that have to fail 








for the system to fail 






Cut Set 


# in Set 


List of Components in 


l Cutset 




J 


COMP (J ,1) 


C0MP(J , 2) . 




up to C0MP(J,1) components 


i 




1 


10 0 0 


0 


0 0 


0 


0 0 




2 




1 


2 0 0 0 


0 


0 0 


0 


0 0 




3 




1 


3 0 0 0 


0 


0 0 


0 


0 0 




4 




1 


4 0 0 0 


0 


0 0 


0 


0 0 




5 




1 


5 0 0 0 


0 


0 0 


0 


0 0 




6 




1 


6 0 0 0 


0 


0 0 


0 


0 0 




7 




1 


7 0 0 0 


0 


0 0 


0 


0 0 




8 




1 


8 0 0 0 


0 


0 0 


0 


0 0 





3. Program Flow and Logic. (NAME2.DEF, PARM2.DEF and RETP2.FOR) 

Input parameters are first read in by the program by calling the INPUT 
subroutine. The program then evoke the SIM subroutine which generates the 
random failure times and compute the key statistics required in the procedure. The 
next subroutine EVAL determines the measures of accuracy for run. REPORT is 
the subroutine which generates the output file for the run OUT2.DAT. 

The variables in the program RETP2.FOR are described in the file 
NAME2.DEF as listed below. 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 



This file contains the declaration for input and output variables 
used in the the RETP2 model. (NAME2.DEF) 20 Jun 91 



Input Variables. 

ISEED « initial random seed selected. 

SEED = current random seed. 

RS * true overall series system reliability. 

ALPHA = level of significance desired. 

NREP - number of replications desired for the simulation. 

TPN = test plan number (1). 

TCN ■= test case number (1, 2, 3 or 4) . 

NCOMP - total number of components in the system. 

NEXP * number of components with EXP failure times. 

NWEI * number of components with WEI failure times. 

NGEO = number of components with GEO failure times. 

TOL - desired tolerance for MLE of iJEI shape parameter. 

Distribution: Exponential WEIbull GEOmetric 

TY(I) = type: 1 23 

PARM(1,I): Scale(l/hr) Scale(l/hr) Prob 

PARM(2 , I) : - Shape 

UT(I) = utilization time (hrs) for component i (EXP and WEI) . 
UC(I) ■= utilization cycles for component i (GEO only). 

NC(I) -= number of test samples (sample size) for cmponent i. 
NF(I) - desired number of failures in test for component i. 

NCS «= number of cut- sets for the system. 

C0MP(J,K) = kth parameter of cut-set j (first being the no. of 
components belonging to the cut-set) . 



Assumed Variables. 



MAXCOMP = maximum number of components allowed in the system. 
MAXREP - maximum number of replications permitted. 

MAXCUT = maximum number of cut-sets. 

Program and Output Variables. 



RS = true overall system reliability. 

TT(I) = total accumulated failure time (hr) for component i 
(EXP and WEI only) . 

TC(I) = total accumulated cycles to failure (incl. failure cycle) 
for component i (GEO only) . 

EBETA(I) *= estimate for shape parameter of component i (if Weibull) . 
REL1(J) = actual reliability for cut-set j. 

REL2(J) « computed reliability for cut-set j for current replication. 
ELM(I) = estimated component failure rate (1/hrs) for component i. 
ELMAX(M) «= max estimated component failure rate for rep m (1/hrs) . 
ER(I) = ratio of estimated failure rate to LMAX. 

ET(I) = same as TT(I) except that these are for Weibull components. 
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C NFC(M) = total number of failed test components. 

C LMU(M) = upper confidence bound for failure rate (1/hrs), 

C RSL(M) - lower confidence limit estimated for system reliability 

C ... for the mth replication. 

C ORSL(M) = ordered RSL(M) (ascending). 

C RSLOW = (l-ALPHA)xlOO percentile of set of RSL(M) . 

C LEVEL = achieved confidence level, ie. proportion of RSL(M) that 
C are lesser than RS (conservative estimate) . 

C 

C-- END OF NAME2.DEF 

C 



Together with the main program in RETP2.FOR are the other subroutines 
needed in the simulation. The declaration of variables is done in the file 
PARM2.DEF. Relevant descriptions are included as comment lines in the source 
code to help explain the program segments. A listing of PARM2.DEF and 
RETP2.FOR is given below. 

c 

C This file contains the declaration for input and output variables 
C used in the the RETP2 model. (PARM2.DEF) 20 Jun 91 

C 

INTEGER MAXCOMP, MAXREP 

PARAMETER ( MAXCOMP - 100 , MAXREP - 1000 , MAXCUT - 20 ) 

REAL* 8 I SEED, SEED 

INTEGER NREP , TCN, NCOMP, NEXP, NWEI , NGEO, NCS , 

* NC (MAXCOMP) , NF (MAXCOMP) , TY (MAXCOMP ) , NFC (MAXREP) , 

* UC (MAXCOMP) , TC (MAXCOMP) , COMP (MAXCUT .MAXCOMP) 

REAL* 8 RS, ALPHA, UT (MAXCOMP) , TT (MAXCOMP) , 

* PARM (2, MAXCOMP ) , ELM (MAXCOMP ) , ER(MAXCOMP) , ET (MAXCOMP) , 

* LMU (MAXREP) , RSL(MAXREP) , ORSL(MAXREP) , 

* ELMAX (MAXREP ) , RSLOW, LEVEL, TOL, EBETA(MAXCOMP) , 

* REL1 (MAXCUT ) , REL2 (MAXCUT) 

C 

C0MM0N/BL0CK1/ISEED , SEED, NREP, TCN, NCOMP, NC, NF, NEXP, NWEI, 

* NGEO, NCS, TY, NFC, UC, TC, COMP 
COMMON/BLOCK2/RS , ALPHA, UT, TT, PARM, ELM, ER, ET, LMU, 

* RSL, ORSL, ELMAX, RSLOW, LEVEL, TOL, EBETA, 

* REL1 , REL2 
C 

C-- END OF PARM2.DEF 

C 

C 

C This file contains the main program and the subroutines 
C for the Reliability Estimation Test Plan 2 (RETP2) model. 

C (RETP2 . FOR) - runs on a IBM PC Compatible . 

C 

C IBM Mainframe Version. 

C 
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C Test Plan 2 : Testing until accumulated time or cycles is achieved 

C - for component i. 

C 

C by Yee Kah-Chee SMC 2802. 

C 20 Jun 91. 

C 

C 1. Main Program (RETP2) . 

C 

PROGRAM RETP2 

Include the declaration files . 

INCLUDE 'NAME 2 DEF' 

INCLUDE ' PARM2 DEF' 

Read in input data. 

CALL INPUT 

Activate simulation . 

CALL SIM 

Process and evaluate output data. 

CALL EVAL 

Generate simulation report. 

CALL REPORT 

STOP 
END 



2. Input Initialisation Subroutine (INPUT). 

SUBROUTINE INPUT 

This subroutine reads in the inputs for the RETP2 model. 

Include the declaration file. 

INCLUDE ' PARM2 DEF' 

INTEGER I, J, K, DUM2(11) 

REAL* 8 DUM1 ( 7 ) 

Read data from 'IN2.DAT' designated as logic unit 1. 
0PEN(UNIT=1) 
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c 



READ(l.lO) 

10 FORMAT (IX ,/////) 

READ(1 , *) I SEED 
READ(1,*) NCOMP 
READ(1 ,*) NEXP 
READ(1 , *) NWEI 
READ(1 , *) NGEO 
READ(1 ,*) TOL 
READ(1,*) ALPHA 
READ(1 ,*) NREP 
READ(1 ,*) TCN 
READ(1 , 20) 

20 FORMAT (IX,///) 

READ ( 1 , * ) NCS 
C 

READ(1 , 30) 

30 FORMAT (IX,////////) 

C 

DO 50 K - 1, NCOMP 
READ(1 ,*) DUM1 
C 

I - NINT(DUM1(1) ) 

TY(I) - NINT(DUM1(2) ) 

C 

IF (TY(I).EQ.l) THEN 
PARM (1,1) = DUM1 ( 3 ) 
PARM(2 , I) = DUMl (4) 
UT(I) = DUMl ( 5 ) 

TT(I) - DUMl ( 6 ) 

NC(I) - NINT(DUMl(7)) 
EBETA(I) - 0 

ELSEIF ( T Y ( I ) . EQ . 2 ) THEN 
PARM (1,1) - DUMl ( 3 ) 
PARM(2 , I) - DUMl (4) 
UT(I) - DUMl (5) 

TT(I) - DUMl ( 6 ) 

NC(I) - NINT ( DUMl ( 7 ) ) 
ELSEIF (TY(I) .EQ.3) THEN 
PARM (1,1) - DUMl (3) 
PARM(2 , I) = DUMl ( 4 ) 
UC(I) = NINT(DUM1(5)) 
UT(I) - DUMl ( 5 ) 

TT(I) = DUMl ( 6 ) 

NC(I) - NINT (DUMl (7 ) ) 
END IF 

50 CONTINUE 
C 

READ(1 , 60) 

60 FORMAT (IX, ///////////) 
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DO 80 I - 1, NCS 
READ(1 , *) DUH2 
J - DUM2 ( 1 ) 

C0MP(J , 1) - DUM2 ( 2 ) 

DO 70 K = 1, C0MP(J , 1) 
C0MP(J,K+1) = DUM2 (K+2) 
70 CONTINUE 
80 CONTINUE 

CLOSE (UNIT-1) 

RETURN 

END 



3. Subroutine (SIM) 



SUBROUTINE SIM 

This subroutine simulates NREP possible outcomes of the test plan 
desired in order to obtain the raw estimates of LMU(M) and RSL(M) 
for each of the replication. 

Include the declaration file 
and declare local variables. 

INCLUDE ' PARM2 DEF' 

INTEGER I, J, K, M, ISUM, KEY, ICOUNT 

INTEGER NCYC(MAXCOMP) , NSYS 

REAL* 8 UNI 

REAL* 8 SUM, PROD, 

* FT (MAXCOMP , 200) , 0FT(MAXC0MP , 200) , 

* DFT(200) , D0FT(200) 

LOGICAL FLAG 



SEED = ISEED 



Compute overall true system reliability RS. 

RS - 1.0 
DO 30 J = 1, NCS 
PROD -1.0 

DO 20 I = 1, C0MP(J , 1) 

K - C0MP(J , 1+1) 

PROD = PROD* ( 1 - SURV(TY(K) , PARM(1 , K) , PARM(2 , K) ,UT (K) ) ) 
20 CONTINUE 

REL1(J) - 1.0 - PROD 
RS - RS * REL1(J) 

30 CONTINUE 
C 

C Start of Simulation. 

C (Initialize replication counter M) . 

C 



69 



o o 



M = 1 

DO WHILE (M.LE.NREP) 

C 

C Test Plan : Sample and determine unknown NF(I) 

C (with known NC(I), Weibull case) until TT(I) is reached. 

C 

C Generate NC(I) failure times, put them in ascending order 
C with the smallest failure time on the top of the list. 

C 

DO 60 I = 1, NCOMP 
C 

C Exponential components. Test each component until it fails, check 
C to see if TT(I) is exceeded, if not, carry on testing. 

C 

I F ( TY ( I ) . EQ . 1 ) THEN 
SUM =0.0 
L = 1 

DO WHILE ( SUM . LE . TT ( I ) ) 

CALL LRNDPC( SEED, UNI, 1) 

FT (I ,L) = -L0G(UNI)/PARM(1 , I) 

SUM = SUM + FT(I,L) 

L = L + 1 
ENDDO 

NF(I) = L - 2 
C 

C Weibull components. Generate NC(I) failure times, put them in 
C an ascending order with the smallest failure time on top and 
C determine NF(I) . 

C 

ELSEIF(TY(I) . EQ. 2) THEN 
DO 40 K = 1, NC(I) 

CALL LRNDPC (SEED, UNI, 1) 

FT(I,K) = (1 .0/PARM(l , I) )*( -LOG(UNI) )**(1 . 0/PARM(2 , I) ) 
DFT(K) = FT(I,K) 

40 CONTINUE 

CALL BUBBLE(NC(I) ,DFT,D0FT) 

DO 42 K - 1, NC(I) 

OFT (I ,K) = DOFT(K) 

42 CONTINUE 

SUM =0.0 
NF(I) = 0 
DO 45 K - 1, NC(I) 

IF (OFT (I ,K) . LT. TT(I) ) THEN 
NF(I) = NF(I) + 1 
END IF 

45 CONTINUE 



C 



GEOMETRIC COMPONENTS. DETERMINE NF(I). 
ELSEIF(TY(I) . EQ. 3) THEN 
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NF ( I ) - 0 

DO 50 K - 1, TT(I) 

CALL LRNDPC( SEED, UNI, 1) 

IF (UNI . GT . PARM( 1 , I) ) THEN 
NF(I) = NF(I) + 1 
END IF 

50 CONTINUE 

C 

END IF 

60 CONTINUE 
C 

C Determine total number of failed test components as well as 
C checking for zero component failure or just failures from a 
C SINGLE component (FLAG will be set to .TRUE, if so). 

C 

I SUM = 0 
ICOUNT - 0 
DO 70 I = 1, NCOMP 
IF (NF(I).GT.O) THEN 
ICOUNT - ICOUNT + 1 
END IF 

I SUM = I SUM + NF(I) 

70 CONTINUE 
C 

NFC(M) = I SUM 
C 

IF ( ICOUNT. LE.l) THEN 
FLAG = .TRUE. 

ELSE 

FLAG = .FALSE. 

ENDIF 



Case A : More than ONE component type experienced failures 
in the test. 

IF (.NOT. FLAG) THEN 



Estimate the failure rate of each component. 

DO 90 I - 1, NCOMP 

Exponential components . 

IF (TY(I).EQ.l) THEN 
ET(I) - TT( I) 

ELM(I) - FLOAT (NF(I))/TT(I) 



Weibull components. 

ELSEIF (TY(I) . EQ. 2) THEN 
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DO 75 K = 1, NC ( I ) 

DOFT(K) - OFT(I ,K) 

75 CONTINUE 

CALL MLESHAPE(DOFT,NC(I) , IMAX(1 ,NF(I) ) ,TOL, 1 . DO , EBETA(I) ) 
EBETA(I) - BN(NC(I))*EBETA(I) 

SUM =0.0 

DO 80 K = 1, NF(I) 

SUM = SUM + OFT ( I , K) **EBETA ( I ) 

80 CONTINUE 

ET(I) = FLOAT (NC (I) -NF(I) )*TT(I)**EBETA(I) + SUM 
ELM(I) = FLOAT (NF ( I ) ) /ET ( I ) 

Geometric components. 

ELSEIF (TY(I).EQ.3) THEN 
ET(I) = TT(I) 

ELM(I) = FLOAT(NF(I))/TT(I) 

END IF 

90 CONTINUE 



Determine the maximum failure rate estimate 
and identify that component. 

ELMAX(M) =0.0 
KEY = 0 

DO 100 I = 1, NCOMP 

IF ( ELM ( I ) . GT . ELMAX (M) ) THEN 
ELMAX(M) = ELM(I) 

KEY = I 
ENDIF 

100 CONTINUE 



Compute the ratios of the failure rate estimates to their maximum. 

DO 110 I = 1, NCOMP 

ER(I) = ELM(I) / ELMAX(M) 

110 CONTINUE 



Determination of LMU(M) 

SUM = 0.0 

DO 120 I = 1, NCOMP 

SUM = SUM + (ER(I)*ET(I) ) 

120 CONTINUE 

LMU(M) = CHISQD( 1-ALPHA, NINT(1 . 3*2*(1+NFC(M) ) ) )/ (2*SUM) 



Compute estimate of overall reliability RSL(M) for the system. 
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RSL(M) =1.0 
DO 140 J = 1, NCS 
PROD =1.0 

DO 130 I = 1, COMP(J.l) 

K = COMP( J , 1+1) 

IF (TY(K).EQ.l) THEN 

PROD = PROD* ( 1 - SURV(TY (K) , LMU (M)*ER(K) , EBETA(K) ,UT(K) ) ) 
ELSEIF (TY(K) . EQ . 2) THEN 

PROD = PROD*(l - SURV(TY(K) , (LMU(M)*ER(K) )**(1 . /EBETA(K) ) , 
* EBETA(K) ,UT(K) ) ) 

ELSEIF (TY(K) . EQ . 3) THEN 

PROD = PROD*(l - SURV (TY(K) , 1 . DO-LMU(M)*ER(K) , 0 . DO ,UT(K) ) ) 
END IF 

130 CONTINUE 

REL2(J) = 1.0 - PROD 
RSL(M) = RSL(M) * REL2(J) 

140 CONTINUE 
C 

ENDIF 

C 

C Case B : Where there are at most ONE component which experienced 

C failures during the test. 

C 

IF (FLAG) THEN 
C 

C Determine number of complete systems (NSYS) implied by the test 
C based on cut- set information after first determining the number 
C of mission cycles tested for each component (to the nearest 
C integer) (NCYC(I)). 

C 

I SUM = 0 

DO 150 I = 1, NC0MP 

NCYC(I) = INT(TT(I)/UT(I)) 

ISUM = I SUM + NCYC(I) 

ELM(I) = 0 
ER(I) = 0 
EBETA(I) =1.0 
150 CONTINUE 

ELMAX(M) = 0 
LMU(M) = 0 
C 

NSYS = ISUM 
ISUM = 0 

DO 170 J = 1, NCS 
ISUM = 0 

DO 160 I = 1, COMP ( J , 1 ) 

ISUM = ISUM + NCYC(COMP(J ,1+1) ) 

160 CONTINUE 

ISUM = ISUM / C0MP(J , 1) 



73 



non non oooooooooo non o non o non non 



NSYS = IMIN(NSYS , ISUM) 

170 CONTINUE 

For zero failures in the test. 

IF ( I COUNT . EQ . 0) THEN 

RSL(M) = ALPHA** (1 . O/FLOAT (NSYS) ) 

For failures experienced by a particular component type. 

ELSE 

CALL GETP (NSYS, ALPHA, RSL(M)) 

ENDIF 

END IF 



Increment replication counter. 
M = M + 1 



ENDDO 

END 



4. Subroutine (EVAL) . 



SUBROUTINE EVAL 

This subroutine calls BUBBLE to sort the array RSL(NREP) in 
ascending order to get an ordered array ORSL(NREP) . It also 
determine the estimate for RSLOW at the specified significance 
level ALPHA and the value of LEVEL in which ORSL(LEVEL) is closest 
to the true reliability RS . 



Include the declaration files 
and declare the local variables. 

INCLUDE ' PARM2 DEF' 

INTEGER INDEX 
REAL* 8 DIFF 

Order the array RSL(NREP) in ascending order. 

DO 10 M - 1, NREP 
ORSL(M) - RSL(M) 

10 CONTINUE 

Bubble Sort. Sink the larger of the pair. 
CALL BUBBLE (NREP, RSL,ORSL) 
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C Determine the (1-ALPHA)% lower confidence bound for the system 
C reliability. 

C 

RSLOW = ORSL(NINT(NREP*( 1-ALPHA) ) ) 

C 

C Finding the % confidence level for the true reliability RS . 

C (ie. the proportion of RSL(M) lesser than RS) 

C 

DIFF - 1.0 
INDEX - 0 

DO 200 M = 1, NREP 

IF (ABS (ORSL(M) -RS) .LT.DIFF) THEN 
DIFF - ABS(0RSL(M) -RS) 

INDEX - M 
END IF 

200 CONTINUE 
C 

LEVEL = FLOAT (INDEX) /NREP 



Record evaluated parameters in RAW2.DAT (unit 2). 

0PEN(UNIT=2) 

WRITE(2 , 300) 

300 FORMAT (IX, ' M LMU(M) ELMAX(M) RSL(M) ' , 

* ' ORSL(M) NFC(M) ' ) 

DO 500 M = 1, NREP 

WRITE(2 ,400) M,LMU(M) ,ELMAX(M) ,RSL(M) ,0RSL(M) ,NFC(M) 
400 FORMAT (IX, I6,2F12.7,2F12.7,I10) 

500 CONTINUE 

CLOSE (UNIT-2) 

RETURN 

END 



5. Report Generation Subroutine (REPORT). 



SUBROUTINE REPORT 

This subroutine record the simulation results into the '0UT2.DAT' 
file as logic unit 3. 

Include the declaration files 
and declare local variables. 

INCLUDE ' PARM2 DEF' 

INTEGER I, J, K, DUM(10) 

Write to output file '0UT2.DAT' designated as logic unit 3. 
0PEN(UNIT=3) 
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WRITE (3 , 10) 

WRITE(3 , 20) NREP 
WRITE(3 , 25) 

WRITE (3 , 26) 

WRITE(3 , 30) 

WRITE (3 ,40) 

WRITE (3,50) ISEED , NCOMP , ALPHA , TOL, NCS , TCN 
C 

WRITE(3 , 60) 

DO 200 I = 1, NCOMP 

WRITE(3 , 70) I , TY ( I ) , PARM (1,1) , PARM (2,1) ,UT(I) ,TT(I) ,NC(I) ,NF(I) 
200 CONTINUE 
WRITE(3 , 80) 

WRITE(3 , 90) 

DO 300 I - 1, NCOMP 

WRITE(3 , 100) I,NF(I) ,ET(I) ,ELM(I) ,ER(I) , EBETA(I) 

300 CONTINUE 

WRITE(3 , 110) 

WRITE(3 , 120) 

DO 500 J = 1, NCS 
DO 400 K = 1, 10 
DUM(K) - COMP( J , K) 

400 CONTINUE 

WRITE(3 , 130) J , DUM.RELl (J) ,REL2(J) 

500 CONTINUE 

WRITE(3 , 140) 

WRITE(3 , 150) RS , ELMAX(NREP) , LMU(NREP) .RSLOW, LEVEL 



10 


FORMAT ( IX , ' OUT2 . DAT 


: Output File 


of the RETP2 simulation') 


20 


FORMAT (IX, ' 


after ' , 15 , 


' replications',/) 




25 


FORMAT ( IX , ' COMMENTS 


: 8 COMPONENT 


IN SERIES 


') 


26 


FORMAT (IX, ' 


DF = NINT 


(1.3 * 2 * (1 + NFC)) 


',/) 



30 FORMAT ( IX, ' Input Parameters:',/) 

40 FORMAT (IX, ' ISEED NCOMP ALPHA TOL NCS TCN',/) 

50 FORMAT ( IX, F10 .1,I8,F8.4,F8.5,2I6,/) 

60 FORMAT (IX, ' I TY(I) PARMl(I) PARM2(I) UT(I) TT(I) NC(I)', 

* ' NF(I) ' ,/) 

70 F0RMAT(1X, 13,16, 2F9 . 5 , 2F9 . 2 , 218) 

80 FORMAT ( IX ,/,' Output Parameters for the LAST Replication:',/) 

90 FORMAT (IX, ' INF(I) ET(I) ELM(I) ER(I)', 

* ' EBETA(I) ' ,/) 

100 FORMAT ( IX , 13 , 16 , E16 , 7 , 2F14 . 7 , F14 . 7 ) 

110 F0RMAT(1X,/, 'Cut-Set Data: ' ,/) 

120 F0RMAT(1X, ' J NUM Component List ' , 

* ' REL1 REL2(M) ' , /) 

130 F0RMAT(1X, 13,15, 913 , 2F12.9) 

140 FORMAT (IX,/, ' RS ELMAX(M) LMU(M) ' , 

* ' RSLOW LEVEL',/) 

150 FORMAT ( IX, 5F12 . 7 , /) 
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c 

CLOSE (UNIT-3) 

C 

RETURN 

END 

C 

C This portion of the file contains functions and subroutines 
C used in the RETP2 model. 

C - 20 Jun 91 

C - by Yee Kah-Chee SMC 2802 

C 

C A. Random Number Generating Subroutine (LRNDPC) . 

C (Courtesy of Mr. David Lim Hung-Heng) 

C 

SUBROUTINE LRNDPC (DSEED, U,N) 

INTEGER N , I 

REAL* 8 U(N) 

REAL*8 D31M1 , DSEED , D31 

D31M1— 2**31 - 1 
D31 -2**31 

DATA D31M1/2147483647 . DO/ 

DATA D31 /2147483648 . DO/ 

DO 5 1=1 ,N 

DSEED = DMOD(950706376.D0*DSEED,D31Ml) 

DSEED = DM0D( 16807 . D0*DSEED , D31M1) 

5 U(I) = DSEED / D31 
RETURN 
END 



B. Survivability Function. 



FUNCTION SURV (TYPE, PARI, PAR2, UTIL) 

C 

C This function returns the survival probability of the component of 
C different types (TYPE) with scale (PARI) and shape (PAR2) parameters 
C given the specified utilization times or cycles (UTIL) . 

C 

INTEGER TYPE, N 
REAL*8 PARI, PAR2 , UTIL 
C 

IF (TYPE.EQ.l) THEN 

SURV = EXP ( - ( PAR1*UTIL) ) 

ELSEIF (TYPE. EQ. 2) THEN 

SURV = EXP( - ( (PAR1*UTIL)**PAR2) ) 

ELSE 

N - NINT(UTIL) 

SURV - PAR1**N 
END IF 
C 

END 
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C. Bubble Sort Routine in ASCENDING Order. 



SUBROUTINE BUBBLE (N , LIST , OLIST) 



This subroutine performs a bubble sort in increasing order (ie. sink 
the greater numeral) for the first N terms in an array LIST and 
returns the result in OLIST. 



LOGICAL DONE 
INTEGER N, K, PAIR 
REAL* 8 LIST(*) , OLIST(*) 



Sink the larger of the pair. 

DO 50 K - 1, N 

OLIST(K) - LIST (K) 

50 CONTINUE 

PAIR = N - 1 
DONE = .FALSE. 

DO WHILE (.NOT. DONE) 

DONE = .TRUE. 

DO 100 K - 1, PAIR 

IF (OLIST (K) .GT.0LIST(K+1)) THEN 
TEMP = OLIST (K) 

OLIST(K) - 0LIST(K+1) 
0LIST(K+1) = TEMP 
DONE - .FALSE. 

END IF 

100 CONTINUE 

PAIR - PAIR - 1 
ENDDO 
END 



D. Unbiasing Factor for Biased MLE for Weibull Shape Parameter. 



FUNCTION BN(I) 

C 

C This functon returns the value of the unbiased factor for the biased 
C maximum likelihood estimate of the shape parameter of a Weibull 
C distribution with a sample size of N. 

C 

INTEGER I 
IF (I.LE.5) THEN 

BN = (1*0 . 699)/(5 . 0) 

ELSEIF (I.EQ.6) THEN 
BN = 0.752 

ELSEIF (I.EQ.7) THEN 
BN = 0.786 

ELSEIF (I.EQ.8) THEN 
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BN = 


0.82 




ELSEIF 


(I.EQ.9) 


THEN 


BN - 


0.8395 




ELSEIF 


(I.EQ.10) 


THEN 


BN = 


0.859 




ELSEIF 


(I.EQ.ll) 


THEN 


BN = 


0.871 




ELSEIF 


(I.EQ.12) 


THEN 


BN = 


0.883 




ELSEIF 


(I.EQ.13) 


THEN 


BN = 


0.892 




ELSEIF 


(I.EQ.14) 


THEN 


BN = 


0.901 




ELSEIF 


(I.EQ.15) 


THEN 


BN = 


0.9075 




ELSEIF 


(I.EQ.16) 


THEN 


BN - 


0.914 




ELSEIF 


(I.EQ.17) 


THEN 


BN = 


0.9185 




ELSEIF 


(I.EQ.18) 


THEN 


BN = 


0.923 




ELSEIF 


(I.EQ.19) 


THEN 


BN - 


0.927 




ELSEIF 


(I.EQ.20) 


THEN 


BN = 


0.931 




ELSEIF 


(I.LE.25) 


THEN 


BN = 


0. 931+(I- 


20)*0. 014/5.0 


ELSEIF 


(I.LE.30) 


THEN 


BN = 


0.945+(I- 


25)*0. 01/5.0 


ELSEIF 


(I.LE.40) 


THEN 


BN = 


0. 955+(I- 


30)*0. 011/10.0 


ELSEIF 


(I.LE.60) 


THEN 


BN = 


0 . 966+(I - 


40)*0. 012/20.0 


ELSEIF 


(I.LE.80) 


THEN 


BN = 


0. 978+(I- 


60)*0. 006/20.0 


ELSEIF 


(I.LE.100) THEN 


BN - 


0.984+(I- 


80)*0. 003/20.0 


ELSEIF 


(I.LE.120) THEN 


BN = 


0 . 987+(I - 


100)*0. 003/20.0 


ELSE 






BN - 


1.0 




ENDIF 






RETURN 






END 







E. Biased MLE of Weibull Shape Parameter. 



SUBROUTINE MLESHAPE (T ,N,R, DEL , B , BNEW) 

This subroutine returns a biased estimator (BNEW) for a Weibull 
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c 

c 

c 

c 

c 



c 



c 

c 



c 



c 



c 



c 

c 

c 

c 

c 

c 

c 



shape parameter using the Newton-Raphson' s Method of Successive 
Approximation. The data parameters consist of an ascending ordered 
list of failure times (T) , sample size (N) , number of failed samples 
(R) , tolerance for convergence (DEL) and an initial estimate of the 
shape parameter (B) . 

LOGICAL DONE 

INTEGER N, R, I 

REAL* 8 GFUNCT , GPRIME , B, BOLD, BNEW, T(*) , DEL, 

* TERM1, TERM2, TERM 3 , SUM1, SUM2 , SUM3 , SUM4 , STEP 

BNEW - B 

DONE = .FALSE. 



DO WHILE (.NOT. DONE) 

DONE = .TRUE. 

TERM1 = FLOAT (N-R)*(T(R)**BNEW) 

TERM2 = FLOAT (N-R)*(T (R)**BNEW)*LOG(T (R) ) 

TERM3 = FLOAT(N-R)*(T(R)**BNEW)*LOG(T(R) )*LOG(T(R) ) 

SUMl - 0.0 

SUM2 =0.0 

SUM3 =0.0 

SUM4 =0.0 



DO 50 I = 1, R 

SUMl = SUMl + T(I)**BNEW 

SUM2 = SUM2 + (T(I )**BNEW)*LOG(T(I ) ) 

SUM3 = SUM3 + (T(I)**BNEW)*LOG(T(I) )*LOG(T(I) ) 
SUM4 = SUM4 + L0G(T(I)) 

50 CONTINUE 



GFUNCT 

* 



( SUM2+TERM2 ) / ( SUM1+TERM1 ) - (1.0/BNEW) 
- ( 1 . 0/FLOAT (R) ) *SUM4 



GPRIME = (1 . 0/(SUMl+TERMl)**2)*( (SUMl+TERMl)*(SUM3+TERM3) 

* - (SUM2+TERM2)**2 ) 

* + (1.0/BNEW** 2) 



PRINT 60, GFUNCT, GPRIME, BNEW 

60 FORMAT ( IX , ' GFUNCT =', F8 . 3 , ' GPRIME =', F8 . 3 , ' BNEW=',F8.3) 



Control magnitude of the marching step towards convergence 
as no more than 0.1. 

IF ( (GFUNCT. LT.0) .AND. (GPRIME. GT. 0) ) THEN 
STEP = VMAX( - . 1D0 , (GFUNCT/GPRIME) ) 

ELS El F ( (GFUNCT. GT.0) .AND. (GPRIME .LT. 0) ) THEN 
STEP = VMAX(- .1D0, (GFUNCT/GPRIME)) 

ELSE 
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STEP = VMIN ( . IDO , (GFUNCT/G PRIME) ) 

ENDIF 

C 

BOLD = BNEW 
BNEW - BNEW - STEP 
C 

Check for convergence of the MLE for the shape parameter B. 



IF (ABS(BOLD-BNEW) .GT.DEL) THEN 
DONE = .FALSE. 

ENDIF 

Avoid overflow error due to large MLE value caused by small 
GPRIME (slope) as GFUNCT approaches to near zero. 

STOP WHEN MAGNITUDE OF BNEW EXCEEDS 7. 

IF (BNEW. GT. 7.0) THEN 
BNEW - BOLD 
DONE = .TRUE. 

ENDIF 

ENDDO 

RETURN 

END 



E. Chi-Square Quantile Function. 



FUNCTION CHISQD(P.N) 



Modified version of Algorithm 451 from Comunications of the ACM 
Aug 1977 Vol.16 No. 8 . 

This function evaluates the quantile at the probability level P 
(left tail area) for the Chi-square distribution with 
N degrees of freedom. 



REAL*8 P 
REAL X 
INTEGER IF 

DIMENSION C(2l) , A(19) 
DATA C/ 1 . 565326E-3 , 

* 1 . 060438E-3 , 

* -6 . 950356E- 3 , 

* -1. 323293E-2, 

* 2 . 277679E- 2 , 

* -8 . 986007E- 3 , 

* -1 . 513904E-2 , 

* 2 . 530010E- 3 , 

* - 1 . 450117E- 3 , 

* 5 . 169654E- 3 , 
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* 


- 1 . 153761E- 2 , 


* 


1 . 128186E-2 , 


* 


2 . 607083E- 2 , 


* 


-0.2237368, 


k 


9 . 780499E-5 , 


k 


-8 . 426812E-4 , 


k 


3 . 125580E-3 , 


k 


-8 . 553069E-3 , 


* 


1 . 348028E-4 , 


* 


0.4713941, 


* 


1.0000886 / 


DATA A/ 


1 . 264616E- 2 , 


* 


- 1 . 425296E- 2 , 


* 


1 . 400483E-2 , 


* 


- 5 . 886090E-3 , 


* 


-1.091214E-2, 


* 


-2 . 304527E-2 , 


* 


3.135411E-3, 


* 


- 2 . 728484E-4 , 


* 


-9 . 699681E-3 , 


* 


1 . 316872E-2 , 


* 


2 . 618914E- 2 , 


* 


-0.2222222, 


* 


5 . 406674E-5 , 


* 


3.483789E-5, 


k 


- 7 . 274761E-4 , 


* 


3 . 292181E-3 , 


k 


- 8 . 729713E-3 , 


k 


0.4714045, 


k 


1. / 



IF (N-2) 10, 20, 30 



10 CALL XFR0MP( .5*(1. -P) ,X,IF) 

CHISQD = X 

CHISQD - CHI SQD*CHI SQD 
RETURN 

20 CHISQD - -2 . *LOG(l . - P) 

RETURN 
30 F = N 

FI = l./F 

CALL XFROMP(P,X, IF) 

T = X 

F2 - SQRT(F1)*T 

IF (N.GE. (2+INT(4.*ABS(T)))) GO TO 40 

CHISQD ■= (((((((C(1)*F2+C(2))*F2+C(3))*F2+C(4))*F2 

* +C(5))*F2+C(6))*F2+C(7))*F1+((((((C(8)+C(9)*F2)*F2 

* +C(10) )*F2+C(11) ) *F2+C ( 1 2 ) ) *F2+C ( 13 ) ) *F2+C ( 14 ) ) ) *F 1+ 

* (((((C(15)*F2+C(16))*F2+C(17))*F2+C(18))*F2 

* +C(19))*F2+C(20))*F2+C(21) 

GO TO 50 

40 CHISQD=(((A(1)+A(2)*F2)*F1+(((A(3)+A(4)*F2)*F2 
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* +A(5))*F2+A(6)))*F1+(((((A(7)+A(8)*F2)*F2+A(9))*F2 

* +A(10))*F2+A(11))*F2+A(12)))*F1+(((((A(13)*F2 

* +A(14))*F2+A(15))*F2+A(16))*F2+A(17))*F2*F2 

* +A(18) )*F2+A(19) 

50 CHISQD - CHISQD*CHISQD*CHISQD*F 
RETURN 
END 

F. Standard Normal Variate Computation Subroutine. 

SUBROUTINE XFROMP (P , X , I FAULT) 

Algorithm AS 24 J .R. STAT. SOC . C. (1969) Vol.18. No. 3. 

This subroutine computes the standard normal deviate X for 
the specified left tail area P. 

REAL* 8 P 
DIMENSION A(5) 

DIMENSION CONNOR (17), HSTNGS(6) 

DATA CONNOR/ 8 . 0327350124E- 17 , 



C 

C 

C 

C 

C 

C 



C 

C 

C 

C 



* 


1.4483264644E-15, 


* 


2.4668270103E-14, 


* 


3 . 9554295164E-13 , 


* 


5 . 9477940136E-12 , 


* 


8 . 3507027951E- 11 , 


* 


1 . 089222 1037E- 9 , 


* 


1 . 3122532964E-8 , 


* 


1.4503852223E-7, 


* 


1.4589169001E-6, 


* 


1 . 3227513228E-5 , 


* 


1 . 0683760684E-4 , 


* 


7 . 5757575758E-4 , 


* 


4 . 6296296296E- 3 , 


* 


2 . 3809523810E- 2 , 


* 


0.1, 


* 


0.33333333333 / 



DATA RTHFPI / 1.2533141373 / 
DATA RRT2PI / 0.3989422804 / 
DATA TERMIN / 1.0E-11 / 

DATA HSTNGS / 2.515517, 

* 0.802853, 

* 0.010328, 

* 1.432788, 

* 0.189269, 

* 0.001308 / 
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I FAULT = 1 

IF ((P.LE.0.0) .OR. (P.GE.1.0)) GO TO 100 
I FAULT - 0 

Get first approximation XO to deviate by Hastings' formula 
B = P 

IF(B . GT . 0 . 5) B = 1.0 - B 

F = - LOG(B) 

E - SQRT(F+F) 

XO - -E + ( (HSTNGS(3)*E+HSTNGS(2) )*E+HSTNGS (1) )/ 

* ( ( (HSTNGS ( 6 ) *E+HSTNGS ( 5 ) ) *E+HSTNGS (4 ) ) *E+1 . 0 ) 

IF (XO.LT.O.O) GO TO 1 
XO - 0.0 
PO - 0.5 
XI = -RTHFPI 
GO TO 7 



Find the area PO corresponding to XO 

1 Y - X0**2 

IF (X0.LE.-1.9) GO TO 3 
Y - - 0 . 5*Y 

(1) series approximation 

PO - CONNOR(l) 

DO 2 L=2 , 17 

2 PO - PO*Y + CONNOR (L) 

P0 = (PO*Y+1.0)*XO 
XI = - (PO+RTHFPI)*EXP( -Y) 

PO = P0*RRT2PI +0.5 
GO TO 7 

(2) continued fraction approximation 

3 Z = 1.0/Y 
A(2) = 1.0 
A(3) - 1.0 
A(4) = Z + 1.0 
A(5) = 1.0 
W - 2.0 

C 

4 DO 6 L=1 ,3,2 
DO 5 J-1,2 

K = L + J 
KA = 7 - K 
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5 A(K) - A(KA) + A(K)*W*Z 
C 

6 W - W + 1.0 
APPRXU - A(2)/A(3) 

APPRXL - A(5)/A(4) 

C - APPRXU - APPRXL 
IF (C.GE.TERMIN) GO TO 4 
XI = APPRXL/XO 
PO - -X1*RRT2PI*EXP( -0 . 5*Y) 

C 

C Get accurate value of deviate by Taylor Series 
C (XI, X2, X3 are derivatives for the Taylor Series 
C 

7 D - F + LOG(PO) 

X2 - X0*X1*X1 -XI 

X3 - Xl**3 + 2 . 0*XO*X1*X2 -X2 
X - ( (X3*D/3 . 0+X2)*D/2 . 0+Xl)*D + XO 
IF (P.LE.0.5) GO TO 100 
X = -X 
100 RETURN 
END 



G. MAXIMUM FUNCTIONS. 



FUNCTION VMAX(X, Y) 
REAL* 8 X, Y 
IF (X.GT.Y) THEN 
VMAX - X 
ELSE 

VMAX - Y 
END IF 
RETURN 
END 

FUNCTION IMAX(X, Y) 
INTEGER X, Y 
IF (X.GT.Y) THEN 
IMAX - X 
ELSE 

IMAX - Y 
ENDIF 
RETURN 
END 



H. Minimum Functions. 



FUNCTION VMIN(X.Y) 
REAL* 8 X, Y 
IF (X.LT.Y) THEN 
VMIN = X 
ELSE 
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VMIN = Y 
END IF 
RETURN 
END 

FUNCTION IMIN(I , J) 
INTEGER I, J 
IF (I.LT.J) THEN 
IMIN = I 
ELSE 

IMIN = J 
END IF 
RETURN 
END 



I. ROUTINE TO EVALUATE VALUE OF P. 



SUBROUTINE GETP(N, ALPHA, NEWP) 

INTEGER N 

REAL* 8 ALPHA, OLDP, NEWP, TOL 
LOGICAL DONE 

OLDP - ALPHA** (1.0/FLOAT(N)) 

NEWP =1.0 
TOL = 0.0001 
DONE = .FALSE. 

DO WHILE (.NOT. DONE) 

GFUNCT = N*OLDP**(N-l) - (N-l)*OLDP**N - ALPHA 
GPRIME = N*(N-l)*OLDP**(N-2) - N*(N-l)*OLDP**(N-l) 

NEWP = OLDP - (GFUNCT/GPRIME) 

IF ( (ABS (NEWP-OLDP) . LE. TOL) .OR. (ABS (GFUNCT) . LE . TOL) ) THEN 
DONE = .TRUE. 

END IF 

OLDP = NEWP 
END DO 
RETURN 
END 

4. Program Output. (OUT2.DAT) 

The result for the simulation run based on the input parameters specified in 
IN2.DAT are computed and written to the file OUT2.DAT. A sample of this file is 
as follows. 

OUT2.DAT : Output File of the RETP2 simulation 
after 1000 replications 

COMMENTS : 8 COMPONENT IN SERIES 

DF = NINT (1.3 * 2 * (1 + NFC)) 
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Input Parameters: 





ISEED 


NCOMP 


ALPHA 


TOL 


NCS TCN 






16807.0 


8 


0.2000 0, 


.01000 


8 3 






I 


TY(I) 


PARMl(I) 


PARM2 ( I ) 


UT(I) 


TT(I) 


NC(I) 


NF(I) 


1 


1 


0.00500 


1.00000 


5.00 


5400.00 


20 


36 


2 


1 


0.00500 


1.00000 


5.00 


5400.00 


20 


23 


3 


1 


0.00500 


1.00000 


5.00 


5400.00 


20 


28 


4 


1 


0.00500 


1.00000 


5.00 


5400.00 


20 


25 


5 


2 


0.01000 


2.00000 


15.00 


2700.00 


20 


20 


6 


2 


0.01000 


2.00000 


15.00 


2700.00 


20 


20 


7 


2 


0.01000 


2.00000 


15.00 


2700.00 


20 


20 


8 


2 


0.01000 


2.00000 


15.00 


2700.00 


20 


20 



Output Parameters for the LAST Replication: 



I 


NF(I) 


ET(I) 


ELM(I) 


ER(I) 


EBETA(I) 


1 


36 


0 . 5400000E+04 


0.0066667 


1.0000000 


0.0000000 


2 


23 


0 . 5400000E+04 


0.0042593 


0.6388889 


0.0000000 


3 


28 


0 . 5400000E+04 


0.0051852 


0.7777778 


0.0000000 


4 


25 


0 . 5400000E+04 


0.0046296 


0.6944444 


0.0000000 


5 


20 


0 . 2378985E+05 


0.0008407 


0.1261042 


1.5441284 


6 


20 


0 . 4825114E+06 


0.0000414 


0.0062175 


2.1482593 


7 


20 


0 . 2539978E+05 


0.0007874 


0.1181112 


1.5843514 


8 


20 


0 . 3708000E+07 


0.0000054 


0.0008091 


2.5186896 



Cut-Set Data: 



J 


NUM 




Component List 








RELl 


REL2 (M) 


1 


1 


1 


0 


0 


0 


0 


0 


0 


0 


0 


0.975309908 


0.955163479 


2 


1 


2 


0 


0 


0 


0 


0 


0 


0 


0 


0.975309908 


0.971117675 


3 


1 


3 


0 


0 


0 


0 


0 


0 


0 


0 


0.975309908 


0.964950144 


4 


1 


4 


0 


0 


0 


0 


0 


0 


0 


0 


0.975309908 


0.968645990 


5 


1 


5 


0 


0 


0 


0 


0 


0 


0 


0 


0.977751195 


0.927053154 


6 


1 


6 


0 


0 


0 


0 


0 


0 


0 


0 


0.977751195 


0.981007099 


7 


1 


7 


0 


0 


0 


0 


0 


0 


0 


0 


0.977751195 


0.923940539 


8 


1 


8 


0 


0 


0 


0 


0 


0 


0 


0 


0.977751195 


0.993218899 






RS 




ELMAX(M) 




LMU(M) 


RSLOW 


LEVEL 


0 


.8269590 


0.0066667 




0.0091745 0.7754698 0.9940000 
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APPENDIX D : Evaluation of Subroutines and Functions 



RANDOM NUMBER GENERATOR 
(LRNDPC) 

Evaluation 

One thousand uniform random real numbers between 0 and 1 are generated 
using the random number generating routine LRNDPC. From these uniformly 
distributed numbers, 1000 exponential (with scale parameter 1) numbers and 1000 
Weibull (with scale parameter 1 and shape parameter 2) numbers were generated. 

Uniform Random Variate 



Plot of U(N+1)vsU(N) 

for 1000 Uniform RVs from LRNDPC 




0 0.1 02 0.3 0.4 0.5 0.6 0.7 0.8 0.8 1 



U(N) 

Figure 1 : Uniform RVs generated by LRNDPC 



Figure 1 above shows a plot of 1000 uniform real numbers against their 
predecessors. The uniformity of the distribution of points over the state space 
confirms LRNDPC’s adequacy in generating uniform random numbers. 
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Exponential and Weibull real numbers were generated using these 1000 
Uniform(0,l) random numbers. The cumulative histograms of these resultant 
random variates were compared with their respective theoretical cumulative 
distribution functions ( cdfs ). 



Exponential Random Variate 



Distribution Fitting 

tot 1000 Exp(1) RV* by LRNDPC 




Figure 2 : Exponential RVs generated by LRNDPC 



Figure 2 shows the close distribution fit between the theoretical cdf (line) and 
the cumulative distribution of exponential RV generated using LRNDPC. 

F(t) = exp (~\t) 
t = -hn{F(t)) 
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Weibull Random Variate 



Distribution Fitting 

(of 1000 Wel(1,2) RVs by LRNDPC 
(Je. Seal© - 1 , Shape - 2) 




Figure 3 : Weibull RVs generated by LRNDPC 



Figure 3 shows the close distribution fit between the theoretical cdf (line) and 
the cumulative distribution of Weibull RV generated using LRNDPC. 

F(t) = exp{-(Ar) p } 

I = -l[ln(F(())P 
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Plot of Unbiasing Factor B(N) vs N 
for Weibull Shape Parameter Estimation 

N = Test Sample Size 
B(N) = Unbiasing Factor for MLE 



Plot of BN vs N 

(Weibull Shape Parameter Estimation) 




The function BN(N) returns the linear-interpolated values of the unbiasing 
factor for the raw MLE $ for both RETP1 and RETP2. 



Evaluation of Subroutine CHISQD and XFROMP 

The x 2 statistics for 1 to 499 degrees of freedom for a values of 0.1 and 0.2 are 
generated using the routines CHISQD and XFROMP. These outputs matched those 
tabulated in the mathematical tables of any general textbook on statistics. 
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Prob [ Xd/ < table value ] = 1 - a =0.9 



df 




0 




1 




2 




3 




4 




5 




6 




7 




8 




9 


0 






2 . 


71 


4 . 


61 


6 . 


25 


7 . 


78 


9 . 


24 


10 . 


64 


12 . 


02 


13 . 


36 


14 . 


68 


1 


15 . 


99 


17 . 


27 


18 . 


55 


19 . 


81 


21 . 


06 


22 . 


31 


23 . 


54 


24 . 


77 


25 . 


99 


27 . 


20 


2 


28 . 


41 


29 . 


62 


30 . 


81 


32 . 


01 


33 . 


20 


34 . 


38 


35 . 


56 


36 . 


74 


37 . 


92 


39 . 


09 


3 


40 . 


26 


41 . 


42 


42 . 


58 


43 . 


75 


44 . 


90 


46 . 


06 


47 . 


21 


48 . 


36 


49 . 


51 


50 . 


66 


4 


51 . 


81 


52 . 


95 


54 . 


09 


55 . 


23 


56 . 


37 


57 . 


51 


58 . 


64 


59 . 


77 


60 . 


91 


62 . 


04 


5 


63 . 


17 


64 . 


30 


65 . 


42 


66 . 


55 


67 . 


67 


68 . 


80 


69 . 


92 


71 . 


04 


72 . 


16 


73 . 


28 


6 


74 . 


40 


75 . 


51 


76 . 


63 


77 . 


75 


78 . 


86 


79 . 


97 


81 . 


09 


82 . 


20 


83 . 


31 


84 . 


42 


7 


85 . 


53 


86 . 


64 


87 . 


74 


88 . 


85 


89 . 


96 


91 . 


06 


92 . 


17 


93 . 


27 


94 . 


37 


95 . 


48 


8 


96 . 


58 


97 . 


68 


98 . 


78 


99 . 


88 


100 . 


98 


102 . 


08 


103 . 


18 


104 . 


28 


105 . 


37 


106 . 


47 


9 


107 . 


57 


108 . 


66 


109 . 


76 


110 . 


85 


111 . 


94 


113 . 


04 


114 . 


13 


115 . 


22 


116 . 


32 


117 . 


41 


10 


118 . 


,50 


119 . 


,59 


120 . 


68 


121 . 


77 


122 . 


,86 


123 . 


95 


125 . 


,04 


126 . 


,12 


127 . 


,21 


128 . 


.30 


11 


129 . 


39 


130 . 


47 


131 . 


56 


132 . 


,64 


133 . 


.73 


134 . 


,81 


135 . 


,90 


136 . 


98 


138 . 


,07 


139 . 


,15 


12 


140 , 


,23 


141 . 


32 


142 . 


40 


143 . 


,48 


144 . 


,56 


145 . 


,64 


146 . 


.72 


147 . 


80 


148 . 


,89 


149 . 


.97 


13 


151 . 


,05 


152 . 


12 


153 . 


,20 


154 , 


,28 


155 . 


,36 


156 . 


,44 


157 . 


.52 


158 . 


,60 


159 , 


.67 


160 . 


.75 


14 


161 . 


.83 


162 . 


,90 


163 . 


.98 


165 . 


,06 


166 . 


,13 


167 . 


,21 


168 , 


.28 


169 . 


,36 


170 , 


.43 


171 . 


,51 


15 


172 . 


.58 


173 . 


.66 


174 . 


.73 


175 . 


.80 


176 . 


.88 


177 , 


.95 


179 , 


.02 


180 , 


.09 


181 , 


.17 


182 , 


.24 


16 


183 , 


.31 


184 . 


,38 


185 . 


.45 


186 . 


.52 


187 , 


.60 


188 , 


.67 


189 


.74 


190 , 


.81 


191 , 


.88 


192 , 


.95 


17 


194 . 


.02 


195 . 


.09 


196 . 


, 16 


197 , 


.23 


198 , 


.29 


199 , 


.36 


200 


.43 


201 , 


.50 


202 , 


.57 


203 , 


.64 


18 


204 , 


.70 


205 , 


.77 


206 , 


.84 


207 . 


.91 


208 , 


.97 


210 , 


.04 


211 


.11 


212 


.17 


213 , 


.24 


214 , 


.31 


19 


215 . 


.37 


216 . 


.44 


217 , 


.50 


218 , 


.57 


219 , 


.63 


220 , 


.70 


221 


.76 


222 


.83 


223 


.89 


224 


.96 


20 


226 


.02 


227 , 


.09 


228 


.15 


229 , 


.21 


230 


.28 


231 , 


.34 


232 


.40 


233 


.47 


234 


.53 


235 


.59 


21 


236 


.65 


237 . 


.72 


238 , 


.78 


239 , 


.84 


240 


.90 


241 , 


.97 


243 


.03 


244 


.09 


245 


.15 


246 


.21 


22 


247 


.27 


248 , 


.33 


249 


.40 


250 , 


.46 


251 


.52 


252 , 


.58 


253 


. 64 


254 


.70 


255 


.76 


256 


.82 


23 


257 


.88 


258 . 


.94 


260 


.00 


261 


.06 


262 


.12 


263 


.18 


264 


.24 


265 


.29 


266 


.35 


267 


.41 


24 


268 


.47 


269 , 


.53 


270 


.59 


271 


.65 


272 


.70 


273 


.76 


274 


.82 


275 


.88 


276 


.94 


277 


.99 


25 


279 


.05 


280 , 


.11 


281 


.16 


282 


.22 


283 


.28 


284 


.34 


285 


.39 


286 


.45 


287 


.51 


288 


.56 


26 


289 


.62 


290 , 


.67 


291 


.73 


292 


.79 


293 


.84 


294 


.90 


295 


.95 


297 


.01 


298 


.07 


299 


.12 


27 


300 


.18 


301 , 


.23 


302 


.29 


303 


.34 


304 


.40 


305 


.45 


306 


.51 


307 , 


.56 


308 


.61 


309 


.67 


28 


310 


.72 


311 , 


.78 


312 


.83 


313 


.89 


314 


.94 


315 


.99 


317 


.05 


318 


.10 


319 


.15 


320 


.21 


29 


321 


.26 


322 


.31 


323 


.37 


324 


.42 


325 


.47 


326 


.53 


327 


.58 


328 


.63 


329 


.68 


330 


.74 


30 


331 


.79 


332 , 


.84 


333 


.89 


334 


.95 


336 


.00 


337 


.05 


338 


.10 


339 


.15 


340 


.20 


341 


.26 


31 


342 


.31 


343 


.36 


344 


.41 


345 


.46 


346 


.51 


347 , 


.56 


348 


.62 


349 


.67 


350 


.72 


351 


.77 


32 


352 


.82 


353 


.87 


354 


.92 


355 


.97 


357 


.02 


358 


.07 


359 


.12 


360 


.17 


361 


.22 


362 


.27 


33 


363 


.32 


364 , 


.37 


365 


.42 


366 


.47 


367 


.52 


368 


.57 


369 


.62 


370 


.67 


371 


.72 


372 


.77 


34 


373 


.82 


374 


.87 


375 


.92 


376 


.96 


378 


.01 


379 


.06 


380 


.11 


381 


.16 


382 


.21 


383 


.26 


35 


384 


.31 


385 


.35 


386 


.40 


387 


.45 


388 


.50 


389 


.55 


390 


.60 


391 


.64 


392 


.69 


393 


.74 


36 


394 


.79 


395 


.84 


396 


.88 


397 


.93 


398 


.98 


400 


.03 


401 


.07 


402 


.12 


403 


.17 


404 , 


.21 


37 


405 


.26 


406 


.31 


407 


.36 


408 


.40 


409 


.45 


410 


.50 


411 


.54 


412 


.59 


413 


.64 


414 , 


.68 


38 


415 


.73 


416 


.78 


417 


.82 


418 


.87 


419 


.92 


420 


.96 


422 


.01 


423 


.05 


424 


.10 


425 , 


.15 


39 


426 


.19 


427 


.24 


428 


.28 


429 


.33 


430 


.38 


431 , 


.42 


432 


.47 


433 


.51 


434 , 


.56 


435 , 


.60 


40 


436 


.65 


437 


.69 


438 


.74 


439 


.78 


440 


.83 


441 , 


.88 


442 


.92 


443 


.96 


445 , 


.01 


446 , 


.05 


41 


447 


.10 


448 


.14 


449 


.19 


450 


.23 


451 


.28 


452 


.32 


453 


.37 


454 , 


.41 


455 


.46 


456 , 


.50 


42 


457 


.54 


458 


.59 


459 


.63 


460 


.68 


461 


.72 


462 


.77 


463 


.81 


464 , 


.85 


465 , 


.90 


466 , 


.94 


43 


467 


.98 


469 


.03 


470 


.07 


471 


.12 


472 


.16 


473 


.20 


474 


.25 


475 


.29 


476 , 


.33 


477 , 


.38 


44 


478 


.42 


479 


.46 


480 


.51 


481 


.55 


482 


.59 


483 


.63 


484 


.68 


485 , 


.72 


486 , 


.76 


487 . 


.81 


45 


488 


.85 


489 


.89 


490 


.93 


491 


.98 


493 


.02 


494 


.06 


495 


.10 


496 


.15 


497 , 


.19 


498 , 


.23 


46 


499 


.27 


500 


.32 


501 


.36 


502 


.40 


503 


.44 


504 , 


.49 


505 


.53 


506 , 


.57 


507 , 


.61 


508 . 


,65 


47 


509 


.69 


510 


.74 


511 


.78 


512 


.82 


513 


.86 


514 


.90 


515 


.95 


516 


.99 


518 , 


.03 


519 , 


.07 


48 


520 


.11 


521 


.15 


522 


.19 


523 


.23 


524 


.28 


525 


.32 


526 


.36 


527 


.40 


528 , 


.44 


529 , 


.48 


49 


530 


.52 


531 


.56 


532 


.60 


533 


.65 


534 


.69 


535 


.73 


536 


.77 


537 


.81 


538 


.85 


539 , 


.89 
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Prob[ Xdf < table value ] = 1 - a =0.8 



df 




0 




1 




2 




3 




4 




5 




6 




7 




8 




9 


0 






1 . 


64 


3 . 


22 


4 . 


64 


5 . 


99 


7 . 


29 


8 . 


56 


9 . 


80 


11 . 


03 


12 . 


24 


1 


13 . 


44 


14 . 


63 


15 . 


81 


16 . 


98 


18 . 


15 


19 . 


31 


20 . 


47 


21 . 


61 


22 . 


76 


23 . 


90 


2 


25 . 


04 


26 . 


17 


27 . 


30 


28 . 


43 


29 . 


55 


30 . 


68 


31 . 


79 


32 . 


91 


34 . 


03 


35 . 


14 


3 


36 . 


25 


37 . 


36 


38 . 


47 


39 . 


57 


40 . 


68 


41 . 


78 


42 . 


88 


43 . 


98 


45 . 


08 


46 . 


17 


4 


47 . 


27 


48 . 


36 


49 . 


46 


50 . 


55 


51 . 


64 


52 . 


73 


53 . 


82 


54 . 


91 


55 . 


99 


57 . 


08 


5 


58 . 


16 


59 . 


25 


60 . 


33 


61 . 


41 


62 . 


50 


63 . 


58 


64 . 


66 


65 . 


74 


66 . 


82 


67 . 


89 


6 


68 . 


97 


70 . 


05 


71 . 


13 


72 . 


20 


73 . 


28 


74 . 


35 


75 . 


42 


76 . 


50 


77 . 


57 


78 . 


64 


7 


79 . 


71 


80 . 


79 


81 . 


86 


82 . 


93 


84 . 


00 


85 . 


07 


86 . 


13 


87 . 


20 


88 . 


27 


89 . 


34 


8 


90 . 


41 


91 . 


47 


92 . 


54 


93 . 


60 


94 . 


67 


95 . 


73 


96 . 


80 


97 . 


86 


98 . 


93 


99 . 


99 


9 


101 . 


05 


102 . 


12 


103 . 


18 


104 . 


24 


105 . 


30 


106 . 


36 


107 . 


43 


108 . 


49 


109 . 


55 


110 . 


61 


10 


111 . 


67 


112 . 


73 


113 . 


79 


114 . 


84 


115 . 


90 


116 . 


96 


118 . 


02 


119 . 


08 


120 . 


14 


121 . 


19 


11 


122 . 


25 


123 . 


31 


124 . 


36 


125 . 


42 


126 . 


48 


127 . 


53 


128 . 


59 


129 . 


64 


130 . 


70 


131 . 


75 


12 


132 . 


81 


133 . 


86 


134 . 


91 


135 . 


97 


137 . 


02 


138 . 


08 


139 . 


13 


140 . 


18 


141 . 


24 


142 . 


29 


13 


143 . 


34 


144 . 


39 


145 . 


44 


146 . 


50 


147 . 


55 


148 . 


60 


149 . 


65 


150 . 


70 


151 . 


75 


152 . 


80 


14 


153 . 


85 


154 . 


90 


155 . 


95 


157 . 


00 


158 . 


05 


159 . 


10 


160 . 


15 


161 . 


20 


162 . 


25 


163 . 


30 


15 


164 . 


35 


165 . 


40 


166 . 


45 


167 . 


49 


168 . 


54 


169 . 


59 


170 . 


64 


171 . 


69 


172 . 


73 


173 . 


78 


16 


174 . 


83 


175 . 


88 


176 . 


92 


177 . 


97 


179 . 


02 


180 . 


06 


181 . 


11 


182 . 


15 


183 . 


,20 


184 . 


25 


17 


185 . 


29 


186 . 


,34 


187 . 


38 


188 , 


,43 


189 . 


47 


190 . 


52 


191 . 


56 


192 . 


61 


193 . 


,65 


194 . 


70 


18 


195 . 


74 


196 . 


,79 


197 . 


83 


198 . 


,88 


199 . 


92 


200 . 


96 


202 . 


01 


203 . 


05 


204 . 


,10 


205 . 


,14 


19 


206 . 


18 


207 . 


23 


208 . 


27 


209 . 


,31 


210 . 


,35 


211 . 


,40 


212 . 


44 


213 . 


,48 


214 . 


.52 


215 . 


,57 


20 


216 . 


,61 


217 . 


.65 


218 . 


69 


219 . 


,73 


220 . 


,78 


221 . 


,82 


222 . 


,86 


223 . 


,90 


224 , 


.94 


225 . 


,98 


21 


227 . 


03 


228 . 


,07 


229 . 


,11 


230 . 


,15 


231 . 


,19 


232 . 


23 


233 , 


.27 


234 . 


,31 


235 , 


.35 


236 , 


.39 


22 


237 . 


43 


238 . 


.47 


239 . 


,51 


240 . 


,55 


241 . 


,59 


242 . 


.63 


243 , 


.67 


244 . 


.71 


245 , 


.75 


246 . 


.79 


23 


247 . 


,83 


248 . 


.87 


249 . 


.91 


250 . 


.95 


251 , 


.99 


253 . 


.02 


254 . 


,06 


255 . 


.10 


256 , 


.14 


257 , 


.18 


24 


258 . 


22 


259 , 


.26 


260 , 


.29 


261 . 


.33 


262 , 


.37 


263 . 


.41 


264 . 


.45 


265 , 


.49 


266 


.52 


267 , 


.56 


25 


268 . 


.60 


269 , 


.64 


270 . 


.67 


271 , 


.71 


272 , 


.75 


273 , 


.79 


274 , 


.82 


275 


.86 


276 


.90 


277 , 


.93 


26 


278 . 


.97 


280 , 


.01 


281 . 


.05 


282 


.08 


283 


.12 


284 , 


.16 


285 . 


.19 


286 , 


.23 


287 


.27 


288 


.30 


27 


289 . 


.34 


290 


.37 


291 


.41 


292 


.45 


293 


.48 


294 


.52 


295 


.55 


296 , 


.59 


297 


.63 


298 


.66 


28 


299 , 


.70 


300 


.73 


301 , 


.77 


302 , 


.80 


303 


.84 


304 


.87 


305 


.91 


306 


.94 


307 


.98 


309 


.02 


29 


310 , 


.05 


311 , 


.09 


312 


.12 


313 


.15 


314 


.19 


315 


.22 


316 


.26 


317 , 


.29 


318 


.33 


319 


.36 


30 


320 , 


.40 


321 


.43 


322 


.47 


323 


.50 


324 


.53 


325 


.57 


326 


.60 


327 


.64 


328 


.67 


329 


.70 


31 


330 


.74 


331 


.77 


332 


.81 


333 


.84 


334 


.87 


335 


.91 


336 


.94 


337 


.97 


339 


.01 


340 


.04 


32 


341 


.07 


342 


.11 


343 


.14 


344 


.17 


345 


.21 


346 


.24 


347 


.27 


348 


.31 


349 


.34 


350 


.37 


33 


351 


.40 


352 


.44 


353 


.47 


354 


.50 


355 


.54 


356 


.57 


357 


.60 


358 


.63 


359 


.67 


360 


.70 


34 


361 


.73 


362 


.76 


363 


.79 


364 


.83 


365 


.86 


366 


.89 


367 


.92 


368 


.95 


369 


.99 


371 


.02 


35 


372 


.05 


373 


.08 


374 


.11 


375 


.15 


376 


.18 


377 


.21 


378 


.24 


379 


.27 


380 


.30 


381 


.34 


36 


382 


.37 


383 


.40 


384 


.43 


385 


.46 


386 


.49 


387 


.52 


388 


.55 


389 


.59 


390 


.62 


391 


.65 


37 


392 


.68 


393 


.71 


394 


.74 


395 


.77 


396 


.80 


397 


.83 


398 


.86 


399 


.89 


400 


.93 


401 


.96 


38 


402 


.99 


404 


.02 


405 


.05 


406 


.08 


407 


.11 


408 


.14 


409 


.17 


410 


.20 


411 


.23 


412 


.26 


39 


413 


.29 


414 


.32 


415 


.35 


416 


.38 


417 


.41 


418 


.44 


419 


.47 


420 


.50 


421 


.53 


422 


.56 


40 


423 


.59 


424 


.62 


425 


.65 


426 


.68 


427 


.71 


428 


.74 


429 


.77 


430 


.80 


431 


.83 


432 


.86 


41 


433 


.89 


434 


.91 


435 


.94 


436 


.97 


438 


.00 


439 


.03 


440 


.06 


441 


.09 


442 


.12 


443 


.15 


42 


444 


.18 


445 


.21 


446 


.24 


447 


.26 


448 


.29 


449 


.32 


450 


.35 


451 


.38 


452 


.41 


453 


.44 


43 


454 


.47 


455 


.50 


456 


.52 


457 


.55 


458 


.58 


459 


.61 


460 


.64 


461 , 


.67 


462 


.70 


463 , 


.72 


44 


464 


.75 


465 


.78 


466 


.81 


467 


.84 


468 


.87 


469 


.89 


470 


.92 


471 , 


.95 


472 


.98 


474 , 


.01 


45 


475 


.03 


476 


.06 


477 


.09 


478 


.12 


479 


.15 


480 


.17 


481 , 


.20 


482 , 


.23 


483 


.26 


484 , 


.29 


46 


485 


.31 


486 


.34 


487 


.37 


488 


.40 


489 


.42 


490 


.45 


491 , 


.48 


492 , 


.51 


493 


.54 


494 


.56 


47 


495 


.59 


496 


.62 


497 


.64 


498 


.67 


499 , 


.70 


500 


.73 


501 , 


.75 


502 , 


.78 


503 


.81 


504 , 


.84 


48 


505 


.86 


506 


.89 


507 


.92 


508 


.94 


509 


.97 


511 


.00 


512 


.03 


513 , 


.05 


514 


.08 


515 


.11 


49 


516 


.13 


517 


.16 


518 


.19 


519 


.21 


520 


.24 


521 


.27 


522 


.29 


523 , 


.32 


524 


.35 


525 , 


.37 
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APPENDIX E : Tabulated Run Results for RETP1 
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Table 1A : 



8 Exp in Series, RS = 0.931 (Hi) 

min X = 0.0002 f/hr, max X = 0.0016 f/hr, UT = 10 hrs 



S/N 


Test 

Plan 


Deg of 
Freedom 


a 


Measures of Accuracy 


RSLOW 


LEVEL 


1 


Test 5 until 
5 failed. 

NFC=40 


2*NFC 

(80) 


0.1 


0.919 


0.982 


0.2 


0.919 


0.960 


2*(NFC+ 

NCOMP) 

(96) 


0.1 


0.906 


1.000 


0.2 


0.905 


0.999 


2*NFC- 

NCOMP 

(72) 


0.1 


0.927 


0.949 


0.2 


0.927 


0.880 


2* (NFC- 
NCOMP) 
(64) 


0.1 


0.934 


0.821 


0.2 


0.930 


0.702 


2 


Test 15 until 
15 failed. 

NFC =120 


2*NFC 

(240) 


0.1 


0.928 


0.955 


0.2 


0.927 


0.908 


2*(NFC+ 

NCOMP) 

(256) 


0.1 


0.923 


0.999 


0.2 


0.923 


0.975 


2*NFC- 

NCOMP 

(232) 


0.1 


0.930 


0.916 


0.2 


0.934 


0.833 


2*(NFC- 

NCOMP) 

(224) 


0.1 


0.932 


0.844 


0.2 


0.932 


0.747 


, 3 


Test 15 until 
11 failed. 

NFC=88 


2*NFC 

(176) 


0.1 


0.927 


0.955 


0.2 


0.926 


0.916 


2*(NFC+ 

NCOMP) 

(192) 


0.1 


0.921 


0.9i6 


0.2 


0.920 


0.988 


2*NFC- 

NCOMP 

(168) 


0.1 


0.930 


0.916 


0.2 


0.927 


0.88V 


2*(NFC- 

NCOMP) 

(160) 


0.1 


0.933 


0.843 


0.2 


0.932 


0.735 
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Table 1A : 



8 Exp in Series, RS = 0.931 (Hi) (Cont...) 

min X = 0.0002 f/hr, max X = 0.0016 f/hr, UT = 10 hrs 



S/N 


Test 

Plan 


Deg of 
Freedom 


a 


Measures of Accuracy 


RSLOW 


LEVEL 


4 


Test 15 until 
7 failed. 

NFC =56 


2*NFC 

(112) 


0.1 


0.923 


0.970 


0.2 


0.923 


0.931 


2*(NFC+ 

NCOMP) 

(128) 


0.1 


0.915 


0.998 


0.2 


0.913 


0.994 


2*NFC- 

NCOMP 

(104) 


0.1 


0.929 


0.919 


0.2 


0.929 


0.853 


2*(NFC- 

NCOMP) 

(96) 


0.1 


0.934 


0.835 


0.2 


0.933 


0.720 


5 


Test 15 until 
3 failed. 

NFC =24 


2*NFC 

(48) 


0.1 


0.915 


0.986 


0.2 


0.912 


0.975 


2*(NFC+ 

NCOMP) 

(64) 


0.1 


0.891 


1.000 


0.2 


0.888 


1.000 


2*NFC- 

NCOMP 

(40) 


0.1 


0.927 


0.944 


0.2 


0.926 


0.860 


2* (NFC- 
NCOMP) 
(32) 


0.1 


0.939 


0.753 


0.2 


0.939 


0.634 
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Table IB : 8 Exp in Series, RS = 0.803 (Lo) 

min X = 0.001 f/hr, max X = 0.0045 f/hr, UT = 10 hrs 



S/N 


Test 


Deg of 




Measures of Accuracy 


Plan 


Freedom 


a 


RSLOW 


LEVEL 


1 


Test 5 until 


2*NFC 


0.1 


0.773 


0.986 




5 failed. 


(80) 


0.2 


0.793 


0.963 




NFC=40 


2*(NFC+ 

NCOMP) 

(96) 


0.1 


0.738 


1.000 






0.2 


0.736 


0.999 






2*NFC- 

NCOMP 

(72) 


0.1 


0.79* 


0.85* 






0.2 


0.792 


0.892 






2* (NFC- 
NCOMP) 
(64) 


0.1 


0.811 


0.831 






0.2 


0.812 


0.713 


2 


Test 15 until 


2*NFC 


0.1 


0.794 


0.962 




15 failed. 


(240) 


0.2 


0.793 


0.916 




NFC =120 


2*(NFC+ 

NCOMP) 

(256) 


0.1 


0.783 


0.993 






0.2 


0.79* 


0.981 






2’" NFC- 
NCOMP 
(232) 


0.1 


0.800 


0.924 






0.2 


0.799 


0.840 






2*(NFC- 

NCOMP) 

(224) 


0.1 


0.806 


0.858 










0.2 


0.805 


0.755 


3 


Test 15 until 


2*NFC 


0.1 


0.792 


0.986 




11 failed. 


(176) 


0.2 


0.790 


0.921 




NFC=88 


2*(NFC+ 

NCOMP) 

(192) 


0.2 


0.776 


0.996 






0.2 


0.774 


0.981 






2*NFC- 

NCOMP 

(168) 


0.1 


0.800 


0.922 






0.2 


0.798 


0.855 






2*(NFC- 


0.1 


0.806 


0.855 






NCOMP) 

(160) 


0.2 


0.806 


0.746 1 
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Table IB : 



8 Exp in Series, RS = 0.803 (Lo) (Cont...) 

min A = 0.001 f/hr, max A = 0.0045 f/hr, UT = 10 hrs 



S/N 


Test 

Plan 


Deg of 
Freedom 


a 


Measures of Accuracy 


RSLOW 


LEVEL 


4 


Test 15 until 
7 failed. 

NFC =56 


2*NFC 

(112) 


0.1 


0.785 


0.974 | 


0.2 


0.782 


0.938 


2*(NFC+ 

NCOMP) 

(128) 


0.1 


0.700 


0.998 


0.2 


0.756 


0.996 


2*NFC- 

NCOMP 

(104) 


0.1 


0.798 


0.925 


0.2 


0.795 


0.860 


2*(NFC- 

NCOMP) 

(96) 


0.1 


0.811 


0.841 






0.2 


0.808 


0.727 


5 


Test 15 until 
3 failed. 

NFC =24 


2*NFC 

(48) 


0.1 


0.759 


0.977 


0.2 


0.755 


0.977 


2*(NFC+ 

NCOMP) 

(64) 


0.1 


0.700 


1.000 


0.2 


0.793 


1.000 




2*NFC- 

NCOMP 

(40) 


0.1 


0.791 


0.949 


0.2 


0.789 


0.872 


2*(NFC- 

NCOMP) 

( 32 ) 


0.1 


0.825 


0.763 


0.2 


0.825 


0.642 



98 



Table 2A : 



8 Wei in Series, RS = 0.980 (Hi) 

min \ = 0.001 f/hr, max \ - 0.008 f/hr, UT = 10 hrs 



S/N 


Test 

Plan 


Deg of 
Freedom 


a 


Measures of Accuracy 


RSLOW 


LEVEL 


1 


Test 5 until 
5 failed. 

NFC =40 


2*NFC 

(80) 


0.1 


0.977 


0.992 


0.2 


0.930 


0.989 


2*(NFC + 
NCOMP) 
(96) 


0.1 


0.937 


0.994 


0.2 


0.943 


0.993 


2*NFC- 

NCOMP 

(72) 


0.1 


0.951 


0.989 


0.2 


0.937 


0.8#G 


2*(NFC- 

NCOMP) 

(64) 


0.1 


0.956 


0.985 


0.2 


0.943 


0.981 


2 


Test 15 until 
15 failed. 

NFC =120 


2*NFC 

(240) 


0.1 


0.97# 


0.918 


0.2 


0.974 


0.913 


2*(NFC+ 

NCOMP) 

(256) 


0.1 


0.977 


0.931 


0.2 


0.972 


0.924 


2*NFC- 

NCOMP 

(232) 


0.1 


0.979 


0.914 


0.2 


0.975 


0.901 


2*(NFC- 

NCOMP) 

(224) 


0.1 


0.9#§ 


0.904 


0.2 


0.975 


0.86£ 


3 


Test 15 until 
11 failed. 

NFC=88 


2*NFC 

(176) 


0.1 


0.982 


0.876 


0.2 


0.977 


0.860 


2*(NFC+ 

NCOMP) 

(192) 


0.1 


0.9?C 


0.9S.4 


0.2 


0.979 


0.882 


2*NFC- 

NCOMP 

(168) 


0.1 


0.983 


0.861 


0.2 


0.978 


0.839 


2*(NFC- 

NCOMP) 

(160) 


0.1 


0.97* 


0.840 


0.2 


0.979 


0.819 
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Table 2A : 



8 Wei in Series, RS = 0.980 (Hi) (Cont...) 

min X = 0.001 f/hr, max X = 0.008 f/hr, UT = 10 hrs 



S/N 


Test 

Plan 


Deg of 
Freedom 


a 


Measures of Accuracy 


RSLOW 


LEVEL 


4 


Test 15 until 
7 failed. 

NFC =56 


2*NFC 

(112) 


0.1 


0.987 


0.800 


0.2 




D.981 


0.779 


2*(NFC+ 

NCOMP) 

(128) 


0.1 


0.985 


0.839 


0.2 


0.978 


0.824 


2*NFC- 

NCOMP 

(104) 


0.1 


0.988 


0.776 


0.2 


0.982 


0.753 


2* (NFC- 
NCOMP) 
(96) 


0.1 


0.989 


0.746 


0.2 


0.989 


0.732 


5 


Test 15 until 
3 failed. 

NFC =24 


2*NFC 

(48) 


0.1 


0.991 


0.621 


0.2 


0.991 


0.584 


2*(NFC+ 

NCOMP) 

(64) 


0.1 


0.993 


0.705 


0.2 


0.988 


0.685 


2*NFC- 

NCOMP 

(40) 


0.1 


0.995 


0.548 


0.2 


0.992 


0.514 


2*(NFC- 

NCOMP) 

(32) 


0.1 


0.996 


0.468 


0.2 


0.993 


0.417 
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Table 2B : 



8 Wei in Series, RS = 0.832 (Lo) 

min \ = 0.003 f/hr, max A = 0.024 f/hr, UT = 10 hrs 



S/N 


Test 

Plan 


Deg of 
Freedom 


a 


Measures of Accuracy 


RSLOW 


LEVEL 


1 


Test 5 until 
5 failed. 

NFC=40 


2'" NFC 
(80) 


0.1 


0.736 


0.985 


0.2 


0.688 


0.983 


2*(NFC+ 

NCOMP) 

(96) 


0.1 


0.696 


0.992 


0.2 


0.641 


0.991 


2*NFC- 

NCOMP 

(72) 


0.1 


0.757 


0.980 


0.2 


0.713 


0.973 


2* (NFC- 
NCOMP) 
(64) 


0.1 


0.778 


0.968 


0.2 


0.739 


0.954 


2 


Test 15 until 
15 failed. 

NFC =120 


2*NFC 

(240) 


0.1 


0.834 


0.876 


0.2 


0.812 


0.876 


2*(NFC + 
NCOMP) 
(256) 


0.1 


0.825 


0.920 


0.2 


0.801 


0.968 


2*NFC- 

NCOMP 

(232) 


0.1 


0.837 


0.882 


0.2 


0.817 


0.858 


2*(NFC- 

NCOMP) 

(224) 


0.1 


0.844 


0.866 


0.2 


0.823 


0.838 


1 3 


Test 15 until 
11 failed. 

NFC=88 


2*NFC 

(176) 


0.1 


0.854 


0.838 


0.2 


0.831 


0.8S8 


2*(NFC+ 

NCOMP) 

(192) 


0.1 


0.842 


0.882 


0.2 


0.817 


0.861 


2*NFC- 

NCOMP 

(168) 


0.1 


0.860 


0.809 


0.2 


0.837 


0.776 


2*(NFC- 

NCOMP) 

(160) 


0.1 


0.865 


0.777 


0.2 


0.844 


0.734 
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Table 2B : 



8 Wei in Series, RS = 0.832 (Lo) (Cont...) 

min X = 0.003 f/hr, max X = 0.024 f/hr, UT = 10 hrs 



S/N 


Test 

Plan 


Deg of 
Freedom 


a 


Measures of Accuracy 


RSLOW 


LEVEL 


4 


Test 15 until 
7 failed. 

NFC =56 


2*NFC 

(112) 


0.1 


0.878 


0.738 




0.2 


-0.857 


i 

0.708 


2*(NFC+ 

NCOMP) 

(128) 


0.1 


0.864 


0.804 


0.2 


0.839 


0.780 


2*NFC- 

NCOMP 

(104) 


0.1 


0.886 


0.702 


0.2 


0.864 


0.655 


2*(NFC- 

NCOMP) 

(96) 


0.1 


0.894 


0.642 






0.2 


0.875 


0.587 


5 


Test 15 until 
3 failed. 

NFC =24 


2*NFC 

(48) 


0.1 


0.910 


0.611 


0.2 


0.888 


A can 

0.560 


2*(NFC+ 

NCOMP) 

(64) 


0.1 


0.885 


0.752 


0.2 


0.856 


0.719 


2*NFC- 

NCOMP 

(40) 


0.1 


0.923 


0.515 


0.2 


0.905 


0.462 


2*(NFC- 

NCOMP) 

(32) 


0.1 


0.936 


0.390 


0.2 


0.922 


0.306 || 
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Table 3A : 



4 Exp and 4 Wei (Mixed) in Series, RS = 0.980 (Hi) 
min A = 0.002 f/hr, max A - 0.008 i/hr, UT = 10 hrs 



S/N 


Test 

Plan 


Deg of 
Freedom 


a 


Measures of Accuracy 


RSLOW 


LEVEL 


1 


Test 5 until 
5 failed. 

NFC=40 


2*NFC 

(80) 


0.1 


0.979 


0.942 


0.2 


0.978 


0.905 


2*(NFC+ 

NCOMP) 

(96) 


0.1 


0.975 


0.987 


0.2 


0.9*1 


0.976 


2*NFC- 

NCOMP 

(72) 


0.1 


0.981 


0.881 


0.2 


0.980 


0.69* 


2*(NFC- 

NCOMP) 

(64) 


0.1 


0.98# 


0.771 


0.2 


0.982 


0.684 


2 


Test 15 until 
15 failed. 

NFC =120 


2*NFC 

(240) 


0.1 


0.981 


0.863 


0.2 


0.980 


0.800 


2*(NFC+ 

NCOMP) 

(256) 


0.1 


0.9*1 


0.771 


0.2 


0.9*9 


0.898 


2*NFC- 

NCOMP 

(232) 


0.1 


0.981 


0.881 


0.2 


0.980 


0.805 


2*(NFC- 

NCOMP) 

(224) 


0.1 


0.982 


0.725 


0.2 


0.981 


0.631 


3 


Test 15 until 
11 failed. 

NFC=88 


2*NFC 

(176) 


0.1 


0.981 


0.864 


0.2 


0.981 


0.801 


2*(NFC+ 

NCOMP) 

(192) 


0.1 


0.979 


0.951 


0.2 


0.978 


0.907 


2*NFC- 

NCOMP 

(168) 


0.1 


0.982 


0.69* 


0.2 


0.98* 


0.698 


2*(NFC- 

NCOMP) 

(160) 


0.1 


0.982 


0.702 


0.2 


0.982 


0.591 
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Table 3A : 



4 Exp and 4 Wei (Mixed) in Series, RS = 0.980 (Hi) (Cont...) 
min \ = 0.002 f/hr, max A = 0.008 f/hr, UT = 10 hrs 



S/N 


Test 

Plan 


Deg of 
Freedom 


a 


Measures of Accuracy 


RSLOW 


LEVEL 


4 


Test 15 until 
7 failed. 

NFC =56 


2*NFC 

(112) 


0.1 


0.981 


0.865 


0.2 


0.980 


0.787 


2*(NFC+ 

NCOMP) 

(128) 


0.1 


0.978 


0.952 


0.1 


0.978 


0.920 


2*NFC- 

NCOMP 

(104) 


0.1 


0.982 


0.769 


0.2 


0.975 


0.676 


2* (NFC- 
NCOMP) 
(96) 


0.1 


0.983 


0.644 






0.2 


0.983 


0.523 


5 


Test 15 until 
3 failed. 

NFC =24 


2*NFC 

(48) 


0.1 


0.981 


0.843 


0.2 


0.981 


0.762 


2*(NFC+ 

NCOMP) 

(64) 


0.1 


0.976 


0.970 


0.2 


0.975 


0.941 


2*NFC- 

NCOMP 

(40) 


0.1 


0.984 


0.684 


0.2 


0.984 


0.580 


2*(NFC- 

NCOMP) 

(32) 


0.1 


0.987 


0.459 


0.2 


0.987 


0.356 
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Table 3B : 



4 Exp and 4 Wei (Mixed) in Series, RS = 0.809 (Lo) 
min A = 0.002 f/hr, max \ = 0.008 f/hr, UT = 10 hrs 



S/N 


Test 

Plan 


Deg of 
Freedom 


a 


Measures of Accuracy 


RSLOW 


LEVEL 


1 


Test 5 until 
5 failed. 

NFC=40 


2*NFC 

(80) 


0.1 


0.788 


0.961 


0.2 


0.777 


0.930 


2*(NFC+ 

NCOMP) 

(96) 


0.1 


0.754 


0.996 


0.2 


0.794 


0.987 


2*NFC- 

NCOMP 

(72) 


0.1 


0.805 


0.910 


0.2 


0.796 


0.874 


2* (NFC- 
NCOMP) 
(64) 


0.1 


0.823 


0.838 


0.2 


0.815 


0.741 


2 


Test 15 until 
15 failed. 

NFC =120 


2*NFC 

(240) 


0.1 


0.808 


0.909 


0.2 


0.805 


0.842 


2*(NFC + 
NCOMP) 
(256) 


0.1 


0.797 


0.962 


0.2 


0.794 


0.917 


2*NFC- 

NCOMP 

(232) 


0.1 


0.814 


0.854 


0.2 


0.819 


0.778 


2*(NFC- 

NCOMP) 

(224) 


0.1 


0.819 


0.787 


0.2 


0.817 


0.7*7 


3 


Test 15 until 
11 failed. 

NFC =88 


2*NFC 

(176) 


0.1 


0.811 


0.885 


0.2 


0.807 


0.820 


2*(NFC+ 

NCOMP) 

(192) 


0.1 


0.797 


0.962 


0.2 


0.792 


0.925 


2*NFC- 

NCOMP 

(168) 


0.1 


0.811 


0.821 


0.2 


0.814 


0.9*7 


2*(NFC- 

NCOMP) 

(160) 


0.1 


0.826 


0.741 


0.2 


0.822 


0.647 
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Table 3B : 



4 Exp and 4 Wei (Mixed) in Series, RS = 0.809 (Lo) (Cont...) 
min \ = 0.002 f/hr, max A = 0.008 f/hr, UT = 10 hrs 



S/N 


Test 

Plan 


Deg of 
Freedom 


a 


Measures of Accuracy 


RSLOW 


LEVEL 


4 


Test 15 until 
7 failed. 

NFC =56 


2*NFC 

(112) 


0.1 


0.814 


0.872 


0.2 


.0.810 


0.792 


2*(NFC+ 

NCOMP) 

(128) 


0.1 


0.792 


0.963 


0.2 


0.787 


0.931 


2*NFC- 

NCOMP 

(104) 


0.1 


0.825 


0.775 


0.2 


0.822 


0.685 


2* (NFC- 
NCOMP) 
(96) 


0.1 


0.836 


0.663 


0.2 


0.834 


0.550 


5 


Test 15 until 
3 failed. 

NFC = 24 


2*NFC 

(48) 


0.1 


0.825 


0.836 


0.2 


0.815 


0.755 


2*(NFC+ 

NCOMP) 

(64) 


0.1 


0.780 


0.970 


0.2 


0.766 


0.940 


2*NFC- 

NCOMP 

(40) 


0.1 


0.849 


0.970 


0.2 




0.842 


0.553 


2*(NFC- 

NCOMP) 

(32) 


0.1 


0.874 


0.437 


0.2 


0.869 


0.338 
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APPENDIX F : Tabulated Run Results for RETP2 
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Table 4A : 8 Exp in Series, RS = 0.961 (Hi) 

X = 0.001 f/hr, UT = 5 hrs 



S/N 


Degrees 


K / E[NFC] 




Measures of Accuracy 


of 

Freedom 


(TT) 


a 


RSLOW 


LEVEL 


1 


2*(1+NFC) 


0.25 / 1.2 


0.1 


0.950 


0.851 






(225) 


0.1 


.0.935 


0.851 






0.5 / 2.4 


0.1 


0.957 


0.857 






(450) 


0.2 


0QS4 


0.857 






1.0 / 4.8 


0.1 


0.957 


0.941 






(900) 


0.2 


0.957 


0.850 






2.0 / 9.6 


0.1 


0.958 


0.916 






(1800) 


0.2 


0.960 


' 

0.850 






3.0 / 14.4 
(2700) 


0.1 


0.959 


0.916 






0.2 


0.959 


0.809 






4.0 / 19.2 
(3600) 


0.1 


0.959 


0.937 






0.2 


0.960 


0.843 






5.0 / 24 


0.1 


0.960 


0.826 






(4500) 


0.2 


0.961 


0.914 






10.0 / 48 


0.1 


0.960 


0.924 






(9000) 


0.2 


0.960 


0.809 






20.0 / 96 


0.1 


0.960 


0.914 






(18000) 


0.2 


0.961 


0.820 






30.0 / 144 


0.1 


0.961 


0.906 






(27000) 


0.2 


0.961 


0.804 
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Table 4A : 8 Exp in Series, RS = 0.961 (Hi) (Cont...) 

\ = 0.001 f/hr, UT = 5 hrs 



S/N 


Degrees 


K / E[NFC] 




Measures of Accuracy 


of 

Freedom 


(TT) 


a 


RSLOW 


LEVEL 


2 


1.3* 


0.25 / 1.2 


0.1 


0.950 


0.851 




2*(1 + NFC) 


(225) 


0.2 


0.935 


0.851 






0.5 / 2.4 


0.1 


0.95# 


0.857 






(450) 


0.2 


0.941 


0.857 






1.0 / 4.8 


0.1 


0.946 


0.8S1 






(900) 


0.2 


0.935 


0.941 






2.0 / 9.6 


0.1 


0.949 


0.989 






(1800) 


0.2 


0.950 


0.948 






3.0 / 14.4 
(2700) 


0.1 


0.949 


0.99* 






0.2 


0.949 


0.969 






4.0 / 19.2 
(3600) 


0.1 


0.948 


0.998 






0.2 


0.950 


0.987 






5.0 / 24 


0.1 


0.949 


0.998 






(4500) 


0.2 


0.949 


0.985 






10.0 / 48 


0.1 


0.949 


1.000 






(9000) 


0.2 


0.949 


0.997 






20.0 / 96 


0.1 


0.949 


1.000 






(18000) 


0.2 


0.950 


1.000 






30.0 / 144 


0.1 


0.950 


1.000 






(27000) 


0.2 


0.949 


1.000 



109 



Table 4B : 8 Exp in Series, RS = 0.819 (Lo) 

X = 0.005 f/hr, UT = 5 hrs 



S/N 


Degrees 

of 

Freedom 


K / E[NFC] 




Measures of Accuracy 


(TT) 


a 


RSLOW 


LEVEL 


1 


2*(1 + NFC) 


0.25 / 1.2 


0.1 


0.774 


0.851 






(45) 


0.2 


.0.702 


0.851 






0.5 / 2.4 


0.1 


0.801 


0.857 






(90) 


0.2 


0.788 


0.857 






1.0 / 4.8 


0.1 


0.801 


0.919 






(180) 


0.2 


0.803 


0.850 






2.0 / 9.6 


0.1 


0.807 


0.916 






(360) 


0.2 


0.814 


0.847 






3.0 / 14.4 
(540) 


0.1 


0.812 


0.916 






0.2 


0.812 


0.809 






4.0 / 19.2 
(720) 


0.1 


0.809 


0.923 






0.2 


0.817 


0.840 






5.0 / 24 


0.1 


0.813 


0.925 






(900) 


0.2 


0.816 


0.814 






10.0 / 48 


0.1 


0.816 


0.919 






(1800) 


0.2 


0.817 


0.809 






20.0 / 96 


0.1 


0.817 


0.814 






(3600) 


0.2 


0.818 


0.820 






30.0 / 144 


0.1 


0.817 


0.907 






(5400) 


0.2 


0.818 


0.804 | 
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Table 4B : 



8 Exp in Series, RS = 0.819 (Lo) (Cont...) 
X = 0.005 f/hr, UT = 5 hrs 



S/N 


Degrees 

of 

Freedom 


K / E[NFC] 
(TT) 


Q 


Measures of Accuracy 


RSLOW 


LEVEL 


2 


1.3* 

2*(1+NFC) 


0.25 / 1.2 
(45) 


0.1 


0.774 


0.851 


0.2 


0.702 


0.851 


0.5 / 2.4 
(90) 


0.1 


0.801 


0.857 


0.2 


0.736 


0.857 


1.0 / 4.8 
(180) 


0.1 


0.759 


0.981 


0.2 


0.753 


0.941 


2.0 / 9.6 
(360) 


0.1 


0.762 


0.989 


0.2 


0.771 


0.948 


3.0 / 14.4 
(540) 


0.1 


0.770 


0.991 


0.2 


0.768 


0.969 


4.0 / 19.2 
(720) 


0.1 


0.766 


0.997 


0.2 


0.772 


0.969 


5.0 / 24 
(900) 


0.1 


0.769 


0.998 


0.2 


0.771 


0.985 




10.0 / 48 
(1800) 


0.1 


0.771 


1.000 


0.2 


0.771 


0.997 


20.0 / 96 
(3600) 


0.1 


0.771 


1.000 


0.2 


0.772 


1.000 


30.0 / 144 
(5400) 


0.1 


0.773 


1.000 


0.2 


0.772 


1.000 
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Table 5A : 8 Wei in Series, RS = 0.956 (Hi) (*) 

X = 0.005 f/hr, UT = 15 hrs 



S/N 


Degrees 

of 

Freedom 


K / E[NFC] 
(TT) 


a 


Measures of Accuracy 


RSLOW 


LEVEL 


1 


2*(1 + NFC) 


0.25 / 1.2 
(45) 


0.1 


1.000 


0.186 


0.2 


1.000 


0.158 


0.5 / 2.4 
(90) 


0.1 


0.986 


0.501 


0.2 


0.979 


0.458 


1.0 / 4.8 
(180) 


0.1 


0.967 


0.767 


0.2 


0.960 


0.732 


2.0 / 9.6 
(360) 


0.1 


0.957 


0.879 




0.2 


0.952 


0.854 


3.0 / 14.4 
(540) 


0.1 


0.952 


0.934 


0.2 


0.946 


0.922 


4.0 / 19.2 
(720) 


0.1 


0.952 


0.934 


0.2 


0.946 


0.928 


5.0/24 

(900) 


0.1 


0.952 


0.940 


0.2 


0.946 


0.928 


10.0 / 48 
(1800) 


0.1 


0.952 


0.940 


0.2 


0.946 


0.928 


20.0 / 96 
(3600) 


0.1 


0.952 


0.940 


0.2 


0.946 


0.928 


30.0 / 144 
(5400) 


0.1 


0.952 


0.940 


0.2 


0.946 


0.928 



(*) 20 test items for each Weibull component. 
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Table 5A : 



8 Wei in Series, RS = 0.956 (Hi) (*) (Cont...) 
X = 0.005 f/hr, UT = 15 hrs 



S/N 


Degrees 

of 

Freedom 


K / E[NFC] 
(TT) 


a 


Measures of Accuracy 


RSLOW 


LEVEL 


2 


1.3* 

2*(1 + NFC) 


0.25 / 1.2 
(45) 


0.1 


1.000 


0.258 


0.2 


0.999 


0.224 


0.5 / 2.4 
(90) 


0.1 


0.983 


0.635 


0.2 


0.913 


0.593 


190 / 4.8 
(180) 


0.1 


0.958 


0.884 


0.2 


0.930 


0.866 


2.0 / 9.6 
(360) 


0.1 


0.946 


0.963 


0.2 


0.939 


0.956 


3.0 / 14.4 
(540) 


0.1 


0.939 


0.984 


0.2 


0.930 


0.976 


4.0 / 19.2 
(720) 


0.1 


0.938 


0.987 


0.2 


0.930 


0.981 


5.0 / 24 
(900) 


0.1 


0.938 


0.987 


0.2 


0.930 


0.981 


10.0 / 48 
(1800) 


0.1 


0.938 


0.987 


0.2 


0.946 


0.981 


20.0 / 96 
(3600) 


0.1 


0.938 


0.987 


0.2 


0.930 


0.981 


30.0 / 144 
(5400) 


0.1 


0.938 


0.987 


0.2 


0.930 


0.981 
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Table 5B : 8 Wei in Series, RS = 0.835 (Lo) (*) 

A = 0.01 f/hr, UT = 15 hrs 



S/N 


Degrees 


K / E[NFC] 




Measures of Accuracy 


of 

Freedom 


(TT) 


a 


RSLOW 


LEVEL 


1 


2*(1 + NFC) 


0.25 / 1.2 


0.1 


0.966 


0.222 I 


; 




(22.5) 


0.2 


.0.963 




0.161 | 

] 






0.5 / 2.4 


0.1 


0.932 


0367 






(45) 


0.2 


0.913 


0.319 






1.0 / 4.8 


0.1 


0.874 


0.704 






(90) 


0.2 


0.858 




0.658 | 






2.0 / 9.6 


0.1 


0.842 


0.851 






(180) 


0.2 


0.832 


0.819 






3.0 / 14.4 
(270) 


0.1 


0.829 


0.924 










0.901 








0.2 


0.814 






4.0 / 19.2 
(360) 


0.1 


0.827 


0.924 






0.2 


0.813 




0.908 






5.0 / 24 


0.1 


0.827 


0.928 






(450) 


0.2 


0.813 


0.908 






10.0 / 48 


0.1 


0.827 


0.928 






(900) 


0.2 


0.813 




0.908 






20.0 / 96 


0.1 


0.827 


0.928 






(1800) 


0.2 


0.813 


0.908 






30.0 / 144 


0.1 


0.827 


0.928 






(2700) 


0.2 


0.813 


0.908 | 
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Table 5B : 



8 Wei in Series, RS = 0.835 (Lo) (*) 
\ = 0.01 f/hr, UT = 15 hrs 



S/N 


Degrees 

of 

Freedom 


K / E[NFC] 
(TT) 


a 


Measures of Accuracy 


RSLOW 


LEVEL 


2 


1.3* 

2*(1+NFC) 


0.25 / 1.2 
(22.5) 


0.1 


0.957 


0.326 


0.2 


0.953 


0.268 


03 / 2.4 
(45) 


0.1 


0.914 


0.572 


0.2 


0.890 


0.522 


1.0 / 4.8 
(90) 


0.1 


0.842 


0.881 


0.2 


0.821 


0.862 


2.0 / 9.6 
(180) 


0.1 


0.802 


0.968 


0.2 


o.im 


0.962 


3.0 / 14.4 
(270) 


0.1 


0.786 


0.986 


0.2 


0.766 


0.983 


4.0 / 19.2 
(360) 


0.1 


0.784 


0.989 


0.2 


0.786 


0.986 


5.0 / 24 
(450) 


0.1 


0.784 


0.989 


0.2 


0.766 


0.986 


10.0 / 48 
(900) 


0.1 


0.784 


0.989 


0.2 


0.766 


0.986 


20.0 / 96 
(1800) 


0.1 


0.784 


0.989 


0.2 


0.766 


0.986 


30.0 / 144 
(2700) 


0.1 


0.784 


0.989 


0.2 


0.766 


0.986 
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Table 6A : 



4 Exp and 4 Wei (Mixed) in Series, RS = 0.958 (Hi) (*) 
A(exp) = 0.001 f/hr, UT(exp) = 5 hrs 
A(wei) = 0.005 f/hr, UT(wei) = 15 hrs 



S/N 


Degrees 

of 

Freedom 


K / E[NFC] 




Measures of Accuracy 


TT(exp) 

TT(wei) 


a 


RSLOW 


LEVEL 


1 


2*(1 + NFC) 


0.25 / 1.2 
(225) 
(45) 


0.1 


1.000 


0.620 






0.2 


0.995 


0.451 






0.5 / 2.4 
(450) 
(90) 


0.1 


0.969 


0.623 






0.2 


0.975 


0.573 






1.0 / 4.8 
(900) 
(180) 


0.1 


0.971 


0.736 






0.2 


0.965 


0.684 






2.0 / 9.6 
(1800) 
(360) 


0.1 


0.9S6 


0.803 






0.2 


0.960 


0.765 






3.0 / 14.4 
(2700) 
(540) 


0.1 


0.960 


0.874 






0.2 


0.956 


0.841 






4.0 / 19.2 
(3600) 
(720) 


0.1 


0.960 


0.873 






0.2 


0.959 


0.839 






5.0 / 24 
(4500) 
(900) 


0.1 


0.959 


0.887 






0.2 


0.956 


0.861 






10.0 / 48 
(9000) 
(1800) 


0.1 


0.959 


0.891 






0.2 


0.959 


0.862 






20.0 / 96 
(18000) 
(3600) 


0.1 


0.959 


0.892 






0.2 


0.955 


0.867 






30.0 / 144 
(27000) 
(5400) 


0.1 


0.9S6 


0.892 






0.2 


0.956 


0.862 
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Table 6 A : 



4 Exp and 4 Wei (Mixed) in Series, RS = 0.958 (Hi) (*) (Cont...) 
A(exp) = 0.001 f/hr, UT(exp) = 5 hrs 
A(wei) = 0.005 f/hr, UT(wei) = 15 hrs 



S/N 


Degrees 

of 

Freedom 


K / E[NFC] 
TT(exp) 
TT(wei) 


a 


Measures of Accuracy 


RSLOW 


LEVEL 


2 

1 


1.3* 

2*(1 + NFC) 


0.25 / 1.2 
(225) 
(45) 


0.1 


1.000 


0.723 


0.2 


0.994 


0.697 


0.5 / 2.4 
(450) 
(90) 


0.1 


0.977 


0.730 


0.2 


0.969 


0.681 


1.0 / 4.8 
(900) 
(180) 


0.1 


0.964 


0.853 


0.2 


0.955 


0.837 


2.0 / 9.6 
(1800) 
(360) 


0.1 


0.954 


0.934 


0.2 


0.949 


0.921 


3.0 / 14.4 
(2700) 
(540) 


0.1 


0.949 


0.961 


0.2 


0.944 


0.978 


4.0 / 19.2 
(3600) 
(720) 


0.1 


0.948 


0.975 


0.2 


0.949 


0.966 


5.0 / 24 
(4500) 
(900) 


0.1 


0.948 


0.984 


0.2 


0.944 


0.978 


10.0 / 48 
(9000) 
(1800) 


0.1 


0.948 


0.97* 


0.2 


0.943 


0.965 


20.0 / 96 
(18000) 
(3600) 


0.1 


0.949 


0.9?* 


0.2 


0.942 


0.979 


30.0 / 144 
(27000) 
(5400) 


0.1 


0.949 


0.978 


0.2 


0.944 


0.971 
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Table 6B : 



4 Exp and 4 Wei (Mixed) in Series, RS = 0.827 (Lo) (*) 
A(exp) = 0.005 f/hr, UT(exp) = 5 hrs 
A(wei) = 0.010 f/hr, UT(wei) = 15 hrs 



S/N 


Degrees 

of 

Freedom 


K / E[NFC] 




Measures of Accuracy 


TT(exp) 

TT(wei) 


Q 


RSLOW 


LEVEL 


1 


2*(1 + NFC) 


0.25 / 1.2 
(45) 
(22.5) 


0.1 


0.971 


0.691 






0.2 


0.938 


0.685 






0.5 / 2.4 
(90) 
(45) 


0.1 


0.916 


0.600 






0.2 


0.893 


0.532 






1.0 / 4.8 
(180) 
(90) 


0.1 


0.876 


0.684 






0.2 


0.855 


• 

0.638 






2.0 / 9.6 
(360) 
(180) 


0.1 


0.849 


0.787 






0.2 


0.836 


0.750 






3.0 / 14.4 
(540) 
(270) 


0.1 


0.838 


0.855 






0.2 


0.824 


0.817 






4.0 / 19.2 
(720) 
(360) 


0.1 


0.83S 


0.870 






0.2 


0.829 


0.819 






5.0 / 24 
(900) 
(450) 


0.1 


0.829 


0.890 






0.2 


0.822 


0.838 






10.0 / 48 
(1800) 
(900) 


0.1 


0.827 


0.900 






0.2 


0.820 


0.863 






20.0 / 96 
(3600) 
(1800) 


0.1 


0.826 


0.902 






0.2 


0.819 


0.871 






30.0 / 144 


0.1 


0.831 


0.684 






(5400) 

(2700) 


0.2 


0.821 


0.854 | 
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Table 6B : 



4 Exp and 4 Wei (Mixed) in Series, RS = 0.827 (Lo) (*) (Cont...) 
A(exp) = 0.005 f/hr, UT(exp) = 5 hrs 
X(wei) = 0.010 f/hr, UT(wei) = 15 hrs 



S/N 


Degrees 

of 

Freedom 


K / E[NFC] 
TT(exp) 
TT(wei) 


a 


Measures of Accuracy 


RSLOW 


LEVEL 


2 


1.3* 

2*(1 + NFC) 


0.25 / 1.2 
(45) 
(22.5) 


0.1 


0.964 


0.718 


0.2 


0.922 


0.700 


0.5 / 2.4 
(90) 
(45) 


0.1 


0.895 


0.9*8 


0.2 


0.866 


0.688 


1.0 / 4.8 
(180) 
(90) 


0.1 


0.846 


0.864 


0.2 


0.818 


0.9*8 


2.0 / 9.6 
(360) 
(180) 


0.1 


0.811 


0.939 


0.2 


0.79* 


0.921 


3.0 / 14.4 
(540) 
(270) 


0.1 


0.798 


0.95* 


0.2 


0.780 


0.952 


4.0 / 19.2 
(720) 
(360) 


0.1 


0.772 


0.9S9 


0.2 


0.788 


0.972 


5.0 / 24 
(900) 
(450) 


0.1 


0.787 


0.985 


0.2 


0.777 


0.98* 


10.0 / 48 
(1800) 
(900) 


0.1 


0.78* 


0.688 


0.2 


0.774 


0.983 


20.0 / 96 
(3600) 
(1800) 


0.1 


0.783 


0.993 


0.2 


0.772 


0.988 


30.0 / 144 
(5400) 
(2700) 


0.1 


0.788 


0.996 


0.2 


0.775 


0.994 
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